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SUMMARY 


Dynamically  scaled  model  helicopter  rotor  tests  were  carried  out  to  collect 
data  for  the  verification  and  improvement  of  helicopter  rotor  and  blade  re¬ 
sponse  calculations. 

Four  fully  articulated  model  rotor  blade  configurations  were  fabricated  and 
tested.  The  configurations  included  variation  in  built-in  twist ,  airfoil 
camoer,  and  stiffness.  The  dynamically  scaled  configurations  were  based  on 
H-3^  rotor  blades  previously  tested  at  full  scale.  The  geometry,  mass,  and 
stiffness  properties  of  the  model  blades  are  provided,  along  with  two- 
dimensional  aerodynamic  test  data  for  the  blade  airfoil  sections.  The  rotor 
operating  regime  tested  for  the  four  configurations  included  lightly  loaded 
and  stalled  conditions,  variations  in  flapping  trim,  and  rotor-wake  inter¬ 
ference  conditions. 

The  complete  body  of  rotor  performance  and  blade  response  time  history  data 
is  preserved  in  digital  form  on  magnetic  tapes.  A  computer  program  to 
sample  the  data  tapes  and  perform  the  most  common  types  of  processing  is 
also  available  for  convenient  future  use. 

Comp  arise".  8  between  model  and  full-scale  performance  and  blade  response 
data  show  qualitative  and  quantitative  similarities.  It  is  therefore  con¬ 
cluded  that  full-scale  rotor  performance  and  blade  response  calculation 
methods  can  be  verified  or  improved  by  correlating  with  model  data. 

Sample  model  test  data  are  compared  with  blade  response  and  rotor  per¬ 
formance  calculations.  The  model  rotor  provided  more  lift  than  could  be 
expected  from  calculations  employing  static  two-dimensional  airfoil  test 
data  and  the  consideration  of  full-scele  unsteady  aerodynamics.  The 
accuracy  of  the  blade  response  and  performance  calculation  method  varied 
between  the  blade  configurations.  Good  accuracy  was  obtained  with  the 
conventional  amorjvt  of  built-in  twist. 

Correlation  of  model  results  obtained  herein  with  previous  full-scale  test 
data  should  be  extended  for  further  evaluation  of  the  relationship  between 
model  and  full-scale  results. 

Rotor  performance  and  blade  response  calculations  should  generally  include 
the  effects  of  wake-induced  inflow.  The  analytical  representation  of  rotor 
blade  excitation  by  wake  vortices  should  be  improved. 

The  analytical  representation  and  data  for  considering  unsteady  aerodynamics 
in  rotor  performance  and  blade  reeponse  calculations  should  be  augmented  to 
Include  characteristics  of  various  symmetrical  and  unsyumetrical  airfoils, 
other  Reynolds  number  variations,  and  any  appreciable  effects  of  yawed 
flow. 
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PREFACE 


This  program  was  conducted  for  the  Eustis  Directorate,  U.  S.  Army  Air 
Mobility  Research  and  Development  Laboratory,  under  Contract  DAAJ02-72-C- 
0026  (Pro.lect  1F263211D157) .  Technical  monitors  for  the  Government  were 
Mr.  Patrick  Cancro  and  Mr.  Edward  Austin. 

Mr.  Edmond  F.  Kiely  was  Task  Manager  for  this  work  at  Sikorsky  Aircraft. 

He  also  supervised  the  design,  fabrication,  and  pretest  checkout  of  model 
components.  Mr.  Charles  F.  Niebanck  supervised  the  analytical  calculations, 
performed  their  correlation  with  test  data,  and  provided  specifications  for 
data  reduction  programs.  Mr.  Robert  Murrill  served  as  the  Sikorsky  test 
engineer  during  the  pretest  activities  and  the  wind  tunnel  operations. 

Mr.  Miles  Waugh  wrote  the  computer  programs  which  stored  the  test  data  on 
tape,  and  which  permit  their  retrieval  for  future  use.  Mr.  John  Corrigan 
used  the  Y-200  version  of  the  Normal  Modes  Aeroelastic  Analysis  to  provide 
the  calculated  blade  response  data.  Technical  supervision  was  provided  by 
Dr.  Raymond  Carlson,  Mr.  Peter  Arcidiacono,  and  Mr.  Edward  S.  Carter. 

This  volume  covers  the  testing  and  correlation  with  calculation  for  four 
model  blades  based  on  the  H-31*  rotor.  Under  this  program,  similar  work 
vas  carried  out  for  the  more  advanced  Sikorsky  IRB  blade  design.  The  data 
and  results  pertaining  to  this  configuration  are  contained  in  a  separate 
report  and  are  provided  with  limited  rights  for  the  use  of  the  Government 
only. 
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INTRODUCTION 


Experimental  data  are  required  to  verify  and  suggest  means  of  improving 
mathematical  analyses  of  helicopter  rotor  performance  and  blade  response. 
The  dynamically  scaled  model  rotor  is  an  excellent  source  of  such  data. 
Dynamically  scaled  model  wind  tunnel  data  for  the  verification  of  analyti¬ 
cal  calculations  can  be  obtained  rapidly  for  an  extended  range  of  accurate¬ 
ly  measured  parameters  at  comparatively  low  cost  and  risk. 

Dynamically  scaled  rotor  models  may  also  play  an  important  direct  role  as 
a  design  tool  for  specific  rotor  systems.  Application  of  dynamic  models 
in  this  way  requires  some  comparative  model  and  full-scale  data  for  an 
existing  configuration,  to  evaluate  the  effects  of  parameters  which  are  not 
simulated  in  the  dynamic  scale  modeling  process.  A  body  of  data  of  this 
type  allows  more  confident  extrapolation  of  the  characteristics  of  new 
full-scale  rotor  designs  from  model  test  results. 

The  primary  objective  of  the  dynamically  scaled  model  rotor  test  documented 
herein  was  the  collection  of  data  for  the  verification  and  improvement  of 
helicopter  rotor  performance  and  blade  response  calculations.  The  four 
model  rotor  blade  configurations  fabricated  and  tested  represented  inde¬ 
pendent  variations  in  twist,  camber  and  stiffness.  The  configurations  fab¬ 
ricated  and  tested  included  dynamically  scaled  fiber  glass  spar  blade  sets 
with  0  degrees  and  -8  degrees  of  built-in  twist,  a  similar  set  with  0 
degrees  twist  and  a  20J 5  chord,  5-degree  deflection  plain  flap,  and  a  fourth 
set  with  a  three-times-scale  stiffness  aluminum  spar.  The  dynamically 
scaled  configurations  were  based  on  the  H-31*  blades  tested  at  full  scale  in 
the  programs  documented  by  Reference  1,  2,  and  3.  Thus,  the  present  data 
may  also  be  useful  for  further  establishing  the  relationship  between  full- 
scale  and  dynamically  scaled  rotor  data. 

A  secondary  objective  of  the  program  was  the  correlation  of  sample  test 
data  with  calculated  results  from  the  Sikorsky-UARL  Normal  Modes  Aeroelastic 
Analysis,  to  support  an  assessment  and  recommendations  for  improvement  of 
that  analysis. 

Test  data  resulting  from  this  program  are  appropriately  sampled  in  this  re¬ 
port  for  presentation  and  correlation  with  calculated  and  full-scale  re¬ 
sults.  The  complete  body  of  rotor  performance  and  blade  response  time  his¬ 
tory  data  are  preserved  in  digital  form  on  magnetic  tapes.  A  computer  pro¬ 
gram  to  sample  the  data  tapes  and  perform  the  most  usual  processing  is  also 
supplied  for  convenient  future  use.  The  magnetic  tapes  are  readable  by 
using  standard  computer  program  routines.  Additional  processing  modes  which 
may  be  devised  for  analysis  and  correlation  of  the  test  data  may  be  applied 
in  the  future  by  providing  appropriate  computer  program  modules. 
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MODEL  DEFINITION 


Data  which  describe  the  various  H-31*  model  blade  configurations,  the  cor¬ 
responding  rotor  system,  and  the  wind  tunnel  model  supporting  the  rotor  are 
supplied  herein. 

These  data  are  required  for  input  to  rotor  response  calculation  methods 
whose  outputs  will  be  compared  to  model  test  results. 

The  data  are  supplemented  with  descriptions  of  how  the  data  were  obtained. 
Where  applicable,  discussions  of  data  verification  and  estimates  of  the 
data  accuracy  are  also  supplied. 

GENERAL  DESCRIPTION  OF  MODEL  AND  BLADE  CONFIGURATIONS  TESTED 

The  wind  tunnel  model  with  H-3b  blades  installed  is  shown  mounted  in  the 
18-foot  test  section  of  the  United  Aircraft  Research  Laboratories  main 
wind  tunnel  in  Figures  1  and  2.  The  compound  helicopter  model  fuselage 
utilized  for  this  test  was  identical  to  that  used  in  the  TRAC  Telescoping 
Rotor  test  program  (Reference  M.  A  six-component  strain  gage  balance  of 
appropriate  capacity  was  installed.  The  110-inch-diameter  four-bladed 
fully  articulated  rotor  is  powered  by  a  variable  frequency  electric  motor. 
The  control  system  provides  independent  collective  pitch,  longitudinal 
cyclic  pitch,  and  lateral  cyclic  pitch.  A  slip  ring  installed  below  the 
swashplate  transfers  the  blade  strain  gage  signals  and  flap  and  lag  po¬ 
tentiometer  signals  to  the  nonrotating  system.  The  strain  gage  balance  for 
measuring  rotor  forces  and  moments  is  mounted  below  the  rotor  transmission 
and  supports  the  entire  drive  train  and  control  system  as  well  as  the  rotor. 
The  grounded  portion  of  the  balance  is  supported  in  a  soft  damped  gimbal 
frame  which  has  been  found  effective  in  eliminating  ground  resonance  type 
mechanical  instabilities. 

The  model  fuselage  consists  of  a  nonstructural  fiber  glass  and  aluminum 
outer  shell  mounted  on  a  steel  frame.  Overall  fuselage  length  is  123 
inches  and  maximum  diameter  is  17.3  inches.  The  general  arrangement  is 
shown  in  Figure  3,  where  the  side  panels  and  nose  cone  have  been  removed 
from  the  model. 

The  four-bladed  fully  articulated  model  rotor  head  is  of  conventional  de¬ 
sign.  The  coincident  flap  and  lag  hinges  are  located  at  the  3-inch  radial 
station  corresponding  to  an  offset  ratio  of  .  05^*5  •  Rotary  viscous  dampers 
equipped  with  3:1  ratio  mechanical  linkages  are  provided  for  each  lag 
hinge.  A  pitch  flap  coupling  ratio  of  0.0  was  used.  Pitch-lag  coupling 
was  measured  to  be  approximately  0.8°  for  lag  motions  from  0°  to  18°  lag. 
Physical  stops  limited  flapping  motion  from  -7.5  to  +21.0°  and  Inplane 
motions  from  7°  lead  to  18°  lag. 

Blade  pitch  control  was  achieved  with  an  electromechanical  swashplate  con¬ 
trol  system.  The  system  used  three  independent  electric  motors  controlled 
from  the  rotor  control  console  to  provide  pure  collective  and  cyclic  con¬ 
trol  inputs  to  the  swashplate.  Required  mixing  is  accomplished 


mechanically  with  three  intermeshed  differential  gearsets.  Mixed  outputs 
are  transferred  to  three  conventionally  positioned  mechanical  actuators  by 
means  of  flexible  shafts.  The  actuators  transfer  the  rotary  motion  of  the 
shafts  into  linear  motion  by  means  of  a  worm  gear  and  Jackscrew.  Control 
angles  are  measured  with  geared  potentiometers  which  effectively  count  the 
number  of  turns  of  the  input  electric  motors.  The  impressed  collective  and 
cyclic  pitch  angles  are  displayed  on  three  digital  meters  in  the  control 
console . 

The  rotor  system  is  powered  by  a  water-cooled  variable  frequency  electric 
mctor,  nominally  rated  at  80  horsepower.  Power  is  transmitted  to  the  rotor 
shaft  by  a  5.25:1  reduction,  two-stage  gearbox.  Gearbox  lubrication  and 
cooling  are  provided  by  an  oil  spray  system. 

The  blades,  rotor  head,  motor,  gearbox,  blade  pitch  control  system,  and 
instrumentation  slipring  comprise  the  internal  model  assembly.  This  por¬ 
tion  of  the  model  is  metric.  It  is  mounted  on  a  six-component  Task  Mk.  XIX 
strain  gage  balance  for  precision  measurement  of  rotor  forces  and  moments. 
The  ground  side  of  the  balance  is  supported  in  a  soft  damped  gimbal  mount. 
The  gimbal  provides  a  pitch  and  roll  freedom  for  the  internal  model  about 
the  axis  formed  by  the  intersection  of  the  rotor  shaft  and  balance  center- 
lines.  Spring  restraint  is  provided  between  the  gimbal  (ground  side  of 
the  balance)  and  the  fuselage.  Cantilever  springs  are  used  in  the  pitch 
direction  and  a  torsion  bar  is  used  in  the  roll  direction.  Damping  is 
provided  by  rotary  viscous  dampers  with  mechanical  linkages  to  increase 
their  effective  output  for  small  gimbal  motions.  One  set  of  dampers  was 
connected  directly  across  the  balance  in  the  roll  direction.  It  was 
demonstrated  experimentally  that  this  soft  grounding  of  the  balance  did 
not  affect  the  steady-state  force  and  moment  measurements.  The  gimbal 
mounting  permitted  the  internal  model  assembly  to  respond  to  destabilizing 
forces  by  rotating  at  the  gimbal  axis  rather  than  storing  the  energy  as 
elastic  deformations  in  the  metric  structure.  The  dampers  dissipate  the 
energy  as  the  gimbal  deflects,  resulting  in  a  stable  system. 

The  wind  tunnel  model  uv  unting  strut  system  shown  in  Figures  1  and  2  was 
developed  in  conjunction  with  the  gimbal  mounting  of  the  rotor  to  provide 
as  nearly  as  possible  a  rigid  support  for  the  external  fuselage.  The  stiff 
mount  was  effective  in  raising  the  primary  natural  frequencies  of  the 
model  installation  to  values  above  those  at  which  mechanical  instabilities 
such  as  ground  resonance  were  likely  to  occur. 

Three  of  the  model  operating  parameters  (forward  speed,  rotor  rpm,  and 
fuselage  angle  of  attack)  were  controlled  from  the  wind  tunnel  facility 
control  console.  The  remainder  of  the  model  controls  were  centralized  in 
the  model  control  console.  The  collective  pitch  and  longitudinal  and 
lateral  cyclic  pitch  settings  were  independently  controlled  by  three-posi¬ 
tion  spring-centered  toggle  switches.  The  impressed  blade  pitch  angles 
were  displayed  on  digital  meters  immediately  above  the  corresponding  toggle 
switches. 

The  four  model  rotor  blade  configurations  fabricated  and  tested  represented 
independent  variations  in  twist,  camber  and  stiffness.  The  configurations 
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fabricated  and  tested  included  dynamically  scaled  fiber  glass  spar  blade 
sets  with  0  degrees  and  -8  degrees  of  built-in  twist,  a  similar  set  with 
0  degrees  twist  and  a  20%  chord,  5-degree  deflection  plain  flap,  and  a 
fourth  set  with  a  three-times -scale  stiffness  aluminum  spar.  The  dynami¬ 
cally  scaled  configurations  were  based  on  the  H-3U  blades  tested  at  full 
scale  in  the  programs  documented  by  References  1,  2,  and  3. 

BLADE  CONSTRUCTION 

An  exploded  view  shoving  typical  construction  details  for  the  H-3^*  dynami¬ 
cally  scaled  blades  is  shown  in  Figure  U. 

The  internal  load-carrying  members  of  the  blade  are  representative  of  full- 
scale  construction,  with  a  load-carrying  D-tube  spar  at  the  leading  edge. 
The  trailing  half  of  the  blade  is  made  up  of  individually  attached  fairing 
segments  or  pockets.  The  pockets  are  fastened  to  the  spar,  but  not  to  each 
other . 

BLADE  STIFFSESS  AND  MASS  DATA 
Presentation  of  Blade  Data 

Inertia  and  chordwise  center-of-gravity  data  for  the  H-3^  model  blade  con¬ 
figurations  are  presented  in  Table  I.  The  blade  stiffness  properties  for 
the  scale  stiffness  fiber  glass  and  the  three-times-scale  stiffness  alumi¬ 
num  spar  configurations  of  the  H-3U  model  blades  are  given  in  Tables  II 
and  III  respectively.  The  actual  distributions  are  represented  by  short 
constant  property  radial  segments  over  the  intervals  given. 

Derivation  of  Blade  Data 


The  blade  distributed  mass  data  outboard  of  radial  station  8.90  arc  deter¬ 
mined  by  dividing  the  full-scale  blade  data  by  the  appropriate  power  of 
the  geometric  scale  factor,  S8  which  is  6.109  in  this  case.  Thus,  v  is 
found  by  dividing  corresponding  full-scale  data  by  s£  and  the  moments  of 
inertia  by  dividing  full-scale  data  by  sjj. 

Blade  mass  properties  Inboard  of  radial  station  8.90  are  not  scaled,  and 
represent  the  properties  of  the  actual  components  provided  under  the  scope 
of  this  investigation. 

Blade  area  moment  of  inertia  and  torsional  constant  data  for  the  fiber  glass 
blade  properties  outboard  of  radial  station  8.90  are  determined  by  dividing 
full-scale  blade  spar  properties  by  S8.  These  values  were  then  adjusted 
in  accordance  vith  stiffness  test  results  for  the  completed  blade.  When 
these  values  are  used  vith  the  values  of  E  and  G  given  below  the  table, 
the  resulting  bending  and  twisting  stiffness  values  EIfl  ,  EIchd,  and  GJ 
for  the  fiber  glass  blades  are  one-fourth  of  those  for  a*nypothetical  model- 
scale  replica  of  the  full-scale  blade  spar. 

Although  the  values  of  Ifl  ,  Ichd»  and  J  based  on  the  actual  crossectional 
area  of  the  fiber  glass-resin  composite  are  not  known  precisely,  the  values 
derived  as  described  above  can  be  considered  representative,  because  of 
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the  thin-walled  D-tube  spar  and  pocket  construction,  which  is  similar  to 
the  actual  blade. 

blade  area  moment  of  inertia  and  torsional  constant  data  for  the  aluminum 
spar  blade  properties  outboard  of  radial  station  8.90  were  based  on  values 
obtained  by  dividing  full-scale  blade  spar  values  by  the  quantity  Usi^/3. 
These  values  were  then  increased  in  accordance  with  stiffness  test  results 
for  the  completed  blade,  including  the  pockets.  When  these  values  are  used 
with  the  values  of  E  and  G  given,  the  resulting  bending  and  twisting  stiff¬ 
ness  values  Elfip,  EIchd,  and  GJ  obtained  are  approximately  three-fourths 
of  those  for  a  hypothetical  xodel-scale  replica  of  the  full-scale  blade. 

Inboard  of  station  8.90,  the  bending  and  twisting  stiffness  distributions 
are  based  on  reasonable  approximations  of  the  actual  cross-sections.  The 
cross-sectional  area  contributions  of  the  various  components  were  weighted 
by  the  ratio  of  their  elastic  or  shear  moduli  (E  or  G)  to  the  reference 
values  given  in  Tables  II  and  III. 

Verification  of  Blade  Data 

Blade  bending  stiffness  is  verified  by  applying  a  known  weight  to  the  blade 
tip  at  the  feathering  axis  and  measuring  the  resulting  deflection  at  twelve 
locations.  This  test  also  verifies  that  the  elastic  axis  and  feathering 
axis  are  coincident  within  the  accuracy  of  measurements.  A  photograph  of 
the  apparatus  used  in  this  test  is  provided  in  Figure  5. 

Blade  twisting  rtiffness  is  verified  by  applying  a  known  torsional  couple 
at  the  blade  tip  and  optically  measuring  the  torsional  deflections  at  six 
spanwlse  stations. 

Results  of  a  set  of  such  measurements  for  the  three-times-scale  stiffness 
aluminum  spar  blade  and  the  scale  stiffness  fiber  glass  spar  blades  are 
shown  in  Figures  6,  T,  8  and  9.  The  test  data  are  compared  with  nominal 
expected  deflections  which  were  calculated  from  beam  theory  for  blade  spar 
properties  alone.  Most  of  the  discrepancy  between  the  calculated  and  ex¬ 
perimental  values  is  believed  to  be  due  to  stiffness  contributions  of  the 
pockets.  This  deviation  is  considered  acceptable. 

Blade  mass  is  verified  by  careful  weighing  of  all  significant  components 
before  assembly.  Upon  completion,  the  total  blade  weight  and  spanwlse 
balance  are  checked. 

EXPERIMENTAL  DETERMINATION  OF  NONROTATING  MODEL  BLADE  NATURAL  FREQUENCIES 

The  nonrotating  natural  frequencies  of  each  blade  set  were  determined  ex¬ 
perimentally.  The  frequencies  of  both  instrumented  blades  from  each  set 
were  checked  as  were  several  of  the  noninstrumented  blades.  The  blades 
were  supported  from  an  articulated  model  rotor  hub  which  was  attached  to 
a  rigid  structure.  A  variable  frequency  air  Jet  was  used  to  excite  the 
blade  elastic  motions.  Strain  gage  outputs  from  the  instrumented  blade 
were  recorded  on  a  direct  writing  oscillograph  for  ease  in  determining 
maximum  response,  mode  shape,  and  for  measuring  decay  rate  for 
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determination  of  structural  damping  coefficients.  The  natural  frequencies 
of  the  noninstrumented  blades  were  determined  by  visual  observations  of 
maximum  blade  excursions. 

The  test  apparatus  utilized  for  these  determinations  is  shown  in  Figure  10. 
The  rotor  blades  were  mounted  from  sleeve-spindle  assemblies  identical  to 
those  used  during  the  wind  tunnel  tests.  The  assembly  was  supported  ver¬ 
tically  from  a  model  rotor  hub.  The  rotor  hub  bearings  provided  a  pin 
connection  at  the  coincident  flap  and  lag  hinge  axis.  The  lag  dampers 
were  disconnected,  and  blade  pitch  freedcm  was  eliminated  by  grounding  the 
pitch  horn  to  the  flapping  pivot. 

The  variable  frequency,  variable  amplitude  air  jet  was  obtained  by  connect¬ 
ing  a  function  generator  to  a  hydraulic  servo  controller.  The  associated 
servo  valve  was  used  to  drive  a  shuttle  valve  vh.cn  provided  modulation  to 
a  "shop  air"  supply  line.  Excitation  frequency  vas  limited  to  100  Hz  where 
compressibility  effects  eliminated  the  modulation  at  the  Jet  nozzle. 

The  experimentally  determined  nonrotating  natural  frequencies  and  critical 
damping  coefficients  derived  from  the  strain  gage  decay  records  are  listed 
in  Table  IV.  The  experimentally  determined  natural  frequencies  are  com¬ 
pared  to  calculated  values  in  Figures  11  and  12.  Reasonable  agreement  is 
further  verification  that  mass  and  stiffness  distributions  are  as  planned. 
Furthermore,  lack  of  significant  coupling  between  flapwise  and  torsional 
moments  in  the  natural  vibration  test  verifies  the  coincidence  of  the  blade 
elastic  axis  and  chordwise  center-of-gravity  position. 

Accuracy  of  Blade  Data 

Blade  mass  and  inertia  properties  are  within  ±5Jf  of  nominal  values.  Vari¬ 
ations  between  individual  blades  are  within  t2%,  except  for  spanwise 
balance,  which  is  held  within  *.5£. 

The  various  blade  stiffness  properties  are  within  a  tolerance  of  ♦lOJ^  to 
-3%  of  design  values.  Variations  between  individual  blades  are  within  t3%. 

FULL-SCALE  EQUIVALENT  STRESS  FOR  FIBER  GLASS  MODEL  BLADES 

The  following  relationships  between  the  bending  and  twisting  moments  ani 
stresses  wtre  obtained  theoretically  for  the  constant  section  part  of  the 
full-scale  H-31*  blades: 


°F  ■  Mrss 

1.53 


°T  "  MTFS 

3.31 
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At  the  .15  radius  station,  the  relation  betveen  twisting  moment  and  stress 
was  similarly  found  to  be 


°T.15R 


gas 

5.93 


The  stress-moment  relations  for  a  hypothetical  aluminum  spar  replica  model 
are  found  by  multiplying  the  factors  in  the  above  equations  by  the  scale 
factor  cubed.  With  a  scale  factor  of  6.11,  the  above  relationships  for 
the  replica  model  cure  as  follows: 

°F  “  1149 


°T  =  69,0  “TOM 
°T.15R  "  38,5  MTRM 

The  fiber  glass  H-34  models  have  one-fourth  the  stiffness  of  a  replica 
model.  At  a  dynamically  similar  operating  condition,  the  fiber  glass  blade 
loads  will  be  one-fourth  those  of  the  replica  model  blade.  Therefore,  the 
following  relationship  is  established  between  fiber  glass  model  blade 
moments  and  full-scale  equivalent  stress,  by  multiplying  the  above  factors 
by  4: 


°FFS  "  595Mp 
°CFS  "  19>*  MC 
°TFS  "  276  **I 
°TFS.15R  "  15U  *T 

ROTOR  SYSTEM  GEOMETRY.  STIFFNESS  AND  DAMPING  DATA 


Remaining  rotor  system  data  required  for  input  to  the  computer  programs  are 
presented  in  Table  V.  Supplementary  discussion  and  estimated  accuracy  of 
these  data  are  contained  in  the  following  paragraphs. 

Geometric  and  dimensional  accuracy  of  the  blade  data  are  as  follows: 


Blade  Length 
Chord 

Lo  il  Twist 
Cumulative  Twist 
Contour 

Shear  Center  Location 
Center-of -Gravity  Location 
Location  of  Blade  Hinges 


±.020  in. 
±.018  in. 
t. 25  deg 
±.50  deg 
±.007  in. 
±.03  in. 
1.03  in. 
±.010  in. 
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The  pitch-flap  coupling  ratio  of  0.0  occurs  when  blade  pitch,  flap,  and  lag 
sure  all  zero. 

The  damper  installation  is  identical  to  that  used  for  the  Reference  5  test, 
with  the  dampers  operating  through  a  3  to  1  linkage.  The  orifices  are 
opened  wide  at  the  17.0  in. -lb-sec  setting,  minimizing  the  in-phase  or 
spring  part  of  the  damper  resistance.  Tests  of  the  dampers  at  the  expected 
range  of  once-per-revolution  frequency  and  motion  showed  that  the  coeffi¬ 
cient  is  within  ±2 Of  of  the  nominal  value,  and  that  the  apparent  damper 
spring  (in  phase  with  displacement)  moment  is  approximately  15Z  of  the 
damper  damping  (out  of  phase  with  displacement)  mouent.  This  implies  that 
the  damper  apparent  spring  constant  is  200  in. -lb/rad  or  less  at  the  once- 
per-revolution  frequency  range. 

Linear  dimensions  defining  the  pitch  control  gometry  axe  accurate  to  1.01 
in.  Angular  dimensions  are  accurate  to  ±.30  deg. 

Although  pitch  control  spring  rate  is  not  accurately  known,  and  consider¬ 
able  azimuthal  variation  in  its  value  occurs,  it  remains  high  enough  so 
that  blade  root  pi+^h  deflections  axe  nepligible.  The  21,000  in. -lb /rad 
value  given  in  Table  V  wuj  calculated  for  a  blade  pushrod  at  its  furthest 
azimuth  distance  from  one  of  the  three  actuators.  It  is  estimated  that 
blade  torsional  root  stiffness  may  differ  from  this  figure  by  ±10,000  in.- 
lb/rad.  Variations  are  caused  by  simplifying  approximations  and  the  rela¬ 
tive  azimutbwise  position  of  the  push  rods  and  actuators.  Blade  response 
observed  during  the  Reference  5  test  supports  the  conclusion  that  blade 
root  pitch  deflections  are  negligible. 

AERODYNAMIC  DATA 


Hub  Aerodynamic  Drag 

The  aerodynamic  drag  of  the  four-blad'd  hub  with  blade  cuffs  installed  was 
measured  in  a  previous  test  program,  reported  in  Reference  6.  The  measure¬ 
ments  are  expressed  as  an  equivalent  area  in  square  feet.  Multiplication 
of  this  equivalent  area  by  dynamic  pressure  in  pounds  oer  square  fooo  gives 
the  aerodynamic  drag  in  pounds. 

The  equivalent  area  varied  between  .28  square  feet  at  a  shaft  angle  of  at¬ 
tack  of  8  degrees  and  .25  square  feet  at  a  shaft  angle  of  attack  of 

-8  degrees.  The  accuracy  of  these  data  are  ±5f. 

The  hub  drag  is  subtracted  from  total  rotor  drag  during  data  reduction. 

This  provides  rotor  drag  in  terms  of  blade  contributions  only.  This  pro¬ 
cedure  is  followed  because  no  special  attempt  has  been  made  to  scale  the 
hub  configuration  within  the  scope  of  this  or  previous  programs. 

Two-Dimensional  Blade  Section  Steady  Aerodynamic  Coefficients 

The  tvo-dimensional  blade  section  steady  aerodynamic  coefficients  were  ob¬ 
tained  in  the  United  Aircraft  Research  Laboratory  Acoustic  Tunnel.  The 
test  airfoils  were  made  from  untwisted  aluminum  spar  model  H-3b  blade 
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sections.  A  steel  spar  was  Inserted  through  the  leading-edge  D-tube  and 
bonded  in  place  to  provide  reinforcement. 

Aerodynamic  load  measurements  of  lift,  drag  and  pitching  moment  were  made 
with  a  strain  gage  balance.  Tunnel  velocity  was  calculated  initially  from 
the  usual  total  and  static  pressure  measurements. 

A  photograph  of  the  wind  tunnel  and  airfoil  section  installation  appears 
in  Figure  13. 

The  data  reduction  calculations  incorporated  balance  axis  interaction,  cor¬ 
rections  for  airfoil  deflection  and  weight,  and  conventional  Jet  boundary 
corrections  from  Reference  7  for  open-top-and-bottom  wind  tunnel  test  sec¬ 
tions. 

The  test  airfoils  had  a  chord  of  2. L9  inches,  and  a  span  of  21.0  inches, 
which  extended  over  the  full  width  of  the  open-top-and-bottom  test  section. 
The  height  of  the  test  section  was  30.88  inches. 

Presentation  of  Two-Dimensional  Blade  Section  Steady  Aerodynamic 
Coefficients 


The  H-3b  model  blade  section  aerodynamic  coefficients  for  the  standard 
NACA  0012  airfoil  and  the  NACA  0012  airfoil  with  a  20X  chord  plain  trailing 
edge  flap  deflected  5  degrees  are  tabulated  in  Tables  VI  through  XV.  Data 
are  provided  for  five  Mach  numbers  ranging  from  .10  to  .56.  Angles  of 
attack  cover  ±180  degrees  at  the  two  lowest  Mach  numbers.  Angle-of-at.tack 
ranges  at  the  higher  Mach  numbers  are  adequate  to  cover  the  expected  regime 
of  operation  of  the  one-half  speed  dynamically-scaled  models. 

The  tabulated  data  are  also  plotted  in  Figures  l4  and  15. 

It  should  be  noted  that  the  data  tabulated  and  plotted  herein  have  been 
faired  from  a  larger  amount  of  test  data.  The  number  of  angle  of  attack 
and  coefficient  values  has  been  chosen  to  fit  the  Sikorsky  Normal  Modes 
Aeroelastlc  Analysis  program  input  tables.  The  data  values  tabulated  and 
plotted  herein  are  the  values  provided  to  that  computer  program  as  steady 
aerodynamic  data.  In  accordance  with  the  assumptions  used  in  that  computer 
program,  the  unflapped  airfoil  is  considered  to  be  completely  symmetrical, 
and  data  are  provided  for  positive  angles  only. 

Data  Reduction  for  Two-Dimensional  Blade  Section  Steady  Aerodynamic 
Coefficients 


This  subsection  documents  the  data  reduction  calculations  for  the  two- 
dimensional  blade  section  steady  aerodynamic  coefficients. 

The  structure  of  the  raw  data  and  the  data  reduction  calculations  them¬ 
selves  reflect  the  usual  wind  tunnel  test  procedure.  A  group  of  operating 
conditions  referred  to  as  a  run  were  preceded  and  followed  by  wind-off 
measurements  called  start  and  end  zeros,  which  provide  zero  airload 
balance  readings  to  which  wind-on  readings  are  referred. 
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The  smallest  complete  unit  of  useful  data  is,  therefore,  a  vind  tunnel  run. 
The  data  reduction  calculations  described  on  the  following  pages  were  re¬ 
peated  for  each  vind  tunnel  run. 

1.  Calculation  of  Aerodynamic  Loads  and  Geometric  Angle  of  Attack. 


The  balance  start  and  end  zeros  from  a  run  were  averaged  from  each  of  the 
three  (lift,  drag,  and  moment)  balance  outputs 


<°L6  +  GLE)/2‘° 


gdo  "  ^gds  +  gde^2*° 


gmo  *  (gms  ♦  GME>/2-° 


For  each  of  the  test  points  in  the  run,  the  difference  between  wind-on 
balance  output  and  average  wind-off  output  for  the  run  was  taken.  For  test 
point  i  the  calculations  are: 

VL(i)  *  °L(i)  "  GL0 

VD(i)  *  °D(i)  ‘  °D0 

VM(i)  *  °M(i)  '  °M0 


The  balance  outputs  were  converted  into  forces  and  moments  for  each  test 
point  through  the  following  matrix  multiplication: 


The  values  of  the  elements  of  the  C _ matrix  were  determined  experimental- 

THn 

ly.  Sample  values  are: 


Cn  ■  .0122  C12  -  .0000202  C13  «  .00023 

C21  ■  .000005  C22  -  .001*73  C23  -  .000137 

c31  ■  .000105  c32  -  .oooo»*88  c33  -  .00859 

The  Indicated  angle  cf  attack  for  each  test  point  was  determined  from  the 
transducer  output  in  a  manner  analogous  to  the  conversion  of  balance  read¬ 
ings  to  force  and  moment: 
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<0AS  *  °A£)/2-° 


VA(1)  *  °A(i)  ‘  °A0 

A(l)  “o 

KA 

The  angle  a0  Is  the  angle  at  which  static  zeros  were  taken.  This  angle 
was  0  or  180  for  most  of  the  runs. 

The  above  calculations  provide  the  aerodynamic  lift  and  drag  forces  in  the 
tunnel  axis  system,  and  the  geometric  angle  of  attack  at  the  point  where 
the  airfoil  is  attached  to  the  pitch  drive  mechanism. 

Because  of  the  flexibility  of  the  model,  its  static  unbalance,  and  the 
large  angle-of -attack  range  to  be  tested,  appropriate  corrections  to  the 
measured  pitching  moment  and  angle  of  attack  were  applied. 

The  torsional  deflection  of  the  airfoil  at  the  mid-span  was  calculated 
from  engineering  beam  theory,  and  used  as  a  correction  to  obtain  the  geo¬ 
metric  angle  of  attack  from  the  indicated  angle  of  attack. 

°G(i)  "  A(  J  +  ®(i) 

For  the  H-31*  model  blade  sections, 

0(i)  -  197  x  10”U  JIR{i)  de« 

The  indicated  lift  and  drag  forces  for  each  test  point  were  resolved  into 
normal  and  chordwise  force  components: 

"(1)  ■  hu)  CO>  °G(i)  *  DR(1)  aG(i> 

°(l)  -LR<1)  sin  *0(1)  *  DR(1)  co*  a0(i) 

The  corrections  applied  to  the  indicated  pitching  moment  included  the  fol¬ 
lowing  at  each  test  point  (i).  The  numerical  coefficients  given  were  ob¬ 
tained  from  engineering  beam  theory,  and  apply  to  the  H-31*  model  two- 
dimensional  sections: 

(a)  A  moment  due  to  the  uniformly  distributed  chordwise  force  on  an 
arm  equal  to  the  average  normal  deflection  caused  by  the  uniform¬ 
ly  distributed  normal  force: 

“nbU)  *  •°0018  C(i)N(i) 

(b)  A  moment  due  to  the  uniformly  distributed  normal  force  acting  on 
an  arm  equal  to  the  average  chordwise  deflection  caused  by  the 
uniformly  distributed  chordwise  force: 


11 


MCB(i)  "  -000379  C(iAi) 

(c)  A  moment  due  to  the  chordvise  component  of  model  veight,  acting 
through  an  arm  equal  to  the  average  normal  deflection  caused  by 
the  uniformly  distributed  normal  force: 

Vi)  -  -00818  Bu)  WH.1..0(1) 

(d)  A  moment  caused  by  airfoil  chordvise  veight  unbalance,  vhich 
varies  from  its  start  and  end  zero  value  because  of  variation  in 
angle  of  attack: 

Vi)  ■  -Vco  [“•  “o  *  cos  °G(  1 )] 

The  chordvise  unbalance  vas  slightly  nose-heavy  for  the  H-3**  model  blade 
tvo-dimensional  sections: 

WMYCG  *  °-oUl85  in. -lb 

With  this  amount  of  unbalance,  the  term  vas  negligible  vith  respect 

to  aerodynamic  pitching  moment.  ' 

The  aerodynamic  part  of  the  pitching  moment  vas  found  by  subtracting  the 
above-defined  corrections  from  the  total  pitching  moment  measured  by  the 
balance: 

MA(i)  "  Vi)  "  “UBU)  ’  M  CB(i)  '  M  NW(i )  ““uwd). 

2.  Calculation  of  Tunnel  Mach  Number.  Density.  Velocity,  and  Dynamic 
Pressure 


The  tunnel  Mach  number,  density,  velocity,  and  dynamic  pressure  vere  de¬ 
termined  initially  from  the  tunnel  total  pressure,  test  section  total  minus 
static  pressure,  and  tunnel  temperature.  The  customary  compressible,  lsen- 
tropic,  adiabatic  flov  and  gas  lav  relationships  vere  used: 


H-p  -  AP 

P  «  1  -  (3-P) 
H  H 


P.  »  H 

*0  (TSC  ♦  ^ 
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V  =  M 


*5 


*VTSC  +  U6°  } 


l  +  xzL 
2 


-1  M2 


In  evaluating  the  above  queintities,  the  following  values  for  y  and  R  were 
used: 


R_  ■  1722.  °F  -  SEC2 
0  FT2 

3.  Calculation  of  Aerodynamic  Coefficients 

The  uncorrected  aerodynamic  coefficients  were  calculated  in  the  usual  man¬ 
ner  from  the  aerodynamic  loads,  model  dimensions  and  dynamic  pressures  for 
each  test  point. 

CKi)  “  ^(i) 
q  Sc 

CdR(i)  "  PR(i) 
q  Sc 

CmR(i)  *  MA(i) 
qSc* 

In  the  above  calculations,  the  model  airfoil  dimensions  are  expressed  in 
feet,  and  the  aerodynamic  pitching  moment  in  foot-pounds. 


U.  Jet  Bound 


Corrections 


The  following  corrections  from  Reference  7  were  applied  to  the  angle  of 
attack  and  the  drag  and  moment  coefficients.  The  corrections  were  theoret¬ 
ically  derived  in  Reference  7  for  the  case  of  an  airfoil  in  an  Infinitely 
long  test  section  with  an  open  top  and  bottom  and  in  lnviscld,  incompressi¬ 
ble  flow. 


“<l)- 

Cd(i)  ■  °«R(l)  '  &  [Cl(1>] 

C«u)  *  W)  -fsrtf]  c«i> 
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5.  Alternate  Calculation  of  Tunnel  Mach  Humber 


The  two-dimensional  coefficient  of  lift  versus  angle-of -attack  relationship 
obtained  for  the  lowest  Mach  number  tested  by  using  the  test  results  and 
the  foregoing  calculations  is  compared  with  similar  data  from  References  8 
and  9  in  Figure  16. 

Figure  16  indicates  satisfactory  agreement  with  these  earlier  sources,  es¬ 
pecially  if  the  effects  of  Reynolds  number  and  Mach  number  are  considered 
in  the  comparison  with  the  Reference  9  data. 

Comparisons  at  higher  Mach  number,  which  were  made  in  the  same  manner,  were 
unsatisfactory,  and  the  problem  was  finally  attributed  to  an  error  in  the 
total  minus  static  pressure  measurements  (AP).  As  the  indicated  AP  in¬ 
creased,  the  data  reduction  calculations  gave  progressively  excessive 
values  for  the  slope  of  coefficient  of  lift  versus  a  and  smaller  than  ex¬ 
pected  values  for  maximum  C^.  In  addition,  there  was  an  excessive  indi¬ 
cated  effect  of  Mach  number  on  the  coefficients  of  drag  and  moment.  Since 
the  data  appeared  satisfactory  at  the  lowest  Mach  number  tested,  an  al¬ 
ternate  method  was  available  for  the  determination  of  test  section  Mach 
number,  from  the  observed  ratio  of  lift  to  unstalled  angle  of  attack,  the 
Prandtl-Glauert  relation,  and  the  adiabatic,  isentropic,  compressible  flow 
relationships . 

dC 

Thus,  if  the  slope  L  of  the  coefficient  of  lift  versus  a  curve  is  known 

da 

below  stall  at  a  low  Mach  number,  the  following  expression  for  lift  over  a 
can  be  written: 


-Y 


The  lift  over  angle-of -attack  ratio  vms  extracted  fror  the  partially  re¬ 
duced  data  at  each  tunnel  speed  setting,  incorporating  an  estimate  of  the 
small  Jet  boundary  correction  in  the  value  of  a  used.  Then  a  plot  of  the 
right-hand  side  of  the  above  expression  was  made  against  the  value  of  the 
single  varying  parameter  M.  (The  value  of  H  was  substantially  constant  for 
all  the  test  points.)  Finally,  the  above  expression  was  solved  graphically 
for  the  Mach  number  that  existed  at  the  tunnel  speed  setting,  and  the  cal¬ 
culations  of  velocity  and  aerodynamic  coefficients  could  proceed  as  before. 

The  above  procedure  for  recalculating  the  tunnel  Mach  number  uses  well- 
established  relationships,  and  its  accuracy  is  consistent  with  the  accuracy 
of  the  test  data  and  its  intended  use.  In  addition  to  its  effect  on  the 
coefficient  of  lift,  the  recalculation  of  velocity  also  improved  the  agree¬ 
ment  of  the  coefficient  of  drag  and  the  coefficient  of  moment  data  vith  the 
trends  to  be  expected  from  previous  data.  The  data  in  Figure  lb  for  the 
Mach  number  M  ■  .33  agree  well  with  the  data  from  Reference  9  which  is 
shown  in  Figure  16.  The  reduction  of  maximum  lift  coefficient  due  to  the 
two  higher  Mach  numbers  in  Figure  lb  is  somewhat  greater  than  the  propor¬ 
tionate  decrease  in  full-scale  data.  These  low  values  do  not,  however* 
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appreciably  affect  the  half-scale  speed  blade  response  calculations  per¬ 
formed  under  this  program.  At  the  two  highest  Mach  numbers,  angles  of 
attack  high  enough  to  reach  the  maximum  lift  coefficients  are  rarely  en¬ 
countered. 

Accuracy  and  Repeatability  of  Aerodynamic  Data 

The  accuracy  of  the  aerodynamic  coefficient  data  is  believed  to  be  governed 
chiefly  by  the  accuracy  of  the  tunnel  velocity  measurement.  The  accuracy 
is  estimated  at  5%, 

Many  of  the  test  points  were  repeated  to  give  an  indication  of  the  repeat¬ 
ability  of  the  data.  The  typical  scatter  band  of  measured  points  around 
the  faired  curves  ranges  from  tS%  at  M  *  .10  to  ±3%  at  M  «  .56  for  the  lift 
data.  The  typical  scatter  band  for  the  drag  coefficients  varies  from  1.005 
at  M  *  .10  to  1.002  at  M  *  .56.  The  scatter  band  for  the  moment  data  is 
typically  1.005  at  M  ■  .10,  and  reduces  to  about  1.001  at  M  *  .56. 

SUPPORT  SYSTEM  DATA 


Natural  Modes 

The  natural  modes  and  frequencies  of  the  support  system  were  calculated 
from  separate  considerations  of  longitudinal,  yaw-lateral,  vertical,  and 
shaft  torsional  motions  of  springs  and  masses  representing  the  system. 

The  longitudinal,  yaw-lateral,  vertical,  and  shaft  torsion  equations  con¬ 
sidered  17,  26,  5«  end  2  degrees  of  freedom  respectively. 

The  undamped  equations  of  motion  of  these  spring-mass  systems  are  of  the 
following  form: 

[A][M]  W  +  {*"}  “  0 

The  natural  frequencies  and  modes  are  then  obtained  from  the  corresponding 
eigen  problem  in  a  routine  manner. 

(M  M  MM*.}- 

Calculation  of  Effective  Mass.  Stiffness.  Damping,  ar.d  Mechanical 
Compliance  at  the  Rotor  Head  From  Natural  Modes 

In  order  to  simulate  mechanical  compliance  at  the  rotor  head  for  analytical 
purposes,  the  head  and  rotor  support  system  was  approximated  as  effective 
mass,  stiffness,  and  damping  values  in  the  lateral  and  fore-aft  directions 
in  the  plane  of  the  rotor,  and  by  similar  quantities  pertaining  to  vertical 
and  shaft  torsional  motions.  This  approximation  permits  the  correct  com¬ 
pliance  for  steady-state  sinusoidal  forcing  to  be  simulated  if  the  effec¬ 
tive  values  of  mass,  stiffness,  and  damping  are  allowed  to  vary  with  fre¬ 
quency.  In  the  present  case,  the  effective  values  were  calculated  for  the 
range  of  four  times  rotor  rotational  frequency  to  be  encountered  during  the 
model  tests.  This  is  expected  to  be  the  major  forcing  frequency  on  the 
rotor  support  and  drive  system  during  the  model  tests. 
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The  method  used  to  calculate  these  effective  values  uses  the  undamped 
natural  modes  and  frequencies  of  vibration,  normalized  to  unit  deflection 
at  the  rotor  head. 

By  applying  the  methods  of  Lagrange,  the  equations  of  motion  for  sinusoidal 
forcing  of  the  damped  linear  system  can  be  expressed  in  terms  of  the  natu¬ 
ral  vibration  modes.  For  convenience,  the  force  is  assumed  to  be  of  unit 
amplitude.  To  provide  the  desired  effective  properties  in  the  various 
directions,  the  calculation  is  repeated  with  the  loads  and  modes  pertaining 
to  those  directions. 

The  equations  of  motion  are  of  the  following  form: 


KP  MKK*8}  ♦  [yny  [k]{>}{4*}  ♦  Wtp][vm]{4m}  »{x2}cos<vr. 

Because  of  orthogonality  properties  of  the  undamped,  structural  modes,  and 
because  the  modes  are  normalized  to  have  a  unit  value  at  the  head,  where 
the  load  is  applied,  the  following  equation  is  obtained: 

[MG1  ]{■!"}  *  [  MGi  “ni  ]{ “"l  *  [Do]  {«*}  ■{I}co'V 

The  elements  of  the  generalized  damping  matrix,  Dq,  include  the  contribu¬ 
tions  of  intentionally  supplied  viscous  dampers  and  the  structural  damping 
naturally  present,  and  are  obtained  from  the  following: 

DGiJ  "  I  CDP  YDpi  YDpj  1  +  i 

DOij  -I  CDp  (YDpl)2  *  2nMGi“nl  1  *  J 
P 

where  Cpp  is  the  coefficient  of  the  pth  viscous  damper,  is  the  deflec¬ 

tion  of  that  damper  in  the  ith  mode,  and  n  is  an  equivalent  structural 
damping,  expressed  as  a  fraction  of  critical  damping. 

For  sinusoidal  forcing  at  constant  amplitude,  the  solutions  to  the  equa¬ 
tions  of  motion  are  also  sinusoidal.  In  other  vords,  the  response  of  the 
1th  mode  is  the  form 

qi  *  aico8u  pt  +  to^sinwpt 

By  equating  sine  and  cosine  terms  in  the  equation  of  forced  vibration,  the 
following  matrix  equation  for  the  a*  and  t>i  is  obtained 


M0i 

1  W- 

N 

m 

f  i  ] 

... 

—  _  — 

Vo 

I  MOi(wnl  ^ 

W 

0 

S.  > 

The  above  equation  is  solved  by  routine  methods  to  obtain  the  a^  and  b^ 
quantities . 
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The  desired  cosine  and  sine  parts  of  the  head  mechanical  compliance  are 
obtained  by  summing  the  response  of  the  individual  modes: 


^  =  Zi  ai 
JCyi  =  Z± 


An  alternative  expression  for  the  above  two  quantities  is  the  resultant 
amplitude  and  phase  angle: 

*ua  *  *ui 

*UA  *  Arot,n  <XUI/XUR> 

The  kinetic  energy  of  the  system  can  be  written  as  follows*  because  of  the 
orthogonality  of  the  modes: 


The  effective  mass  at  the  rotor  head,  to  be  used  in  the  analytical  simula¬ 
tion,  is  taken  to  be  that  mass  which  will  have  the  same  kinetic  energy  due 
to  the  response  under  a  unit  force  as  the  multimodal,  system: 

15  “eff  *ua  “f  ■  T 


m_  .ifn  <»?;  *?> 

£i  *i*  1 

The  apparent  stiffness  and  damping  are  set  at  those  values  which  will 
cause  the  effective  mass  to  have  the  same  response  to  the  vibratory  force 
as  the  multimodal  system. 


KeFF  "(Sjf’)2  +  ^ 


CEFF  “  [SrF  ~  ^EFF  “fJ 

The  above  procedure  is  applied  separately  to  obtain  values  for  longitudi¬ 
nal,  lateral,  vertical,  and  torsional  effective  mass,  stiffness,  damping, 
and  mechanical  compliance.  The  resultant  amplitude  per  unit  force,  phase 
of  the  amplitude  response,  and  effective  mass,  stiffness,  and  damping  are 
supplied  on  Figures  17  and  18. 
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Accuracy  of  Support  System  Compliance  and  Effective  Property  Data  at  the 
Rotor  Head 


Comparison  of  similar  previous  calculations  with  test  results  indicates 
variations  of  as  ouch  as  50jf  between  calculated  and  test  values.  Neverthe¬ 
less,  the  mechanical  compliance  of  the  model  rotor  support  system  is  com¬ 
paratively  low,  and  even  this  variation  is  not  expected  to  have  any  signi¬ 
ficant  effect  on  blade  stress  or  response. 
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MODEL  TEST  DATA 

WIND  TUNNEL  FACILITY 

The  United  Aircraft  Research  Laboratories  main  wind  tunnel  located  in  East 
Heir t ford,  Connecticut,  shown  in  Figure  19,  was  used  for  the  model  rotor 
tests  documented  herein.  This  tunnel  is  an  atmospheric  total  pressure, 
single-return  closed  test  section  type.  The  low-speed  test  section  which 
was  utilized  for  the  program  is  octagonal  in  shape,  approximately  18  ft 
across  the  flats,  and  is  capable  of  speeds  of  approximately  170  knots. 
Remotely  controlled  air  exchangers  located  in  the  low-velocity  portion  of 
the  circuit  are  employed  to  maintain  tunnel  stagnation  temperatures  at 
desired  values.  The  large  area  of  the  exchangers  insures  that  tunnel 
stagnation  pressure  is  atmospheric.  The  tunnel  fan  is  driven  by  a  9000- 
horsepower  electric  motor. 


Wind  tunnel  velocity  was  measured  with  static  pressure  taps  located  in  the 
side  wall  of  the  test  section.  Tunnel  stagnation  temperature  was 
measured  with  a  thermocouple  located  in  the  tunnel  settling  chamber. 

MODEL  INSTRUMENTATION  AND  CALIBRATION 


Blade  Strain  Gages 

Twelve  strain  gage  bending  bridges  were  installed  on  each  instrumented 
blade.  The  bridges  were  arranged  to  measure  flapwise  moments  at  five  radi¬ 
al  stations,  chordwise  moments  at  five  stations  and  torsional  moments  at 
two  stations.  The  scan  wise  location  and  installation  of  the  strain  gages 
are  shown  in  Figure  20. 


I 


All  blade  strain  gage  bridges  were  assembled  with  four  active  arm  350  ohm 
gages  to  provide  maximum  transducer  output  and  optimum  calculation  of 
temperature,  shear,  and  centrifugal  effects.  The  flapwise  and  torsional 
gages  were  installed  at  the  quarter-chord  of  the  blade.  The  chordwise 
gages  were  installed  at  the  flapwise  neutral  axis  on  the  leading  and 
trailing  edge  of  the  D  spar.  All  bridge  wiring  was  routed  chordwise 
from  the  strain  gages  to  the  rear  of  the  D  spar  and  from  there  spanwise  to 
solder  tabs  at  the  root  end  of  the  blade.  Segmented  trailing  edge  pockets 
were  then  installed  on  the  spar  encasing  the  wires.  Thus,  the  instrumented 
blade  airfoil  contour  was  not  disrupted  except  where  the  strain  gage  itself 
was  installed. 

The  model  blade  strain  gage  bridges  were  calibrated  in  a  fixture  which  pro¬ 
vided  cantilever  restraint  at  the  blade  root  end.  Flapwise  and  chordwise 
forces  and  torsional  moments  are  applied  at  the  tip.  A  minimum  of  four 
discrete  loads  were  applied  in  both  the  positive  and  negative  directions. 
All  results  were  referenced  to  a  precision  series  resistance  calibration 
value.  The  fixture  shown  in  Figure  5  was  also  used  for  the  strain  gage 
calibration. 
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Blade  Hinge  Motions 


Instrumented  blade  flap  and  lag  motion  were  measured  by  using  film 
potentiometers.  A  pinion  and  sector  gear  mechanically  amplified  the  hinge 
motion,  providing  a  high  resolution  output  from  the  transducers .  The  flap 
and  lag  hinge  motions  were  calibrated  for  each  wind  tunnel  run  by  manually 
placing  the  blade  on  the  lead-lag  and  flapping  stops,  and  recording  the 
transducer  output  at  each  position. 

Model  Internal  Balance 


Six  components  of  rotor  force  and  moment  were  measured  on  a  Task  Corpora¬ 
tion  MK  XIX  strain  gage  balance.  The  balance  was  calibrated  by  applying 
known  pure  forces  and  moments  in  the  shaft  axis  system  at  the  rotor  head. 
A  minimum  of  five  discrete  loadings  were  applied  in  both  the  positive 
and  negative  directions  for  each  of  the  six  coordinates. 

Control  Positions 


Cyclic  and  collective  pitch  inputs  were  measured  by  potentiometers  which 
were  gear-driven  by  the  electromechanical  control  actuating  system. 

Shaft  Angle 

Rotor  shaft  angle  of  attack  was  measured  by  the  normal  equipment  provided 
for  that  purpose  in  the  wind  tunnel.  The  shaft  angle  is  controlled  by  the 
pitching  strut,  which  is  attached  to  the  aft  part  of  the  fuselage,  as  shown 
in  Figures  1  and  2.  The  vertical  position  of  this  strut  is  remotely 
controllable  from  the  tunnel  console,  causing  the  model  to  pitch  about 
trunnions  atop  the  two  main  struts  visible  in  Figures  1  and  2  .  The  shaft 
angle  readout  on  the  tunnel  console  was  calibrated  by  using  an  inclinometer. 
The  shaft  angle  variation  caused  by  the  internal  gimble  deflection  under 
the  applied  rotor  airloads  and  gravity  was  investigated  and  found  to  be 
negligible. 

Miscellaneous  Model  Instrumentation  and 
Monitoring  Equipment 

Timing  signals  for  defining  rotor  revolutions  and  azimuth  positions  were 
produced  by  a  photoelectric  device.  Once-per-revolution  and  72-per- 
revolution  pulses  were  recorded. 

Four  accelerometers  were  used  for  monitoring  internal  model  roll  motions, 
gearbox  vibration,  and  external  fuselage  yawing  motions. 

Motor  and  gearbox  temperatures  were  monitored.  Motor  electrical  power 
input,  water  and  oil  system  pressure,  and  gearbox  chip  detector  warnings 
were  available. 

Data  Acquisition  and  Monitoring 

The  principal  data  recording  device  used  for  the  test  program  was  a 
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fourteen  track  FM  tape  recorder.  The  system  used  narrow  band  FM  multiplex 
with  standard  subcarrier  oscillators  for  IRIG  channels  7  through  l6,  allow¬ 
ing  ten  channels  of  data  to  be  recorded  on  each  track.  Three  multiplexed 
direct  record  tracks  were  used,  allowing  all  dynamic  measurements  to  be 
recorded  simultaneously,  to  provide  proper  time  correlation  of  the  signals. 
Additionally,  one  track  was  used  for  audio  comments,  one  track  for  rotor 
azimuth  reference,  and  one  track  for  a  60-Hz  standard  time  reference. 

All  signal  conditioning  functions  were  accomplished  using  Sikorsky 
designed  electronic  modules. 

Rotor  forces  and  moments  measured  by  the  strain  gage  balance  were  displayed 
on-line,  using  a  Sikorsky-daveloped  Balance  Axis  Converter.  In  addition  to 
containing  the  signal  conditioning  module  for  each  of  the  bridges  in  the 
six  component  strain  gage  balance,  the  unit  contains  a  small  analog  computer 
which  performs  the  matrix  manipulation  to  convert  the  measured  strain  gage 
element  outputs  co  pure  forces  and  moments  in  the  shaft  axis  coordinate 
system  at  the  rotor  head.  The  unit  utilizes  third-order  Buttervorth  fil¬ 
ters  to  obtain  steady  readings  from  the  vibratory  signals .  Resultant 
forces  and  moments  in  engineering  units  (pounds  and  foot-pounds)  are  dis¬ 
played  on  six  digital  meters.  A  70  mm  camera  was  used  to  obtain  a  perma¬ 
nent  record  of  the  meter  readings  at  each  test  point. 

Critical  transducers  were  monitored  continuously  during  the  test  runs. 

The  instrumented  blade  flap  and  lag  motions  were  displayed  on  oscilloscopes. 
Any  FM  channel  could  be  displayed  on-line  for  study  of  amplitude,  wave 
form,  or  signal  quality.  A  four-trace  oscilloscope  synchronized  to  the 
instrumented  blade  1-per-revolution  signal  was  used  for  monitoring 
representative  blade  strain  gage  bridges. 

An  electronic  flapping  resolver  was  provided  in  the  control  console  for 
on-line  independent  display  of  longitudinal  and  lateral  tip  path  plane  tilt. 
The  unit  responded  to  first-harmonic  content  of  blade  flapping  motion  only. 
The  resolver  was  used  continuously  during  the  test  runs  in  conjunction  with 
the  cyclic  controls  for  trimming  the  rotor  tip  path  plane  either  perpen¬ 
dicular  to  the  rotor  shaft  (trimmed  flapping)  or  at  prescribed  Inclinations 
to  the  shaft. 


WIND  TUNNEL  TEST  PROCEDURE 


Aerodynamic  Tares 


A  series  of  wind  tunnel  runs  were  made  without  rotor  blades  to  measure 
aerodynamic  loadings  on  the  rotor  head  and  other  aerodynamic  ally  loaded 
components  of  the  metric  part  of  the  model.  These  loadings  were  deducted 
from  the  total  rotor  aerodynamic  loads  during  the  data  reduction  process. 
Thus  the  rotor  performance  results  reported  herein  are  due  to  the  blades 
only.  This  procedure  is  followed  because  no  special  attempt  has  been  made 
within  the  scope  of  this  or  previous  programs  to  scale  the  hub  configura¬ 
tion. 


The  aerodynamic  tares  were  determined  for  the  entire  range  of  shaft  angles 
and  forward  speeds  planned  for  the  test. 


Typical  Wind  Tunnel  Run  Procedure 


The  typical  wind  tunnel  run  began  with  a  series  of  calibration  records. 
These  records  were  taken  with  the  rotor  and  tunnel  turned  off  and  are  as 
follows : 


1.  A  "ZCAL"  or  electrically  zero  output  record  was  recorded  for 
each  transducer. 

2.  An  "RCAL"  or  known  electrical  change  from  the  electrical  zero 
condition  was  recorded  for  each  condition.  The  difference  between 
the  "RCAL”  and  "ZCAL”  records  was  electrically  equivalent  to  the 
transducer  response  to  a  certain  amount  of  physical  moment  or 
hinge  motion. 

3.  At  least  two  "XCAL"  records  were  taken.  A  low-lead  record 
was  taken  with  the  blade  manually  placed  on  the  lower  flapping 
and  lead  stops.  Then  a  high-lag  record  was  taken  with  the  blade 
held  on  the  upper  flapping  and  lag  stops.  These  records  were  used 
to  provide  physical  calibrations  for  the  hinge  motion  transducers 
by  using  the  known  blade  hinge  positions  on  the  physical  motion 
stops . 

4.  Three  or  more  "static  zero"  records  were  taken  to  obtain 
gravity  tares  for  the  rotor  balance  at  various  shaft  angles. 

5.  The  rotor  was  then  started  and  brought  up  to  an  operating  speed 
to  be  used  during  the  run.  Then  a  "dynamic  zero”  record  was  taken 
at  zero  lift  in  hover  to  provide  a  zero  reference  for  the  blade 
strain  gages.  It  was  more  useful  to  consider  this  condition  rather 
than  a  static  condition  as  a  moment  reference,  because  of  the 
fairly  large  steady  drift  in  blade  bending  bridge  output  which 
occurs  as  the  rotor  is  brought  up  to  speed. 

6.  Dynamic  zeros  were  repeated  if  more  than  one  rotor  speed  was 
to  be  used  during  the  run. 

7.  The  tunnel  was  brought  up  to  operating  speed  and  the  rotor 
controls  were  set  to  provide  the  desired  operating  condition. 

The  balance  axis  converter  was  used  to  set  the  rotor  at  a 
preselected  lift  for  some  of  the  conditions.  In  others,  the  col¬ 
lective  pitch  was  set  to  a  desired  value.  In  either  case  the 
cyclic  pitch  was  varied  to  provide  zero  or  some  prescribed  amount 
of  first  harmonic  flapping  at  the  hinge.  The  shaft  angle  was 
also  set  at  a  prescribed  value.  When  the  desired  operating 
condition  was  reached,  data  was  recorded. 

8.  When  the  data  points  in  the  run  were  complete,  the  dynamic, 
static,  ZCAL,  RCAL  and  XCAL  records  were  repeated.  These  final 
zero  and  calibration  records  were  used  to  check  for  transducer 
signal  drifting,  which  was  usually  no  more  than  a  few  percent 
of  the  full-scale  reading. 


9.  After  the  data  acquisition  system  vas  connected  to  each 
blade  configuration,  records  were  taken  on  tape  as  the  blade 
vas  manually  bent  and  twisted  in  the  positive  flapvise,  chordvise 
and  torsion  directions.  This  provided  an  independent  check  on 
data  signal  polarity. 

HELICOPTER  MODEL  TEST  DATA  REDUCTION  PROCESS 
Rotor  Performance  Parameters 

The  data  reduction  process  for  rotor  performance  consisted  mainly  of 
applying  certain  corrections  to  data  recorded  on  line,  performing  some 
relatively  simple  conversions  to  coefficient  form,  and  calculating  some 
of  the  more  common  nondimens ioned  rotor  parameters. 

The  quantities  recorded  on-line  for  each  operating  condition  included 
tunnel  velocity  in  knots (V),  rotor  revolutions  per  minute  (N),  shaft 
(or  fuselage)  angle  of  attack  in  degrees  (afr),  tunnel  stagnation 
temperature  in  degrees  farenheit  (Tgg),  and  the  six  balance  axis  converter 
outputs,  Xjjac»  yBAC»  zBAC»  lBAC»  MBAC*  NBAC* 

The  rotor  shaft  load  components  due  to  blade  aerodynamic  forces  were 
obtained  from  the  balance  axis  converter  readings  by  using  the  following 
equations . 

Xs  -  (Xbac  -  HX  -  X0)/NSFX 
Ys  -  (YBAC  -  Hy  -  Y0)/NSfy 
z8  »  (zbac  -  Hz  -  zo)/Nsfz 
Ls  *  (L^q  -  Hl  -  Lq)/SsFL 
Ms  ■  (MflAC  -  hM  -  Mo)/NgFM 
Ns  *  (NfiAC  “  hN  -  No)/NSFn 

The  quantities  Hx  through  Hr  in  the  above  equations  represent  aerodynamic 
tares,  and  gravity  tares  due  to  changes  in  shaft  angle  of  attack.  The 
quantities  Xq»  Yq»  Zo  and  Lq,  Mq,  Nq  are  the  averages  of  balance  axis 
converter  readings  taken  at  zero  angle  of  attack  static  zero  conditions 
before  and  after  each  wind  tunnel  run. 

The  control  angles  were  corrected  for  the  svashplate  phase  angle  ^SP» 
which  vas  16.6  degrees. 

Als  "  cos  (*SP)  +  BlsR  «in  (*sp) 

*Lb  "  BisR  cos  (*sp)  -  AlsR  sin  (♦sp) 


The  advance  ratio  is  calculated  from 

y  *  V  x  1.688 
OR 

The  density  is  calculated  by  using  the  following  equation  for  compressible, 
adiabatic  flow: 


The  rotor  thrust  coefficient  is  calculated  from: 


The  angle- of- attack  tunnel  vail  correction  is  calculated  by  using  an 
expression  from  Reference  10: 

do  ■  -.55  Mi  /- Ct  \ 

*  w 

a8  “  Qfr  ♦  La 

The  Zs  and  Xs  forces  are  resolved  into  the  wind  axis  system: 

Lr  *  -Zg  cos  a8  ♦  Xg  sin  as 
Dr  *  -Zq  sin  a8  -  Xg  cos  a8 

The  rotor  force  and  moment  coefficients  are  calculated  from 
CL  _  Lr 

CD  dr 

-*fT 

Cy  ys 

o  RFn 
Crm  Ds 
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giWott 


where 


CpM  *  Ms 
—  RFq 


fq  *  pnn2R4o 

The  rotor  control  angle  of  attach  is 
ac  *  “s  ”  Bls 

The  values  of  first-harmonic  flapping  components  at  the  hinge  were 
determined  with  the  digital  uata  processing  program  for  operating  condi¬ 
tions  in  which  these  quantities  were  not  set  equal  to  zero.  This  program 
will  be  discussed  in  a  subsequent  section  of  this  report. 

The  blade  tip  Mach  number  at  the  90°  azimuth  is  determined  from: 

OR  ♦  1.688V  , 

m1,90  * 


The  values  of  Vs  and  fisR  are  the  simulated  full-scale  equivalent  speeds. 
The  conversion  factor  2.0  between  the  actual  and  simulated  speeds  is 
appropriate  to  the  1/U  stiffness  models. 

V8  «  2.0  x  V 

flsR  «  2.0  x  OR 


Reduction  of  Blade  Response  Data 


The  blade  response  data  on  the  FM  multiplex  analogue  tapes  which  were 
recorded  during  the  wind  tunnel  runs  were  converted  to  Sikorsky  format 
digital  tapes.  The  digital  data  were  then  placed  temporarily  on  disc  stor¬ 
age  for  computer-aided  editing  and  quality  checking.  Upon  completion  of 
this  task,  the  samples  of  blade  response  data  were  converted  to 
engineering  units  and  placed  on  standard  FORTRAN  readable  tapes.  One  tape 
was  written  for  each  of  the  fou.  blade  configurations.  The  blade  response 
data  consist  of  response  parameter  values  taken  at  5°  azimuth  Intervals 
over  a  forty-revolution  record  length.  Rotor  performance  parameters  for 
each  of  the  operating  conditions  were  also  placed  on  the  tape  for  each 
configuration. 


A  computer  program  was  written  to  extract  the  blade  response  data  from  the 
final  tapes  and  perform  the  most  usual  types  of  processing.  The  retrieval 
computer  program  and  copies  of  the  final  dnta  tapes  are  available  from 
the  Eustls  Directorate  of  the  U.  S.  Army  Air  Mobility  Research  and 
Development  Laboratories. 
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Standard  Sikorsky  procedures  and  equipment  were  used  to  vrite  digital 
tape  records  from  the  FM  multiplex  analogue  tapes  taken  during  the  wind 
tunnel  runs . 


The  digital  records  for  each  rotor  operating  condition  and  transducer 
consisted  of  values  observed  at  azimuth  intervals  of  5  degrees,  starting 
at  zero  azimuth  as  the  instrumented  blade  passed  over  the  tail  of  the 
model,  and  continuing  for  more  than  forty  revolutions.  The  data  were 
expressed  in  noimalized  computer  units.  All  calibration  and  reference 
zero  records  were  converted  in  the  same  manner,  except  that  external  timing 
pulses  vere  used  for  calibration  records  during  which  the  rotor  vas  not 
turning.  For  records  where  the  rotor  was  turning,  the  digital  conversion 
intervals  were  governed  by  the  once-per-revolution  and  the  72-per-revolu  - 
tion  (every  5  degrees)  pulses  provided  from  the  model.  The  normalized  com¬ 
puter  unit  values  on  the  Sikorsky  digital  tape  were  the  equivalent  of 
decimal  integer  values  with  a  range  from  0  to  512,  because  of  the  particu¬ 
lar  word  length  utilized. 

Quality  Checking  and  Editing  of  Blade  Response  Data 

The  blade  response  data  were  placed  on  temporary  disc  storage  for  editing 
and  quality  checking. 

The  basic  quality  check  option  calculated  and  printed  out  the  maximum 
data  value,  minimum  data  value,  average  data  value,  and  root-mean-square 
deviation  from  average  of  the  data  in  a  specified  time  history  record 
from  a  certain  transducer.  These  quantities  were  output  in  normalized 
computer  units  and  in  engineering  units.  Data  discontinuties  were  found 
and  counted  by  comparing  actual  and  interpolated  values,  as  discussed  later. 
The  number  of  data  values  in  each  record  which  reached  the  maximum  and 
minimum  value  vere  also  counted  to  give  an  Indication  as  to  how  much  data 
in  a  questionable  record  might  be  off  scale.  The  quality  check  vas  carried 
out  for  specified  groups  of  records  and  channels.  The  quality  checking 
program  also  contained  options  for  printout  of  sample  time  histories  and 
frequency  analyses,  to  investigate  and  diagnose  any  questionable  data. 

The  computerized  quality  checking  vas  supplemented  with  inspection  of  the 
oscillograph  records  made  when  digital  conversion  took  place. 

The  mathematical  operations  for  the  quality  check  options  are  discussed 
in  the  following  paragraphs. 

Conversion  From  Computer  Units  to  Engineering  Units 


An  average  or  individual  data  value  in  normalized  computer  units  Wp  is 
converted  to  engineering  units  by  using  the  following  transformation: 
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where 


W£  is  the  desired  value  in  engineering  units. 

Be  Is  the  physical  equivalent  of  the  electrical  difference  between 
the  "ZCAL"  and  "RCAL"  records  on  the  tape. 

Wr  is  the  average  computer  unit  value  in  the  "RCA L"  records  on  the 
digital  tape  for  the  vind  tunnel  run. 

Wz  is  the  average  computer  unit  value  in  the  "ZCAL”  records  on 
the  digital  tape  for  the  vind  tunnel  run.  Wo  is  the  average 
computer  unit  value  in  the  zero  reference  records  on  the  digital 
tape  for  the  vind  tunnel  run.  Fc  is  an  arbitrary  offset  constant 
which  is  used  when  the  physical  reference  condition  is  not  at 
a  zero  data  value,  or  when  there  is  a  need  for  a  zero  drift 
correction  or  the  like. 

The  values  Re,  Wr,  Wz,  Wq  and  Fc  remain  constant  for  a  data  time  history 
segment.  Different  values  are  obtained  for  each  transducer,  of  course. 

If  post  run  calibrations  and  zeros  are  invalidated  by  transducer  failure, 
the  prerun  calibrations  and  zeros  only  are  used. 

Detection  of  Discontinuities 


In  order  to  check  data  for  spiking  caused  by  incipient  instrumentation 
failure,  data  values  Xj  in  the  time  history  records  are  compared  with 
values  Xj  Interpolated  from  neighboring  values: 

Xj  -  2  (Xj+1  «■  Xj-i)-l  (Xj+2  ♦  Xj_2) 

3  Z 

A  discontinuity  is  counted  whenever 


XJ  -XJ 


NpD 


The  value  Nrd  was  set  by  trial  and  error  by  raising  the  value  of  Njd  when 
discontinuities  were  counted  in  an  acceptable  record.  A  value  of  8  was 
found  suitable  for  the  blade  bending  and  torsion  data.  A  value  of  16 
was  suitable  for  the  hinge  motions. 


Frequency  Analysis 

The  frequencies  at  which  amplitude  components  are  to  be  found  are  defined 
by  a  starting  frequency,  (Fg),  an  increment  UFrq),  and  the  number  of 
increments,  (Npj  ♦  l) *  Note  that  frequencies  are  expressed  in  cycles 
per  revolution. 

(Frq)j  *  Fs  ♦  (AFrq)(J-1)  i  -  1,2,3  ...(NfI  ♦  1) 
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The  time  history  segment  for  frequency  analysis  begins  with  the  start  of 
the  rotor  revolution  denoted  by  Npg»  counting  from  the  beginning  of  the 
recorded  sample. 

The  time  history  segment  for  frequency  analysis  ends  with  the  end  of  the 
rotor  revolution  denoted  by  Nfe»  counting  from  the  beginning  of  the 
recorded  sample. 

The  index  k  takes  the  values  1,  2,  3  . ..  K»  and  counts  digital  data  values 
in  the  time  history  segment  for  frequency  analysis,  starting  at  azimuth 
zero  of  the  Kfs  revolution. 


The  integer  K  is  the  maximum  value  of  k  and  is  defined  by 
K  *  (l  +  nfe-nfs)nprv  +  1 

where  Nppy  is  the  number  of  samples  per  rotor  revolution  (72).  K  refers 
to  the  first  data  value  in  the  (NpE+i)  th  revolution. 

(WT)kis  the  kth  digital  data  value  in  the  segment  starting  with  the 
Npgth  revolution  of  the  sample. 


The  amplitude  components  A j ,  Bj  and  Rj  desired  are  then  given  by  the 
following  expressions: 

f  r 

K-l 


J  NPrv(1  ♦  Nra  - 


V. 


1  <"r>X  “•  ^ 

k-2,3,U 


<“T>1  *  ~  [('ta)j  1J. 


r 


B 


J  NpRy(l  ♦  Npg  -  Nps) 


(  “z  ) 


K-l 


I  (WT)k  .tn[(r.q)J  g»Qt-l)| 

k-2,3,1*  L  "pRV  J 

(Wj )lt  *in  ["(Frq)j  2ir(  ic  -  l)1 
L  nprv  J 


"j  -  <v  ♦ 


The  above  equations  are  numerically  evaluated  equivalents  of 


a(u1) 


W£(t)cos  (wj*)d<l 
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b(u><)  *  1 
J  T 

0»</> 

r(uj)  «£a2(wj)  +  b2(wj)] 

Where  u«  is  a  frequency  ratio  with  respect  to  rotational  frequency  Q 
Wg  is  tine  history  data  in  engineering  units,  and  is  the  azimuth. 

The  accuracy  of  the  numerical  frequency  analysis  depends  on  the  number  of 
samples  per  revolution  and  on  the  length  of  the  time  history  record  if 
the  frequencies  to  be  analyzed  are  nonharmonic.  For  the  72-per-revolution 
sampling  rate,  nonharmonic  frequencies  beyond  15  cycles  per  revolution 
are  not  recommended.  For  the  UO-revolution  time  history  record  length, 
frequencies  belov  .15  per  revolution  are  not  recommended. 

Generation  of  Final  Data  Tapes 

When  quality  checking  of  the  response  data  was  complete,  data  samples 
for  each  condition  number  vere  specified  and  converted  to  engineering  units 
as  indicated  on  page  26 .  Forty  revolutions  of  time  history  data  for  each 
operating  condition  v<*re  placed  on  the  final  tape. 

The  prerun  and  post  run  calibration  data  vere  normally  averaged  in  the 
conversion  to  engineering  units.  If  a  strain  gage  transducer  failed 
during  a  wind  tunnel  run,  only  the  prerun  calibration  data  vere  used. 

The  electrically  generated  "ZCAL"  and  "RCAL"  records  vere  used  to  cali¬ 
brate  the  strain  gage  transducer  rate  of  change  of  engineering  units 
vith  respect  to  computer  units,  as  it  appears  in  the  equation  on  page  26, 
The  physical  zero  reference  for  the  strain  gage  bridge  measurements  of 
blade  bending  and  tvlsting  vat  the  average  transducer  output  during  the 
zero  lift  hover  condition.  This  physical  zero  condition  avoided  the 
substantial  drift  in  bending  and  torsion  bridge  output  vhich  occurs  as 
the  rotor  is  brought  up  to  operating  rotational  speed.  In  some  cases 
the  chordvlse  bending  bridges  displayed  a  moderate  amount  of  zero  drift 
during  a  vind  tunnel  run.  This  vaa  compensated  by  adding  appropriate 
Fc  values  for  the  operating  conditions  in  the  run. 

Physically  generated  transducer  signals  vere  provided  for  the  calibration 
of  blade  hinge  motions,  by  manually  placing  the  blade  against  the  motion 
stops  at  the  beginning  and  end  of  each  run. 

Data  Retrieval  and  Display  Program 

The  data  retrieval  and  display  program  retrieves  data  from  the  final  tapes 
according  to  operating  condition,  physical  measurement  channel,  end 
rotor  revolution. 


/ 


We(\I>)  sin  )  d<J/ 


-? 


* 
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Present  processing  options  are  direct  time  history  sampling)  frequency 
analy  *s,  maximum  and  minimum  values  summary)  average  revolution,  and,, 
harmonic  - . 

The  data  retrieved,  and  display  program  documentation  and  copies  of  the 
fined,  data  tapes  are  available  from  the  Eustis  Directorate  of  the  U.  S. 

Army  Air  Mobility  Research  and  Development  Laboratory. 

The  final  data  tapes  are  FORTRAN-readable  to  facilitate  the  provision  of 
additional  date  processing  modes  that  may  be  required  in  the  future. 

The  users'  instructions  for  the  data  retrieval  program  are  presented  under 
Appendix  II.  It  is  expected  that  these  instructions  will  be  used  vlth 
the  descriptions  appearing  in  the  next  report  section,  "TEST  DATA  OBTAINED1.' 

The  data  processing  options  are  straightforward  and  require  little  addition¬ 
al  explanation,  since  the  data  on  tape  is  already  in  engineering  units. 

The  frequency  analysis  option  is  identical  to  the  similar  quality  check 
calculation  described  earlier  on  page  28,  except  that  the  terms  Rg/tWR- 
Wz)  are  not  required  to  convert  to  engineering  units. 

TEST  DATA  OBTAINED 


The  reduced  data  tape  for  each  blade  configuration  contains  19  rotor  per¬ 
formance  parameters  and  1*0  revolutions  of  digital  time  history  data  for 
lb  blade  response  transducers  for  each  rotor  operating  condition.  Data 
is  expressed  in  coefficient  form  or  in  engineering  units. 

The  reduced  data  tapes  for  the  dynamically  scaled  fiberglass  blades  contain 
85  conditions  each.  The  tape  for  the  three-times-scale  stiffness  aluminum 
blade  contains  93  conditions.  The  condition  numbers  refer  to  similar 
conditions  on  all  four  H-3b  rotors,  to  facilitate  comparative  studies. 

The  various  test  conditions  for  each  rotor  configuration  are  listed  in 
Table  XVI  with  descriptive  terms  that  indicate  the  purpose  of  each 
condition  or  series  of  conditions.  Note  that  a  particular  phenomenon  being 
explored,  such  as  blade  stall,  does  not  necessarily  occur  at  each 
condition  of  a  series  of  conditions  so  designated.  The  conditions  are 
generally  numbered  in  the  order  in  which  they  were  taken  during  the  wind 
tunnel  runs. 

The  digital  time  history  samples  are  provided  in  conventional  computer 
language  as  discrete  values  which  begin  at  the  zero  azimuth  (blade  over 
the  tail)  and  proceed  at  5-degree  intervals  for  bo  rotor  revolutions. 

The  number  of  unavoidable  strain  gage  transducer  failures  was  minimized 
by  running  the  more  lightly  loaded  conditions  first.  The  conditions  and 
transducer  channels  for  which  the  strain  gages  were  operative  are  given 
in  Tables  XVII  through  XX. 
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The  time  histories  of  blade  response  on  the  reduced  data  tapes  are  1+0 
revolution  samples  of  the  actual  tine  histories.  Noticeable  amounts  of 
nonharmonic  notion  are  present  in  sane  of  the  test  data.  Flapvise  and 
torsion  response  exhibited  an  apparently  randan  nonharmonic  response 
vhen  rotor  stall  conditions  vere  approached  and  exceeded.  For  example,  the 
.65R  flapvise  bending  and  .20R  torsion  responses  of  the  zero-degree  tvist 
dynamically  scaled  fiberglass  spar  blade  contained  variations  of  this  type 
vhich  vere  ±10  percent  of  the  average  half  peak-to-peak  amplitude, 
vhen  operating  at  condition  7.  In  addition  to  the  apparently  random 
response,  the  cbordvise  strain  gages  and  lead-lag  motions  contain  excita¬ 
tions  from  the  rotor  drive  system.  Frequency  analysis  results  for  the  .20R 
chordvise  bending  moment  response  in  zero  lift  hover  are  given  in  Table 
XXI.  These  nonharmonic  components  (except  for  the  small  7-per-revolution 
component)  are  effectively  removed  by  the  revolution-averaging  procedure. 

If  required,  the  unaveraged  time  history  can  be  corrected  for  the  drive 
system  excitation  by  analyzing  'or  the  particular  frequency  components 
appearing  in  Table  XXI,  and  then  subtracting  those  results  from  the  time 
history  as  recorded. 

The  reduced  data  tapes  and  a  computer  program  to  retrieve  the  data  are 
available  from  the  Eustis  Directorate  of  the  U.  S.  Army  Air  Mobility 
Research  and  Development  Laboratory. 

Rotor  Performance  Data 


A  complete  listing  of  all  performance  data  is  provided  for  each  of  the 
four  configurations  in  Tables  XXII  through  XXV.  Sample  performance  data 
are  plotted  in  Figures  21  through  32.  Comparisons  are  made  betveen  the 
performance  of  the  four  H-31*  model  rotors.  Comparisons  are  also  made 
betveen  the  performance  of  the  H-3U  model  rotors  and  the  corresponding 
full-scale  results  reported  in  References  2  and  3.  Correlation  with 
calculated  results  shown  in  these  plots  will  be  discussed  in  a  subsequent 
part  of  this  report. 

Torque  versus  drsg  crossplots  provided  in  Figure  21  show  the  comparative 
performance  levels  of  the  four  model  rotor  configurations.  The  variations 
of  torque  versus  drag  for  the  various  configurations  fulfill  qualitative 
expectations.  The  flapped  airfoil  blade  is  substantially  better  in  its 
performance  characteristics  than  the  other  tvo  dynamically  scaled  blades. 
The  three-times-scale  stiffness  aluminum  blade  performance  becomes 
substantially  worse  than  that  of  the  dynamically  scaled  blades  as  the 
rotor  reaches  a  stalled  condition. 

Sample  comparisons  of  drag  versus  lift  and  torque  versus  lift  for  one 
shaft  angle  and  advance  ratio  are  provided  in  Figures  22  and  23.  The 
model  rotor  data  variations  are  qualitatively  similar  to  the  full-scale 
data,  vhich  vas  obtained  from  Reference  3  •  The  quantitative 
variations  betveen  the  model  and  the  full-scale  performance  data  can  be 
attributed  to  Reynolds  number  effects. 
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Figures  24  through  27  display  model  rotor  performance  at  various  advance 
ratios  and  shaft  angles  included  in  the  test  program.  The  variations 
shown  in  these  plots  are,  in  general,  qualitatively  similar  to  the  data 
from  Reference  3.  The  intersection  or  close  approach  of  the  curves  in 
Figures  24(a)  for  ag  »  2°  and  os  *  0°  at  low  lift  was  typical  for  the 
fiberglass  spar  dynamically  scaled  models  with  this  small  increment  in 
shaft  angle.  The  three- times- scale  stiffness  aluminum  spar  model  did  not 
exhibit  this  tendency.  Figures  24  through  27  can  be  used  to  provide  initial 
trim  values  of  shaft  angle  and  control  position  for  blade  response 
calculation  correlation  studies,  when  a  desired  rotor  lift  and  drag  is 
given.  (Note  that  lift  and  cyclic  pitch  versus  collective  pitch  for  an 
advance  ratio  ofy*.5  are  presented  in  Figures  2d  and  29  for  comparison 
with  full-scale  data). 

Figures  28,  29  and  30  compare  model  rotor  performance  parameter  variation 
versus  collective  pitch  with  corresponding  full-scale  data  from  Reference 
3  ,  and  with  calculated  data.  Plots  of  this  type  are  of  special  interest 
for  the  comparison  of  model  data  with  full-scale  data  and  with  calculated 
data. 

The  comparison  between  model  and  full-scale  lift  versus  collective  data 
shows  that  the  model  develops  more  lift  beyond  stall  than  would  be 
expected  from  a  comparison  of  full-scale  and  model  airfoil  two-dimensional 
steady-flow  airfoil  data.  The  model  airfoil  two-dimensional  steady-flow 
maximum  coefficient  of  lift  Bhown  in  Figure  Ik  is  about  60  percent  of 
full-scale  values.  The  model  rotor  lift  shown  for  a  collective  pitch  of 
10  degrees  in  Figure  28  is,  however,  about  85  percent  of  the  full-scale 
value.  Plots  of  the  drag,  torque  and  cyclic  pitch  data  against  collective 
pitch  display  good  qualitative  and  reasonable  quantitative  agreement. 

The  comparisons  between  model  and  full-scale  results  in  Figures  22,  23, 
and  28  through  30  show  that  model  and  full-scale  results  are  very  similar 
qualitatively,  and  that  the  quantitative  effects  of  Reynolds  number  and 
Mach  numbers  are  milder  than  would  be  expected  from  a  comparison  of  full- 
scale  and  model  steady-flow  two-dimensional  airfoil  data. 

Figures  31  and  32  are  included  here  to  provide  comparative  test  and 
calculated  data  for  the  remaining  two  model  configurations.  The  comparat¬ 
ively  small  value  of  Bls  required  for  the  flapped  blade  is  attributed  to 
aeroelastic  twisting  on  the  advancing  part  of  the  azimuth. 

The  significance  of  the  calculated  data  shown  on  Figures  21  and  28  through 
32  will  be  discussed  under  a  subsequent  part  of  this  report. 


Blade  Response  Sample  Data 


Figure  33  presents  an  example  of  the  variation  of  the  algebraic  maximum 
and  minimum  bending  and  twisting  response,  as  a  function  of  rotor  lift 
for  the  four  blade  configurations. 

Bending  response  of  the  fiberglass  spar  dynamically  scaled  rotors  is  re¬ 
latively  similar.  The  untwisted  aluminum  spar  three- times- scale  stiffness 
blade  has  a  much  higher  bending  moment  response.  This  is  an  indication 
that  the  bending  displacement  of  this  blade  is  much  more  than  one-third 
as  large  as  the  dynamically  scaled  blade  deflections. 

Even  at  the  moderate  scale  forward  speed  of  124  knots,  the  flapped  blade 
has  a  much  larger  torsional  moment  than  the  other  blades,  except  for 
stalled  conditions.  Thus  the  performance  benefits  of  this  blade  would 
have  to  be  considered  as  a  trade-off  against  continuous  high  pushrod 
loadings . 

Figure  34  presents  sample  time  histories  reproduced  from  Reference  1, 
along  with  model  data  from  the  present  investigation,  which  was  converted 
to  equivalent  full-scale  stress  by  using  the  factors  given  on  page  7 
of  this  report. 

The  agreement  between  model  and  the  full-scale  bending  response  is  quite 
good,  considering  the  difference  between  operating  conditions.  The  most 
notable  difference  in  bending  response  is  the  high  peaks  in  flapwise 
bending  at  .80R. 


The  model  torsional  response  at  the  lower  lift  shown  in  Figure  34(c)  has 
some  similarities  to  the  full-scale  data.  A  stall-induced  oscillation 
occurs  on  the  model  retreating  blade  at  the  higher  lift  condition  shown. 

The  retreating  blade  oscillation  does  not  appear  on  the  full-scale  blade, 
presumably  because  of  less  stalling. 

Sample  average  revolution  time  history  data  are  presented  in  Figures  35 
through  46.  Blade  bending  and  torsion  response  are  presented  for  conditions 
a  advance  ratios  of  w*.35»  .40,  and  .50.  Other  blade  response  samples 
are  presented  and  compared  with  calculated  results  under  the  report  section 
dealing  with  the  normal  modes  aeroelastic  analysis.  The  conditions  chosen 
for  these  samples  include  some  of  those  for  which  the  blades  were  expected 
to  be  close  to  or  intersecting  the  rotor  wake. 


Flapwise  bending  moments  response  is  presented  in  Figures  35  through  38. 

The  flapwise  bending  response  for  the  -8°  twist  blade  is  not  consistently 
higher  than  the  response  for  the  dynamically  scaled  untwisted  blade  because 
of  the  moderately  high  lift  coefficient.  This  trend  is  noticeable  in  the 
full-scale  results  presented  in  Figures  67  through  69  of  Reference  2.  The 
flapwise  bending  results  for  the  untwisted  5- degree  trailing -edge  flap 
blade  are  very  similar  to  those  for  the  -8°  twist  blade,  possibly 
because  of  a  steady  aeroelastic  twist.  The  three-time^-scale  stiffness 
aluminum  spar  blade  flapwise  moments  are  much  higher  than  those  on  the 
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dynamically  scaled  fiberglass  spar  blades,  but  with  a  similar  vave  form. 

Chordwise  bending  moment  response  is  sampled  in  Figures  39  through  42. 
Chordwise  bending  moment  comparisons  show  the  reduced  response  expected 
for  the  >8°  twist  blade,  as  compared  to  that  for  the  0°  twist  blade.  The 
trends  are  again  similar  to  the  corresponding  full-scale  results  presented 
in  Figures  70  through  72  of  Reference  2.  The  chordwise  bending  results 
for  the  trailing-edge  flap  blade  are  similar  to  the  results  for  the  -8° 
twist  blade.  The  three-times -scale  stiffness  aluminum  spar  blade  chord- 
vise  moments  are  much  higher  than  the  dynamically  scaled  fiberglass  spar 
blade  moments,  and  with  a  much  different  wave  form. 

Torsional  moment  response  is  sampled  in  Figures  43  through  46.  The  tor¬ 
sional  response  sanple  for  the  untwisted  blade  at  y  ■  .35  in  Figures  43(a) 
shows  an  aggravated  response  in  the  azimuth  region  around  i>  *  340°,  which 
may  be  due  to  the  wake  vortex  passage  frequency,  which  is  approximately 
equal  to  the  torsional  natural  frequency  at  this  advance  ratio.  The  tor¬ 
sional  response  in  the  azimuth  region  around  M  340°  is  much  smaller  at 
advance  ratios  of  U  ■  .40  and  y  »  .30,  as  shown  in  Figures  43(b)  and  8l(a) 
respectively.  At  these  advance  ratios,  the  vortex  passage  frequency  is 
respectively  lower  and  higher  than  the  blade  torsional  frequency 

The  torsional  response  of  the  three-times-scale  stiffness  aluminum  spar 
blade  shows  a  higher  frequency  which  is  approximately  equal  to  the  tor¬ 
sional  natural  frequency,  and  which  is  characteristic  of  a  stall  flutter 
response.  A  lower  frequency  is  also  present,  and  this  may  be  due  to  wake 
vortex  forcing  or  bending-torsion  coupling. 
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NORMAL  MODE  AERO  ELASTIC  ANALYSIS  CALCULATIONS 


GENERAL  DESCRIPTION  OF  NORMAL  MODES  AEROELASTIC  ANALYSIS 


The  normal  mode  aeroelastic  analysis  utilized  in  the  work  described  herein 
is  based  on  a  step-by-step  timevise  integration  of  the  elastic  rotor  blade 
equations  of  motion,  which  are  documented  in  Reference  11.  These  equations 
of  motion  are  expressed  in  terms  of  uncoupled  flapwise  bending,  chordwise 
bending,  and  torsion  natural  vibration  modes.  The  modes  are  calculated 
according  to  the  methods  presented  in  References  12  and  13.  Articulated 
blade  rigid  body  flap  and  lag  hinge  motions  are  also  included  as  separate 
degrees  of  freedom. 

When  a  steady-state  rotor  operating  condition  is  being  analyzed,  the  inte¬ 
gration  of  the  equations  of  nation  proceeds  in  small  but  finite  timevise 
steps  from  arbitrary  starting  values.  After  a  number  of  rotor  revolutions, 
the  calculated  blade  motions  typically  become  cyclic  within  a  desired 
tolerance.  At  this  point,  time  variations  of  blade  bending  and  twisting 
loads  or  stresses  are  calculated  from  the  blade  response  shapes.  Rotor 
performance  parameters  are  calculated  by  integrating  blade  local  aero¬ 
dynamic  loadings  over  the  rotor  disc,  in  a  manner  similar  to  that  used  in 
Reference  lU. 

The  rotor  blade  is  represented  analytically  by  up  to  15  spanvise  segments. 
The  mass  ar  elastic  properties  of  each  segment  are  defined  by  the  local 
properties  of  the  actual  rotor  blade. 

Aerodynamic  loading  calculations  for  the  blades  consider  the  important 
effects  of  unsteady  flow  in  the  manner  reported  in  Reference  15.  The  local 
aerodynamic  lift  and  pitching  moment  coefficients  at  each  segment  and  time 
step  are  determined  by  a  table  look-up  procedure,  with  local  angle  of 
attack,  nondimens ional  time  rate  of  change  of  angle  of  attack,  and  non- 
dimensional  pitch  acceleration  as  parameters.  The  table  of  unsteady  aero¬ 
dynamic  coefficients  was  derived  from  oscillating  airfoil  test  data  pre¬ 
sented  in  References  15  and  16.  A  table  of  static  (steady  flow)  airfoil 
data  is  automatically  used  for  angle  of  attack  and  speed  regimes  for  which 
unsteady  aerodynamic  data  are  unavailable. 

The  calculated  blade  response  cases  presented  herein  utilized  steady  flow 
data  for  local  blade  drag  coefficients. 

As  cm  option,  inflow  velocities  which  vary  over  the  rotor  disc  may  be 
utilized  in  the  normal  modes  analysis.  The  variable  velocity  input  re¬ 
places  a  basic  uniform  inflov  velocity  option  which  utilizes  the  classical 
momentum  theory.  A  method  for  calculating  inflov  velocity  due  to  a  vortex 
wake  vith  a  prescribed  helical  geometry  is  presented  in  Reference  17.  The 
method  used  herein  is  similar,  but  computational  efficiency  is  improved  by 
representing  shed  wake  effects  with  modifications  in  the  shapes  of  local 
lift  curve  shapes.  A  computer  program  to  calculate  inflov  velocity  distri¬ 
bution  in  this  fashion  is  being  provided  to  the  Government  under  the  pro¬ 
visions  of  Contract  NAS  2-6463.  The  work  of  this  contract  is  being  per- 
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formed  by  Sikorsky  Aircraft  and  the  United  Aircraft  Research  Laboratories 
for  the  National  Aeronautics  &  Space  Administration  Ames  Research  Center. 

The  normal  modes  analysis  has  been  extended  to  consider  the  motions  of  a 
rotor  system  with  up  to  seven  blades ,  attached  to  a  rigid  fuselage  with  six 
degrees  of  freedom.  This  modification  is  described  ir  References  13,  18, 
and  19,  and  was  provided  to  the  Government  under  the  provisions  of  Contract 
DAAJ02-71-C-002i* . 

COMPARISONS  OF  NORMAL  MODES  AERO  ELASTIC  ANALYSIS  CALCULATIONS  WITH  H-3U 
MODEL  TEST  RESULTS 


Application  of  Analysis 

For  the  purpose  of  comparison  between  model  test  and  normal  modes  aero- 
elastic  analysis  calculations,  the  following  parameters  wrre  considered  as 
independent  variables  defining  an  analytical  operating  condition: 

1.  Forward  velocity  (V  or  Va) 

2.  Rotational  speed  (QR  or  ftgR) 

3.  Shaft  Angle  (a8) 

1*.  Lift  coefficient  -  solidity  ratio  (CL/o) 

5.  First -harmonic  flapping  components  at  hinge  (als  and  bla) 

At  each  analytical  condition,  rotor  control  inputs  (A^g,  Bis  and  0C)  or 
aerodynamic  scaling  factors  were  manipulated  to  produce  a  reasonable  match 
between  the  analytical  values  of  the  above  parameters  and  chose  existing 
at  some  test  condition. 

Choice  of  the  above  set  of  parameters  caused  the  test  and  analytical  blade 
motions  to  be  geometrically  similar  with  respect  to  the  remote  velocity. 

This  similarity  is  expected  to  simplify  the  diagnosis  of  any  shortcomings 
in  the  analysis. 

The  specific  conditions  chosen  for  correlation  are  given  in  Table  XXVI. 

The  condition  numbers  correspond  to  the  experimental  condition  numbers 
listed  in  Tables  XXII  through  XXV.  The  various  unsteady  airfoil  data 
scaling  and  inflow  calculations  are  discussed  under  the  subheadings  below. 

Unsteady  Airfoil  Data  and  Associated  Scaling  Procedure 

Aerodynamic  loadings  for  the  blades  are  calculated  in  the  normal  modes  aero- 
elastic  analysis  by  using  a  table  look-up  procedure  at  each  local  segment 
and  time  step  of  the  numerical  integration.  The  effects  of  unsteady  flow 
are  considered  by  the  method  presented  in  Reference  15.  This  method  uses 
generalized  data  from  a  sinusoidally  oscillating  airfoil  test.  The  normal 
force  and  pitching  coefficients  are  made  functions  of  angle  of  attack,  a, 
nondimensional  time  rate  of  change  of  angle  of  attack.  A,  and  nondimens ional 
pitch  acceleration,  B.  The  parameters  A  and  B,  as  they  are  used  in  the 
equation  of  motion,  are  defined  by  the  following  equations: 
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B  = 


The  unsteady  aerodynamic  tables  used  herein  vere  versions  of  the  data  used 
in  Reference  15,  with  appropriate  scaling  factors  applied  to  relate  the 
data  to  variois  Mach  numbers  and  model  scale.  Figures  47(a)  and  47(b)  are 
plotted  from  the  unsealed  tables  and  show  the  effects  of  typical  values  of 
the  parameters  A  and  B  for  an  NACA  0012  airfoil  chord  of  2  feet  and  a  Mach 
number  of  .20. 

The  versio  i  of  the  normal  modes  aeroelastic  analysis  used  herein  utilized 
steady-stete  aerodynamic  coefficients  for  drag  and  for  local  blade  condi¬ 
tions  outside  the  angle-of-attack  range  of  the  unsteady  aerodynamic  co¬ 
efficient  tables. 

As  discussed  in  References  15  and  20,  provision  has  been  made  in  the  normal 
modes  analysis  to  scale  the  single  set  of  unsteady  aerodynamic  normal  force 
and  pitching  moment  coefficients  to  various  other  Mach  numbers. 

The  local  value  of  the  unsteady  normal  force  coefficient  Cn(M,  a,  A,  B)  on 
a  blade  segment  at  Mach  number  M  is  found  from  the  tabulated  data 
Cnt(<>T»  A,B)  by  using  the  following  transformation: 

Cnoa(^*)  dent 

The  local  value  of  the  unsteady  moment  coefficient  C  (M,a,  A,  B)  on  a 
blade  segment  at  Mach  number  M  is  found  from  the  tabulated  data  C^fcup.AjB) 
by  using  the  following  transformation: 

CB(M,a,A,B)  *  Cntfc^  =  a,A,B) 

sm'  ' 

where  asn(M)  80,1  asin'M)  8X6  the  r*ap«ctive  angles  of  attack  for  normal  force 
and  pitching  moment  stall  at  Mach  number  M,  and  where  aBnt  and  a8nrt,  cure  the 
similar  quantities  for  the  tabulated  data  at  A»B»0. 

The  original  Intent  of  these  seal  ng  procedures  vas  the  transformation  of 
the  one  set  of  A«B*0  data  for  M*.l  and  c«2ft.,  to  match  other  steady-state 
tvo-dimensional  data  as  closely  as  possible.  It  is  assumed  that  this 
seeding  procedure  also  applies  satisfactorily  to  the  tabulated  unsteady 
data  for  which  Aj<0  and  Bj*0.  In  Reference  15  the  scaling  procedure  was 
developed  and  used  to  account  for  the  effects  of  varying  Mach  number. 

In  the  present  model  case,  the  scaling  procedures  attempted  to  account  for 
a  large  change  in  Reynolds  number  as  well  as  Mach  number.  The  use  of  the 
scaling  transformations  in  this  way  for  the  model  is  considered  to  be  an 
approximation,  and  reflects  the  lack  of  an  unsteady  aerodynamic  coefficient 
tabulation  at  model  scale  for  the  normal  modes  aeroelastic  analysis  Y-200 
program. 

The  transformations  expressed  in  the  above  equations  vere  used  to  scale  the 
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existing  unsteady  data  table  for  use  with  the  model.  The  scaling  para¬ 
meter  inputs  for  the  normal  modes  analysis  Y-200  program  are  given  in  Table 
XXVII.  The  special  low  Mach  number  values  in  Table  XXVII  allov  the  desired 
transformation  to  be  carried  out  with  the  existing  computer  program.  The 
Y-200  computer  program  identifies  values  given  for  M*0  in  the  input  tables 
of  asn(M)  and  a^CM)  versus  Mach  number  as  the  quantities  asnt  and  nSBrt. 
Hence,  machine  interpolation  in  those  tables  will  result  in  the  desired 
values  of  cxgnM^snt  and  aam(M)/a8jnt  for  all  significant  local  velocities. 
The  ratio  Cn88(M)/Cn88.,  is  input  directly  as  a  function  of  Mach  number, 
so  no  special  treatment  is  required. 

Because  the  unsteady  aerodynamic  data  transformation  to  model  scale  is 
approximate,  scaling  parameters  were  considered  as  variables  in  the  normal 
modes  analysis  calculations,  and  were  changed  as  necessary  to  provide  more 
useful  comparisons  between  analytical  results  and  the  test  data. 

Comparison  of  initial  normal  modes  aercelastic  analysis  calculations  for 
the  H-31*  model  blades  with  test  results  from  the  Reference  5  and  21  model 
tests  indicated  that  a  straightforward  application  of  the  scaling  method 
resulted  in  rotor  performance  predictions  that  were  far  too  pessimistic. 
Rotor  performance  parameters  from  the  present  test  later  confirmed  this 
observation.  This  behavior  was  aggravated  by  calculated  angle-of- attack 
ranges  that  exceeded  the  tabulated  unsteady  data  angle- of- attack  ranges. 

In  order  to  obtain  more  useful  results  with  the  existing  Y-200  computer 
program  and  unsteady  data  tables,  three  different  unsteady  aerodynamic 
data  scaling  cases  were  utilized.  The  parameter  of  these  scaling  cases  are 
given  in  Table  XXVII.  The  use  of  each  scaling  cast  is  shown  in  Table  XXVI. 
The  results  of  applying  each  scaling  case  to  the  tabulated  data  are  shown 
in  Figures  U8  through  U9.  The  aerodynamic  coefficients  for  a  2-foot  chord 
and  a  Mach  number  of  .2,  and  for  zero  values  of  the  parameters  A  and  B  are 
plotted  on  these  figures  from  the  data  tables  for  tv.e  Y-200  program.  Also 
shown  in  Figures  U8  through  1»9  are  the  results  of  applying  the  three 
different  scaling  cases  to  the  tabulated  data  for  A»B«0,  and  the  corres¬ 
ponding  two-dimensional  steady-state  test  data. 

Scaling  cases  1  and  2  were  set  up  and  utilized  to  avoid  excessively  pessi¬ 
mistic  performance  calculations,  and  to  stay  within  the  angle- of- attack 
range  of  the  existing  unsteady  data  tables.  The  transformed  A*B*0  coeffi¬ 
cient  of  lift  curves  for  these  cases  compare  reasonably  well  with  previous 
steady-state  data  synthesized  from  model  hover  tests,  and  presented  in 
Figure  11  of  Reference  21.  Scaling  case  3  was  developed  to  match  the  trans¬ 
formed  A«B»0  data  to  the  static  two-dimensional  data  reasonably  well  with 
the  existing  procedures.  This  case  was  used  to  examine  the  effects  of  non- 
uniform  inflow  on  the  pessimistic  performance  predictions  mentioned  earlier. 

Koch  Number  and  A  axle-  o  f-  At  tack  Boundaries  for  Unsteady  Aerodynamic  Utili¬ 
zation 


Provision  has  been  made  in  the  Y-200  normal  nodes  aeroelastic  analysis  com¬ 
puter  program  for  user  control  of  the  Mach  number  and  unsealed  angle- of- 
attack  (a)  region  over  which  unsteady  aerodynamic  coefficients  are  utilized. 
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This  control  is,  of  course,  limited  by  the  angle-of-attack  range  for  vhich 
the  unsteady  aerodynamic  data  are  tabulated.  Outside  the  user  defined  Mach 
number  and  angle-of-attack  region,  the  computer  program  automatically  re¬ 
verts  to  steady  flow  aerodynamic  data  tables.  A  transition  region  is  also 
user  defined,  vithin  vhich  the  computer  program  interpolates  between  the 
unsteady  and  steady  flow  aerodynamic  coefficients. 

The  above  described  unsteady  aerodynamic  utilization  boundaries  for  the 
normal  modes  aeroelastic  analysis  calculations  described  herein  are  given 
in  Table  XXVIII.  If  a  blade  segment  angle  of  attack  at  a  given  local  Mach 
number  is  less  than  the  values  on  the  boundary  defined  by  CNREG1  versus  Mach 
number,  the  local  blade  normal  force  coefficient  is  taken  from  the  table  of 
unsteady  aerodynamic  coefficients  by  the  computer  program.  If  the  angle  of 
attack  is  above  CNREG2,  the  normal  force  coefficient  is  taken  from  the  table 
of  steady  aerodynamics.  If  the  angle  of  attack  lies  betveen  the  two 
boundaries,  an  interpolation  betveen  steady  and  unsteady  data  on  the  re¬ 
spective  boundaries  is  used.  In  order  to  perform  this  interpolation,  co¬ 
efficient  values  are  found  at  the  Mach  number  on  both  the  CNREG1  (unsteady) 
and  the  CNREG2( steady )  boundaries.  The  coefficient  to  be  used  is  then 
found  ty  interpolation  with  respect  to  angle  of  attack.  The  scaling 
techniques  discussed  under  the  previous  subheading  are,  of  course,  utilized 
when  finding  the  unsteady  value  on  CNREG1.  The  coefficient  of  moment  data  is 
handled  in  an  analogous  manner,  using  the  CMREG1  and  CMREG2  boundaries 
also  given  in  Table  XXVIII.  Plots  of  boundaries  of  the  type  discussed  here 
are  given  in  Figure  31  of  Reference  15. 

The  unsteady  data  utilization  boundaries  for  the  flapped  (6p-5°)  blade  were 
specially  set  to  calculate  the  response  of  that  blade  with  the  existing  com¬ 
puter  program  and  set  of  symnetrlcal  unsteady  airfoil  date.  The  calcula¬ 
tions  for  the  flapped  airfoil  used  the  unsymnetrical  steady-state  data  when¬ 
ever  the  local  blade  Mach  number  vas  greater  than  ^.28.  This  utilization 
boundary  for  the  flapped  blade  is  given  in  Table  XXVIII.  Use  of  this 
boundary  allowed  the  effects  of  the  unstalled  coefficient  of  moment  and 
lift  due  to  camber  to  be  considered  in  the  calculation  with  the  existing 
computer  program  and  symmetrical  unsteady  data  tables.  On  the  retreating 
blade  azimuth,  unsteady  aerodynamics  with  its  effects  on  stall  were  con¬ 
sidered.  Figure  50  shows  a  sample  of  steady  and  unsteady  data  utiliza¬ 
tion  regions  for  the  flapped  blade.  This  use  of  steady  and  unsteady  data 
vas  expected  to  produce  a  reasonably  good  approximation  of  the  proper  air¬ 
loads,  since  the  area  ofosteady-state  stall  for  the  blade  extended  from 
about  t  *  190°  to  ^  *  20°.  Above  steady  state  stall  angles  of  attack, 
comparisons  of  Figures  lMa)  and  15(a)  show  the  flapped  and  unflapped 
data  to  be  similar,  at  least  within  the  accuracy  of  the  data  scaling 
techniques  discussed  under  the  previous  subheading. 

Rotor  Inflow  Velocity  Calculations 

The  basic  inflow  velocity  option  for  the  normal  modes  aeroelastic  analysis 
Y-200  computer  program  is  a  uniform  distribution  vhose  nondimens ionali zed 
value  (X)  is  one  of  the  input  parameters.  Upon  completion  of  a  calculation 
of  blade  response  and  rotor  performance  ,  the  value  of  shaft  angle  (ag)  is 
ob tinned  by  the  computer  program  from  the  following  momentun  relation: 


as  »  ARCTAnT-  -  ^L__l 

8  [v  2pV^2j 

In  the  present  case,  the  value  of  as  was  considered  to  be  an  independent 
variable.  The  input  values  of  A  consistent  with  the  desired  values  of 
C^,  a8,  and  y  were  obtained  from  the  following: 

A  -  y  sin  a8  -  Si 

As  an  option,  the  normal  modes  aeroelastic  analysis  Y-200  computer  program 
can  consider  input  values  of  rotor  inflow  velocity  which  vary  over  the  rotor 
disc.  The  variable  inflow  velocities  used  herein  were  considered  to  be  due 
to  blade  vortex  wakes  whose  geometry  was  unaffected  by  their  own  induced 
velocity.  This  type  of  wake,  which  can  be  described  au  a  prescribed  or  un¬ 
distorted  wake  is  formed  by  vortex  elements  which  travel  downstream  from 
the  rotor  with  constant  relative  velocity. 

The  transport  velocity  of  the  prescribed  wake  elements  is  assumed  to  be 
made  up  of  a  component  V  in  the  free-stream  relative  velocity  direction 
and  a  momentum  component 


downwards  normal  to  it. 

Hie  inflow  velocities  in  the  rotor  disc  are  calculated  by  a  method  similar 
to  that  used  in  Reference  17.  In  the  method  used  herein,  computational 
efficiency  is  inproved  by  considering  only  the  trailing  vortex  effects  in 
the  computation  of  inflow  velocity.  Shed  vortex  effects  are  assumed  to  be 
sufficiently  well  represented  by  the  use  of  unsteady  aerodynamic  data  in 
the  blade  response  calculations. 

Figure  51  is  a  simplified  diagram  of  one  of  the  blades  of  a  four-bladed 
rotor  and  the  undistorted  tip  vortex  segments  in  proximity  to  it  as  it 
passes  through  one  revolution.  The  view  is  downward  normal  to  the  rotor 
disc.  The  rotor  hub  is  passing  the  observer  from  left  to  right  with  con¬ 
stant  velocity  V.  The  straight  line  segments  represent  the  blade  at  various 
azimuth  positions.  The  curved  lines  represent  tip  vortices  laid  down  by 
the  blade  tips.  The  vortex  elements  are  moving  normal  to  the  plane  of  the 
paper.  The  induced  velocity  due  to  the  vortex  segments  is  upward  on  the 
convex  side  of  each  sequent  and  downward  on  the  concave  side. 

Samples  of  the  azimuthal  variation  in  inflow  velocity  due  to  the  wake  shown 
in  Figures  51  are  given  in  Figure  52  for  various  constant  radii.  The 
various  positive  and  negative  peaks  in  the  calculated  distribution  do 
correlate,  at  least  approximately  with  the  radial  and  azimuthal  positions 
of  some  of  the  blade-vortex  intersections  on  Figure  51.  The  effects  cf 
variations  in  vortex  strength  with  cyclic  blade  loading,  and  the  effects 
of  the  vertical  position  of  the  vortex  segments  are  also  considered  in  the 
calculation  of  the  induced  velocity  distribution.  Vortex  segments  for  all 


four  blades  from  four  revolutions  are  represented  in  the  induced  velocity 
calculation  as  used  herein. 

Presentation  and  Discussion  of  Comparative  Test  and  Calculated  Data 

Model  rotor  test  data  and  corresponding  data  calculated  vith  the  normal 
modes  aeroelastic  analysis  Y-200  computer  programs  are  compared  herein, 
with  the  objectives  of  finding  out  how  well  the'  analysis  works,  and  of  pro¬ 
viding  information  and  motivation  for  diagnosis  of  shortcomings.  The  extent 
of  the  comparisons  and  diagnostic  data  presented  herein  is  necessarily 
limited,  but  is  sufficient  to  support  an  assessment  of  the  present  status 
of  the  Y-200  computer  program  and  some  recommendations  for  its  improvement. 

The  comparisons  between  calculated  and  test  results  for  rotor  performance 
and  blade  response  should  be  extended  in  future  work  to  include  test  data 
for  a  larger  sample  of  the  operating  conditions  obtained  under  this  program. 
Future  work  should  also  relate  the  effects  of  a  wider  range  of  analytical 
assumptions  and  options  to  the  correlation  between  calculated  and  test  data. 

Comparisons  herein  between  theory  and  test  are  first  carried  out  on  the 
basis  of  the  more  important  rotor  performance  parameters  and  selected  one- 
half  peak-to-peak  blade  bending  moments  for  a  relatively  quick  overview  of 
the  quality  of  correlation  between  analysis  and  test.  These  comparisons 
appear  in  Figures  53  through  58,  and  are  sumnarized  further  in  Tables  XXIX 
through  XXXIII. 

The  simplified  comparisons  of  Figures  53  through  56  are  necessary  for  a 
comprehensive  Judgment  of  the  practiced  value  of  the  analysis,  but  are  not 
sufficient  to  fully  evaluate  a  highly  detailed  mathematical  model  such  as 
the  normal  modes  aeroelastic  analysis. 

In  contrast  to  more  empirical  methods,  the  normal  modes  aeroelastic  analy¬ 
sis  predicts  the  desired  rotor  design  information  by  using  a  highly  de¬ 
tailed  mathematical  simulation  of  the  rotor  blade  and  its  environment  to 
calculate  a  response  time  history.  If  the  simulation  is  successful,  the 
calculated  and  experimental  time  history  will  be  more  or  less  similar,  and 
it  can  then  be  concluded  that  the  analysis  is  producing  a  reliable  predic¬ 
tion  through  its  detailed  consideration  of  many  interacting  factors. 

In  view  of  the  above  considerations,  sample  comparisons  of  calculated  and 
average  revolution  test  data  are  provided  for  each  analytical  condition  in 
Figures  59  through  89.  In  addition  to  their  use  as  a  tool  for  evaluation 
of  the  analysis,  the  time  history  comparisons  are  required  as  a  guide  for 
detailed  diagnosis  of  shortcomings. 

Overview  of  Test  and  Calculated  Data  Correlation 

In  order  to  provide  a  useful  comparison  between  test  and  calculation,  the 
normal  modes  aeroelastic  analysis  was  trimmed  to  essentially  the  same  lift, 
shaft  angle,  first-harmonic  flapping  components,  velocity,  and  rotational 
speed  as  the  corresponding  test  condition.  Choice  of  these  parameters  as 
independent  variables  prevented  the  test  and  analysis  correlation  from 
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being  affected  by  different  dynamic  loadings  due  to  different  amounts  of 
Made  coning  and  flapping,  and  provided  a  similar  geometrical  relationship 
between  the  rotor  and  its  wake  for  the  variable  inflow  calculation. 

The  correlation  between  calculation  and  test  at  constant  lift  on  the  drag 
versus  torque  plots  for  each  rotor  are  shown  in  Figures  21(a)  through  21(c). 
The  calculated  and  test  data  for  rotor  performance  parameters  and  cyclic 
pitch  required  for  trim  versus  collective  pitch  is  presented  in  Figures  29 
through  31.  It  should  be  noted  that  the  uniform  inflow  calculated  data 
presented  in  Figure  28(a)  does  not  imply  a  continued  predicted  linear  trend 
for  rotor  lift  coefficient  versus  collective  pitch.  Additional  uniform 
inflow  cases  would  show  a  decreased  slope  for  higher  collective  pitch,  as 
indicated  by  the  variable  inflow  calculated  result  shown. 

Figures  53  through  58  and  Tables  XXIX  through  XXXIII  have  been  prepared  to 
provide  an  overview  of  blade  bending  moment  correlation  and  a  more  intensive 
examination  of  performance  correlation  between  calculation  and  test.  Separ¬ 
ate  comparisons  between  theory  and  test  values  of  rotor  drag  coefficient, 
torque  coefficient,  and  control  position  are  made  at  essentially  the  same 
lift,  shaft  angle,  and  first-harmonic  flapping. 

The  comparisons  between  these  quantities  may  be  of  secondary  interest  when 
comparing  different  rotor  systems ,  but  are  more  important  for  comparisons 
between  theory  and  test.  The  solutions  obtained  from  a  highly  detailed 
mathematical  model  like  the  normal  modes  aeroelastic  analysis  are  the 
result  of  complicated  interactions  between  many  factors.  In  order  to  have 
confidence  in  the  mathematical  simulation,  it  must  be  verified  that  quanti¬ 
ties  such  as  control  position  and  shaft  angle  are  approximately  correct. 

If  a  need  for  improvement  of  the  analysis  does  become  evident,  comparisons 
of  such  quantities  are  required  to  diagnose  the  source  of  the  inaccuracy. 

On  Figures  53  through  57,  the  four  blade  configurations  are  identified  on 
the  horizontal  axis  of  the  plots,  in  terms  of  the  parameters  varied,  which 
are  blade  twist,  flap  deflection,  and  elastic  stiffness.  Figures  53  through 
57  show  the  effects  of  blade  configuration  properties  on  rotor  performance, 
control  position,  and  one-half  peak-to-peak  blade  bending  moment  amplitude. 
In  addition,  the  effects  of  variable  inflow  are  presented  on  Figures  5*»  and 
55,  and  the  effects  of  unsteady  aerodynamic  scaling  assumptions  are  shown 
on  Figure  5U.  Figure  58  shows  the  results  of  a  further  investigation  of 
the  effects  of  variable  inflow  for  the  zero  degree  twist  dynamically  scaled 
fiber  glass  spar  blades.  The  horizontal  scale  has  been  extended  on  this 
plot  to  provide  more  clarity. 

The  quantitative  data  on  the  correlation  between  test  and  analysis  has  been 
further  summarized  in  Tables  XXIX  through  XXXIII.  It  should  be  noted  that 
this  further  sumnarization  is  meaningful  only  for  the  few  particular  opera¬ 
ting  cases  studied  here. 

Table  XXIX  presents  the  average,  range,  and  (where  applicable)  the  root- 
mean-square  of  the  difference  between  calculated  and  test  values  of  rotor 
drag  and  torque  coefficients  for  the  conditions  analyzed.  The  significance 
of  the  differences  between  calculated  and  test  values  can  be  evaluated  by 
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comparison  to  typical  operating  values  such  as  -.010  for  Cp/a  and  .006  for 
Cq/o. 

Table  XXX  presents  the  average  and  range  of  the  fraction  of  test  vibratory 
moment  predicted  by  the  analysis ,  for  the  conditions  analyzed. 

Tables  XXXI  and  XXXII  present  a  summary  comparison  of  calculated  and  test 
increments  in  rotor  drag  and  torque  coefficients  which  are  due  to  configura¬ 
tion  changes. 

Table  XXXIII  presents  a  summary  comparison  of  calculated  and  test  values  of 
bending  moment  amplitude  for  the  various  configurations,  divided  by  the 
corresponding  amplitudes  for  the  zero-degree  twist,  zero-degree  flap, 
dynamically  scaled  blade. 

The  information  contained  in  Figures  53  through  58  and  21  supports  the 
following  observations: 

Uniform  Inflow  Calculations 

The  magnitude  and  sign  of  deviations  in  rotor  drag  and  torque  between  cal¬ 
culation  and  test  for  the  uncambered  blades  were  inconsistent  but  compensa¬ 
ting.  Predictions  of  the  effects  of  configuration  change  on  rotor  perfor¬ 
mance  were  reasonably  good  for  moderately  loaded  conditions,  but  became  in¬ 
consistent  for  heavily  stalled  conditions.  Performance  prediction  for  the 
cambered  blade  was  poor,  as  shown  in  Figure  21. 

Calculation  of  rotor  pitching  moment  for  the  moderately  loaded  out-of¬ 
flapping-trim  condition  was  reasonably  good,  except  for  an  over-prediction 
on  the  three-times-scale  stiffness  aluminum  spar  blade. 

Prediction  of  control  positions  required  for  the  uncambered  blades  was 
quite  accurate,  except  for  the  heavily  stalled  operating  conditions.  The 
changes  in  control  position  required  due  to  blade  camber  were  poorly  pre¬ 
dicted,  even  at  the  lightly  loaded  conditions. 

The  correlation  of  cadculated  and  test  bending  moment  was  notably  better  on 
the  -8°  twist  blade  than  on  the  other  blade  configurations.  On  the  un¬ 
twisted  blades,  the  use  of  uniform  inflow  resulted  in  a  substantial  under¬ 
prediction  of  flapwise  and  chordwise  bending  moment  amplitude. 

The  blade  torsion  moment  analytical  prediction  was  inconsistent  in  its 
agreement  with  the  comparatively  low  test  values  existing  at  the  moderately 
loaded  conditions.  There  was  a  tendency  to  overpredict  the  higher  torsion 
moments  existing  at  the  more  heavily  stalled  operating  conditions. 

The  qualitative  effects  of  blade  configuration  change  on  blade  flapwise 
and  torsion  moments  were  generally  predicted  by  the  analysis.  The  quanti¬ 
tative  effects  of  blade  twist  were  substantially  overpredicted. 

With  the  uniform  inflow  assumption,  the  effects  of  configuration  change  on 
chordwise  bending  moment  amplitude  were  poorly  predicted. 
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Variable  Inflow  Calculations  and  Aerodynamic  Scaling 


The  variable  inflow  calculations  utilized  a  prescribed  geometry  vortex 
wake,  as  previously  described  on  pages  1*0  and  1*1.  The  effects  of  unsteady 
aerodynamic  scaling  assumptions  are  presented  in  combination  with  variable 
inflow.  These  scaling  assumptions  have  been  described  on  page  36. 

Scaling  assumption  3  matched  the  scaled  unsteady  aerodynamic  table  to  the 
static  aerodynamic  data.  This  resulted  in  a  pessimistic  performance  pre¬ 
diction  and  considerable  error  in  control  position,  as  shown  on  Figure  5**. 
This  was  true  even  though  variable  inflow  was  included. 

Use  of  scaling  assumption  1  in  conjunction  with  variable  inflow  produced 
the  results  shown  in  Figure  51*  for  the  -8-degree  twist  blade,  Figure  58  for 
the  0-degree  twist  blade,  and  Figure  55  for  all  four  configurations. 

The  results  shown  in  Figure  58  indicate  that  the  assumption  of  a  vortex 
core  of  size  .3c  is  better  than  .6c,  and  .3c  was  used  in  the  remaining 
calculations. 

The  use  of  variable  inflow  with  scaling  assumption  1  usually  resulted  in 
significant  improvements  in  performance  and  control  position  correlation 
over  that  obtained  with  the  uniform  inflow  assumption,  except  for  the 
cambered  blade.  The  rotor  pitching  moment  coefficient  correlation  obtained 
for  the  untrimned  flapping  condition  36  was  usually  more  exact  with  the 
uniform  inflow  calculation  than  the  variable  inflow  calculation.  An 
exception  to  this  trend  was  noted  for  the  aluminum  spar  three-times-scale 
stiffness  blade,  as  shown  on  Figure  51*. 

The  use  of  variable  inflow  in  the  calculation  improved  flapwlse  bending 
moment  correlation  over  that  obtained  with  uniform  inflow,  but  increased 
the  tendency  to  overpredict  torsional  moments.  Chordwise  moment  correla¬ 
tion  was  improved  by  variable  inflow,  but  the  magnitude  was  still  generally 
underpredicted. 

Discussion  of  Correlation  Between  Test  and  Calculated  Time  History  Data 

Sample  comparisons  of  average  revolution  time  histories  and  normal  modes 
aeroelastlc  analysis  results  are  provided  in  Figures  59  through  89.  Exami¬ 
nation  of  these  figures  supports  the  discussions  to  follow  which  supplement 
those  given  in  the  previous  subsection.  The  independent  parameters  of  the 
various  operating  conditions  used  for  the  correlation  study  are  given  in 
Table  XXVI.  A  tabulation  of  all  comparative  test  and  calculated  perfor¬ 
mance  parameters  is  given  in  Tables  XXXIV  through  XXXIX . 

The  discussions  and  observations  which  follow  concern  the  effects  on  cor¬ 
relation  between  calculation  and  test  of  the  inflow  and  aerodynamic  scaling 
assumptions  and  the  parameters  of  the  four  blade  configurations. 

Flapwlse  Bonding  Response 

Flapwise  bending  time  histories  calculated  for  the  zero  twist  blades  have 
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some  features  which  are  similar  to  the  test  data,  but  vhich  are  substanially 
underpredicted  with  the  uniform  inflow  calculations.  For  examples  of  this 
see  Figures  59,  60,  63,  64,  65,  66  and  67.  The  underprediction  of  the 
relative  peaks  at  azimuth  angles  of  about  60  degrees  and  300  degrees  on 
the  outboard  part  of  the  blade  is  a  particularly  important  shortcoming. 

The  use  of  variable  inflow  in  the  calculation  generally,  but  not  always, 
improves  flapwise  bending  correlation  for  the  untwisted  blades,  but  ampli¬ 
tudes  are  still  generally  underpredicted.  Examples  of  flapwise  bending 
correlation  with  variable  inflow  are  given  in  Figures  60,  64,  and  67. 

The  correlation  of  flapwise  bending  calculation  and  test  time  history  is 
generally  much  better  for  the  -8-degree  twist  blade  than  it  is  for  the  un¬ 
twisted  blades.  This  can  be  seen  by  comparing  Figures  59  and  6l  for  the 
uniform  inflow  assumption  and  Figures  60  and  62  for  the  variable  inflow 
assumption.  As  with  the  untwisted  blades,  the  use  of  variable  inflow  in¬ 
stead  of  uniform  inflow  sometimes  degrades  the  correlation. 

The  use  of  unsteady  aerodynamic  scaling  case  1  instead  of  scaling  case  3 
improves  flapwise  bending  correlation  by  a  relatively  small  amount.  (A 
discussion  of  these  scaling  cases  is  provided  on  page  36. )  The  effect  of 
the  scaling  assumptions  on  flapwise  bending  correlation  can  be  seen  by  com¬ 
paring  Figures  62(a)  and  62(b). 

The  use  of  an  assumed  vortex  core  size  of  .6  chord  in  place  of  a  .3  chord 
core  size  reduces  the  overall  amplitude  of  the  already  underpredicted 
calculated  flapwise  bending  response,  as  shown  in  Figures  67 Cb)  and  67(c). 

It  also  reduces  the  higher  frequency  parts  of  the  calculated  response 
considerably. 

Chordwise  Bending  Response 

The  calculated  chordwise  bending  response  is  substantially  underpredicted 
for  the  untwisted  blades  with  the  uniform  inflow  assumption.  Some  of  the 
low  frequency  part  of  the  inboard  blade  chordwise  bending  moment  is  pre¬ 
dicted,  but  very  little  of  the  higher  frequency  response.  Examples  of  this 
are  given  in  Figures  68(a)  and  68(e). 

The  use  of  variable  inflow  in  the  calculation  improves  the  chordwise  bend¬ 
ing  amplitude  correlation  with  test  data.  The  time  history  waveform  corre¬ 
lation,  however,  is  sometimes  enhanced,  and  sometimes  degraded.  Examples 
of  the  effect  of  variable  inflow  can  be  seen  by  comparing  Figures  69(a) 
and  66(g)  or  69(b)  and  68(i). 

The  correlation  between  calculated  and  test  chordwise  bending  response  is 
generally  better  for  the  -8-degree  twist  blade  than  for  the  untwisted  blade 
configurations.  The  use  of  variable  inflow  usually  improves  correlation, 
but  sometimes  degrades  it.  Examples  of  the  effects  of  variable  inflow  on 
correlation  of  chordwise  bending  response  for  the  -8-degree  twist  blade  can 
be  seen  by  comparing  Figures  71(b),  and  71(c)  vith  Figures  70(b)  and  70(c) 
respectively. 

The  use  of  unsteady  aerodynamic  scaling  case  1  Instead  of  scaling  case  3 
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improves  chordwise  bending  correlation  by  a  relatively  small  amount.  (A 
discussion  of  the  scaling  cases  is  provided  on  page  36. }  The  effect  of  the 
scaling  assumptions  on  chordvise  bending  correlation  can  be  seen  by  com¬ 
paring  Figures  71(a)  and  71(b). 

The  use  of  an  assumed  vortex  core  size  of  .6  chord  in  place  of  a  .3  chord 
core  size  had  little  effect  on  the  calculated  chordvise  bending  response, 
as  shown  in  Figures  76(b)  and  76(c). 

Effect  of  Trailing-Edge  Flap  on  Bending  Response 

Figures  59  through  76  show  that  the  test  bending  response  for  the  flapped 
blade  configuration  resembles  the  response  for  the  -8-degree  tvist  configur¬ 
ation  more  than  it  resembles  the  response  for  the  other  untwisted  blades. 
This  is  true  for  both  the  flapvise  and  the  chordwise  response,  and  can  be 
attributed  to  an  aeroelastically  induced  tvist.  This  effect  of  the  trailing 
edge  flap  does  not  appear  in  the  calculated  response  plotted  in  the  Figures 
cited  above. 

Blade  Hinge  Motions 

Correlation  between  calculated  and  test  hinge  motions  is  considered  good 
for  most  of  the  sample  operating  conditions.  Examples  from  uniform  inflow 
calculations  are  given  in  Figures  77  and  79*  The  amplitude  and  waveform 
correlations  are  distinctly  better  than  those  for  the  bending  response. 

The  constant  offsets  in  the  lag  response  are  due  to  deviations  between  cal¬ 
culated  and  test  values  of  rotor  torque. 

The  use  of  the  variable  inflow  assumption  improves  the  hinge  motion  corre¬ 
lation  in  some  cases,  and  degrades  it  in  others.  Examples  of  the  effect  of 
the  variable  inflow  assumption  on  the  hinge  motion  correlation  are  given  in 
Figures  78  and  80. 

Torsional  Response 

The  use  of  the  scaled  unsteady  full-scale  aerodynamic  data  for  the  model 
calculations  generally  results  in  an  overprediction  of  retreating  blade 
torsional  response.  Examples  of  the  correlation  between  calculated  and  test 
torsional  response  are  shown  in  Figures  81,  83,  85  and  87  for  the  uniform 
inflow  assumption.  The  use  of  the  variable  inflow  assumption  aggravates 
the  tendency  to  overpredict  the  retreating  blade  torsional  response,  as 
shown  in  Fibres  82,  84,  86  and  88. 

Unsteady  aerodynamic  scaling  case  3  results  in  increased  blade  stalling, 
and  therefore  greater  overprediction  of  the  stall-induced  retreating  blade 
oscillations  than  with  scaling  case  1.  This  is  shovn  by  comparing  Figures 
83(a)  and  83(b). 

The  effects  of  variable  inflov  on  the  torsional  response  of  the  unstalled 
blade  are  shown  in  Figures  89(b)  and  89(c),  for  the  two  vortex  core  sizes. 
The  variable  inflov  distribution  used  resulted  in  an  overprediction  of 
torsional  response.  A  change  in  vortex  core  size  from  .3  chord  to  .6  chord 


had  little  effect  on  the  torsional  response  caused  by  variable  inflow. 

Under  moderate  loading  conditions,  test  data  shovs  that  the  retreating 
blade  torsional  oscillation  takes  place  at  a  frequency  which  is  lower  than 
the  torsional  natural  frequency.  This  is  shown  for  data  from  conditions  7 
and  36  in  Figures  8l  and  83 .  The  torsional  vibration  in  condition  7  is  at 
approximately  5  cycles  per  revolution,  which  is  considerably  lover  than  the 
torsional  natural  frequency  of  approximately  7.h  cycles  per  revolution  for 
the  dynamically  scaled  blades.  The  aluminum  spar  three-times-scale  stiff¬ 
ness  blade  torsional  response  at  conditions  7  and  36  is  shown  on  Figures 
87(a)  and  87(b).  The  torsional  response  for  this  blade  in  the  last  azimuth 
quadrant  exhibits  the  5-cycle-per- revolution  response,  and  another  response 
at  approximately  10  cycles  per  revolution.  Both  of  these  frequencies  are 
lower  than  the  torsional  natural  frequency,  which  is  approximately  13  cycles 
per  revolution  for  this  blade. 

The  torsional  oscillations  predicted  by  the  analysis  are,  in  contrast, 
higher  than  the  torsional  natural  frequency.  The  calculated  frequency  of 
these  oscillations  for  the  dynamically  scaled  blades  is  about  10  cycles 
per  revolution,  and  is  about  16  cycles  per  revolution  for  the  aluminum 
span:  three-times-scale  stiffness  blade.  These  frequencies  also  are  cal¬ 
culated  for  the  heavier  blade  loading  conditions.  At  the  heavier  loading 
conditions,  such  as  condition  71  or  condition  82,  the  retreating  blade 
test  torsional  oscillation  frequency  tends  to  be  close  to  the  torsional 
natural  frequency.  These  data  for  condition  71  are  presented  in  Figures 
8l(c),  83(c),  85(c),  and  87(c).  The  retreating  blade  torsional  response 
frequency  is  about  7.7  cycles  per  revolution  for  the  untwisted  dynamically 
scaled  blade,  vhich  compares  to  a  natural  frequency  of  7.1*  cycles  per  re¬ 
volution.  The  tvisted  dynamically  scaled  blaule  exhibits  a  response  fre¬ 
quency  of  about  6.9  cycles  per  revolution.  The  aluminum  spar  three-times- 
scale  stiffness  blaule  has  a  test  response  frequency  of  about  12  cycles 
per  revolution,  vhich  compares  with  the  torsional  natural  frequency  of  13 
cycles  per  revolution. 

The  correlation  between  calculated  and  test  torsional  response  for  the 
flapped  blade  is  shown  in  Figures  85  and  86.  The  torsion  data  for  the 
flapped  blade  has  a  positive  offset,  because  of  the  use  of  the  zero  lift 
hover  condition  as  a  reference  for  the  strain  gage  instrumentation.  As 
discussed  earlier  on  page  22,  this  blade  had  a  negative  pitching  moment 
at  the  zero  condition,  vhich  is  not  knuvn  accurately,  but  is  estimated  to 
be  about  5*0  inch-pounds. 

The  correlation  between  calculated  and  test  torsional  response  for  this 
blade  reflects  the  use  of  symmetrical  unsteady  data  on  the  retreating  part 
of  the  azimuth,  as  explained  on  page  39.  The  flapped  airfoil  test  data 
shovs  much  less  torsional  response  in  this  region  than  the  other  blades, 
because  of  the  higher  maximum  lift  coefficient  available  before  stalling 
occurs.  The  calculated  response  in  this  region  is  similar  to  that  for  the 
unflapped,  untwisted  blade,  since  the  same  data  is  being  used. 

The  variations  in  flapped  blade  torsional  response  vhich  occur  on  the  un¬ 
stalled  parts  of  the  azimuth  are  underpredicted  by  the  analysis. 


ASSESSMENT  OF  NORMAL  MODES  AEROELASTIC  ANALYSIS 


Previous  sections  of  this  report  present  and  discuss  the  correlation  be¬ 
tween  test  data  generated  under  the  present  program  and  calculated  data 
from  the  normal  modes  aeroelastic  analysis. 

These  comparisons  have  shown  that  satisfactory  predictions  of  rotor  per¬ 
formance  and  blade  response  are  possible  when  the  blades  have  a  conventional 
twist  rate  of  -8  degrees  and  a  synnetrical  airfoil  section. 

For  other  applications,  however,  the  comparisons  demonstrate  a  need  for  the 
improvement  of  the  calculation  method  itself,  or  the  associated  input  data. 
In  this  section,  evident  reasons  for  lack  of  satisfactory  correlation  are 
presented,  with  the  supporting  diagnostic  data.  Suggested  irprovements 
which  are  expected  to  benefit  the  accuracy  of  the  normal  modes  analysis 
are  also  presented. 

Evaluation  of  Unsteady  Aerodynamic  Data  and  Associated  Scaling  Procedures 

The  current  set  of  tabulated  unsteady  aerodynamic  data  for  the  Y-200  ver¬ 
sion  of  the  normal  modes  aeroelastic  analysis  was  obtained  from  oscillating 
NACA  0012  airfoil  tests  for  a  chord  of  2  feet  and  a  Mach  number  of  .2. 

These  unsteady  airfoil  data  are  applied  to  other  flow  conditions  by  multi¬ 
plying  the  tabulated  coefficients  and  angles  of  attack  by  scaling  factors, 
in  the  manner  explained  earlier  on  page  36. 

The  scaling  procedure  appears  to  work  reasonably  well  for  full-scale, 
symmetrical  blade  sections.  The  work  of  Reference  22  supports  this  con¬ 
tention.  Figure  15  of  Reference  22  shows  excellent  correlation  of  cal¬ 
culated  and  test  variations  of  C^/a  with  tip  path  plane  angle  of  attack. 
Except  for  a  constant  offset,  the  test  values  of  C^/a  beyond  stall  were 
obtained  at  the  test  values  of  collective  pitch  and  tip  path  plane  angle 
of  attack. 

The  present  investigation  used  dynamically  scaled  models  of  the  rotors  used 
for  the  Reference  22  correlation.  Figures  l»8(a)  through  U8(e)  show  the  full 
scale  (M».2,c»2  feet)  and  model  two-dimensional  static  coefficient  of  lift 
versus  angle  of  attack.  The  stall  coefficient  of  lift  and  angle  of  attack 
is  much  lower  for  the  model  airfoil.  Scaling  case  3  matched  the  unsteady 
aerodynamic  tables  to  the  static  datr  quite  well. 

Scaling  case  3  was  used  for  variable  inflow  calculations  sumarized  in 
Figure  51*.  The  excessively  high  collective  pitch  and  torque  coefficients 
calculated  to  reach  the  test  values  of  lift  and  flapping  indicated  that  the 
analysis  was  overpredicting  the  amount  of  blade  stall.  A  much  more  favorable 
correlation  was  obtained  for  the  twisted  blade  by  using  scaling  case  1, 
which  is  shown  in  Figure  U8.  The  improved  correlation  can  be  noted  on 
Figure  5** . 

These  results  imply  that  the  current  scaling  procedures  cannot  be  applied 
without  restriction.  The  errors  involved  in  scaling  to  other  full-scale 
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rotor  blade  sections  from  the  .2  Mach  number,  2-ft  chord  condition  will 
probably  be  much  smaller  than  the  substantial  discrepancy  noted  for  the 
model.  Nevertheless,  additional  unsteady  airfoil  data  should  be  obtained 
for  other  Mach  numbers ,  Reynolds  numbers ,  and  airfoil  sections ,  so  that  the 
effects  can  be  properly  evaluated. 

In  the  past  few  years,  new  rotor  blade  designs  have  been  provided  with 
cambered  airfoils.  Current  scaling  procedures  available  with  the  Y-200 
version  of  the  normal  modes  aeroelastic  analysis  cannot  adequately  repre¬ 
sent  the  combined  effects  of  a  cambered  airfoil  and  unsteady  aerodynamics, 
since  the  tabulated  unsteady  aerodynamic  data  is  for  a  symmetrical  airfoil. 

In  the  present  investigation  an  attempt  was  made  to  simulate  the  effects  of 
the  5-degree  trailing-edge  flap  airfoil  and  unsteady  aerodynamics  by  using 
the  flapped  airfoil  steady-state  aerodynamics  data  on  the  advancing  blade 
azimuth  and  unsteady  aerodynamic  data  on  the  retreating  blade  azimuth. 

This  technique  used  the  existing  Y-200  program  input  provisions,  as  des¬ 
cribed  on  page  38.  Useful  performance  predictions  could  not  be  obtained 
with  this  technique,  as  shown  in  Figures  53  through  57  and  21. 

Comparison  of  the  rotor  performance  versus  collective  pitch  data  for  the 
-8-degree  twist  blade  in  Figures  29(a)  through  29(d)  with  the  similar  data 
for  the  flapped  blade  in  Figures  31(a)  through  31(d)  indicates  that  the  sub* 
atantial  calculated  performance  discrepancy  for  the  flapped  blade  is  re¬ 
lated  to  the  underprediction  of  rotor  lift  at  a  given  collective  pitch. 

This  is  the  result  of  using  the  symmetrical  unsteady  data  on  the  retreat¬ 
ing  blade  azimuth.  This  data  does  not  reflect  the  higher  lift  coefficients 
available  with  the  flapped  airfoil.  The  requirement  for  lateral  trim  then, 
in  turn,  also  prevents  the  development  of  the  correct  calculated  lift  on 
the  advancing  blade  azimuth. 

Examination  of  the  various  unsteady  aerodynamic  scaling  cases  for  coeffi¬ 
cient  of  lift  and  moment  in  Figures  U8(a)  through  l*9(e)  indicates  a  limited 
potential  for  simulating  new  airfoil  characteristics  on  advanced  blade 
designs.  For  example,  an  abrupt  stall  of  the  type  shown  on  Figure  U8(c) 
cannot  be  represented.  The  moment  coefficient  scaling  procedure  consists 
of  an  angle-of-attack  factor  only,  which  provides  an  excessively  steep 
moment  stall,  compared  to  the  measured  steady  data  shown,  for  example,  on 
Figure  l»9(a).  The  coefficient  of  moment  below  stall  for  the  scaled  data 
at  zero  pitching  velocity  and  acceleration  is  essentially  zero,  and  there¬ 
fore,  the  steady  pitching  moment  coefficient  effects  of  camber  or  shifts  of 
aerodynamic  center  of  pressure  cannot  be  included. 

Further  generalization  of  the  unsteady  data  would  be  accomplished  by  sub¬ 
tracting  the  tabulated  zero  pitching  velocity  and  acceleration  data  from 
the  remaining  tabulated  data.  This  would  form  a  table  of  coefficient  in¬ 
crements  due  to  unsteady  aerodynamics.  This  table  could  be  modified  in 
some  prescribed  manner  to  reflect  Reynolds  number,  Mach  number,  and  other 
effects  as  they  are  defined  through  future  test  and  analysis. 

In  use,  the  proposed  table  of  unsteady  aerodynamic  increments  could  be 
added  to  any  particular  set  of  steady  aerodynamic  characteristics,  which 
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coild  also  reflect  various  types  of  stall  characteristics  and  the  effects 
of  blade  camber,  Mach  number ,  and  Reynolds  number . 

Further  improvements  in  the  mathematical  model  of  blade  airloads  have  been 
found  advisable  by  other  investigators  in  References  15,  20,  and  22.  These 
Improvements  include  computationally  simpler  analytical  definitions  of  un¬ 
steady  airfoil  behavior  in  place  of  the  current  tabular  data. 

Evaluation  of  Calculated  Blade  Bending  Response 
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Tht  quality  of  the  correlation  between  normal  modes  aeroelastic  analysis 
blaas  bending  response  calculated  results  and  test  data  reported  herein  is 
similar  to  what  has  been  published  previously. 

Most  of  the  previously  published  comparisons  have  displayed  results  for 
blades  with  a  conventional  amount  of  built-in  twist.  The  previous  full- 
scale  twisted  blade  correlation  is  similar  in  quality  to  that  reported 
herein  for  the  -8-degree  twist  model  blade.  Examples  of  blade  bending  cor¬ 
relation  obtained  earlier  for  twisted  blades  are  shown  in  Figures  25  through 
27  and  30  through  32  of  Reference  1  ,  Figures  65  and  66  of  Reference  15, 
Figures  58  and  6l  of  Reference  20,  and  Figure  53  of  Reference  23.  All  of 
these  twisted  blade  results  show  reasonably  good  correlation  between  cal¬ 
culated  and  test  flapwise  bending  time  history,  and  at  least  fair  amplitude 
correlation  for  chordwise  bending  response. 

Correlation  of  the  normal  modes  aeroelastic  analysis  results  with  test 
data  for  eui  untwisted  model  blade  was  published  in  Reference  5. 

The  work  of  Reference  5  was  concerned  primarily  with  transient  conditions, 
but  the  pretransient  steady-state  conditions  are  comparable  to  the  results 
of  the  present  investigation. 

The  results  in  Figure  19,  page  123,  of  Reference  5  show  a  calculated 
underprediction  of  flapwise  bending  response  for  the  uniform  inflow  assump¬ 
tion,  similar  to  the  quality  of  the  untwisted  blade  correlation  for  the 
present  investigation. 


The  chordwise  bending  response  is  also  substantially  underpredicted  by  the 
calculation  for  the  untwisted  blades,  as  shown  in  Figures  68,  69,  72,  73, 
7U ,  75,  and  76.  A  significant  improvement  in  flapwise  correlation  is 
obtained  by  including  variable  inflow  in  the  calculation  for  the  untwisted 
blades,  but  the  effect  on  the  chordwise  response  is  less  prominent. 

The  quality  of  the  correlation  between  calculated  and  test  flap  and  lag 
hinge  motion  data  is  better  than  the  bending  response  correlation,  and  is 
not  strongly  affected  by  blade  twist.  This  can  be  seen  by  comparing 
Figures  77(g)  and  77(h)  to  21(a)  and  21(b)  respectively. 

Figures  53  through  57  show  that  except  for  the  blade  with  the  trailing 
edge  flap,  the  correlation  of  the  calculated  control  positions  with  test 
data  for  the  untwisted  blades  is  about  as  good  as  it  is  for  the  twisted 
blades . 


50 


These  observations  about  the  correlation  of  calculated  and  test  data  provide 
motivation  for  more  detailed  investigations  into  the  reasons  for  the  lack 
of  bending  response  correlation.  The  aspects  of  the  calculation  method 
which  should  be  studied  further  in  this  regard  are  discussed  in  the  follow¬ 
ing  paragraphs. 

It  can  be  expected  that  the  untwisted  blade  is  more  strongly  affected  by 
nonuniform  Inflow,  and  that,  therefore,  inaccuracies  in  the  calculation  of 
its  inflow  are  more  important.  The  untwisted  blade  tip  is  more  highly 
loaded  in  the  lift  direction.  This  higher  tip  loading  will  cause  the  trail¬ 
ing  vorticity  to  be  more  concentrated  at  the  blade  tip.  This,  in  turn , 
causes  higher  tip  inflcw  for  the  untwisted  than  for  the  twisted  blade. 
Consideration  of  this  inflow  results  in  an  Induced  angle-of-attack  distri¬ 
bution  which  is  aerodynamically  equivalent  to  a  certain  amount  of  built-in 
twist.  This  results  in  load  distribution  components  which  tend  to  be  more 
like  those  of  the  twisted  blade.  Consideration  of  variable  inflow  in  the 
untwisted  blade  response  calculations  does  make  improvements  of  this  type 
in  the  correlation  with  test  data,  as  seen  by  comparing  Figures  67(n.), 

67(b),  and  67(c). 

Thus,  since  additional  improvement  of  the  same  type  of  flapwise  bending 
response  i r  still  needed,  more  accurate  calculations  of  wake-induced  inflow 
are  Indicated,  especially  for  the  untwisted  blade.  A  method  for  more 
accurate  calculation  of  the  blade  wake  geometry  has  been  utilized  and  re¬ 
commended,  for  example,  in  References  20  and  22. 

The  induced  inflow  velocity  calculation  utilizes  the  modified  form  of  the 
Blot-Savart  representation  for  a  large  number  of  straight-line  vortex  seg¬ 
ments  which  make  up  the  helical  wake.  This  vortex  theory  states  that  the 
induced  velocity  is  inversely  proportional  to  the  distance  from  the  fila¬ 
ment.  Without  modification,  this  Idealization  becomes  inaccurate  for 
sufficiently  small  distances  from  the  vortex  filament,  since  Indefinitely 
large  velocities  are  Implied.  Since  rotor  blades  often  pass  very  close 
to  the  wake  vortices,  indefinitely  large  inflow  velocity  impulses  would 
appear  in  the  calculation.  These,  in  turn,  would  cause  the  calculation  of 
unrealistically  large  load  impulses.  In  ozdier  to  avoid  this  theoretical 
shortcoming,  and  simulate  the  physical  vortices,  it  is  assumed  that  self- 
induced  velocities  are  zero  inside  of  an  arbitrarily  defined  core  radius 
around  each  vortex  segment. 

The  effects  of  two  different  core  size  assumptions  are  shown  in  Figures 
67(b),  76(b),  69(b),  and  67(c),  76(c),  89(c).  The  ass\nption  of  a  core 
radius  .6  times  blade  chord  produces  too  little  flapwise  blade  response. 

The  assumption  of  a  core  radius  of  .3  times  blade  chord  provides  more  lov^ 
frequency  response,  but  the  higher  frequency  part  of  flapwise  bending 
response  is  overpredicted.  These  results  show  that  for  this  operating 
condition  a  significant  part  of  the  calculated  blade  flapwise  response  is 
due  to  vortex  filaments  passing  closer  to  the  blade  than  .6  times  the  chord. 
Therefore,  the  accuracy  of  the  present  lifting  line  theory  with  nearby 
vortices  should  be  investigated,  and  a  more  refined  theory  used  if  it  is 
Justified. 
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The  quantitative  effects  of  limitations  on  the  number  of  modes  included  in 
the  normal  modes  aeroelastic  analysis  calculations  should  be  investigated. 

In  addition  to  rigid  blade  flapping  and  lagging  motions,  the  calculations 
use  a  limited  number  of  orthogonal  natural  vibration  modes  as  elastic 
degrees  of  freedom.  There  degrees  of  freedom  are  particularly  easy  to  use 
in  the  numerical  integration  of  the  equations  of  motion. 

As  discussed  in  Reference  11,  the  solution  of  the  equations  of  motion  is 
expressed  in  terms  of  a  summation  of  the  modal  motions.  Furthermore,  the 
elastic  bending  moments  or  stresses  are  calculated  from  a  sunn  at  ion  of 
modal  curvatures,  multiplied  by  the  appropriate  constants. 

This  formulation  of  the  problem  in  terms  of  orthogonal  modes  is  theoreti¬ 
cally  rigorous  if  enough  modes  are  included  in  the  solution.  dp  to  four 
elastic  flapvise  and  two  elastic  chordvise  modes  are  permitted  in  the 
current  version  of  the  Y-200  program. 

An  increase  in  the  number  of  modes  will  require  more  computer  time  to  per¬ 
form  the  calculations,  because  of  the  increase  in  the  number  of  equations  of 
motion,  and  because  smaller  azimuth  step  sizes  are  required  with  the  higher 
frequency  modes  to  be  added.  In  addition,  a  larger  number  of  spanvise  seg¬ 
ments  are  indicated,  so  that  the  more  complicated  higher  frequency  mode 
shapes,  slopes,  and  curvatures  can  be  adequately  represented. 

The  comparisons  between  test  and  calculated  data  suggest  that  some  loss  of 
accuracy  may  be  occurring  because  of  approximations  in  the  calculation  of 
blade  elastic  response.  As  pointed  out  previously,  the  correlation  of 
O-degree  twist  blade  hinge  motions  with  test  data  is  much  better  than  the 
corresponding  bending  moment  correlation.  This  suggests  that  blade  aero¬ 
dynamic  blade  loadings  are  calculated  with  reasonable  accuracy,  and  that  the 
elastic  response  to  this  loading  is  underpredicted.  The  improvement  in 
flapvise  correlation  between  the  0- degree  twist  and  the  -8-degree  twist 
blade  elastic  moments  may  be  due  to  greater  compatibility  between  the  -8- 
degree  blade  loading  snd  the  mode  shapes. 

Tables  XL  through  XLV  have  been  prepared  to  show  samples  of  the  contribu¬ 
tion  of  each  of  the  various  flapvise  and  chordvise  modes  to  the  total 
bending  moment.  In  each  case,  the  ccntribution  of  the  final  mode  added  is 
significant,  especially  for  the  chordvise  response.  This  implies  that 
additional  significant  response  could  be  obtained  by  including  additional 
modes. 

In  contrast,  the  modal  contributions  to  hinge  motion  become  progressively 
smaller  with  each  succeeding  mode,  as  shown  in  Table  XL VI . 

The  information  in  the  above  mentioned  tables  suggests  that  the  position  of 
the  blades  in  space  is  determined  with  comparative  accuracy  by  the  calcula¬ 
tion.  The  contributions  of  each  additional  mode  to  blade  position  does 
grow  very  small  after  the  first  few  modes.  This  is  in  contrast  to  the 
contribution  to  blade  curvature  of  each  additional  mode,  which  does  not 
become  negligible  after  the  first  few  modes. 
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In  view  of  the  above,  the  consideration  of  an  alternate  method  of  calcu¬ 
lating  blade  bending  moments  and  stresses  is  suggested.  The  bending  moment 
can  also  be  calculated  by  a  separate  spanvise  integration  of  the  already 
calculated  aerodynamic  and  inertial  loadings.  Table  XLVII  provides  a  sample 
summaticn  of  flapvise  modal  accelerations.  Even  though  the  component  local 
acceleration  and  corresponding  inertial  loading  contributions  do  not  become 
negligible  after  the  first  few  modes,  it  is  at  least  plausible  to  expect 
better  convergence  of  the  spanvise  integrals  which  form  the  shear  and 
bending  moments. 

Evaluation  of  Calculated  Blade  Torsional  Response 

The  correlation  of  calculated  torsional  response  from  the  normal  modes  aero- 
elastic  analysis  with  the  corresponding  test  data  shoved  that  the  calcula¬ 
tion  tended  to  overpredict  retreating  >  i  de  torsional  oscillations.  This 
effect  is  contrary  to  most  of  the  pre  -  .sly  published  similar  comparisons, 
which  typically  show  an  underprediction  of  torsional  response  by  the  cal¬ 
culation. 

Previous  comparisons  of  torsional  response  calculations  and  test  data  for 
full-scale  blades  appear  in  Figures  U8  through  53  and  58  through  62  of 
Reference  15,  and  Figures  50,  53  through  5 7,  59*  and  60  of  Reference  20. 
Previous  comparisons  of  torsional  response  calculations  and  Reference  5 
model  test  data  appear  in  Figures  62  through  6b  of  Reference  20. 

The  use  of  variable  inflow  tended  to  substantially  aggravate  the  tendency 
of  the  analysis  to  overpredict  retreating  blade  torsional  response.  This 
tendency  was  not  apparent  in  the  calculations  shown  in  Figure  30  of  Refer¬ 
ence  15.  Some  tendency  of  the  variable  inflow  to  cause  overprediction  of 
torsional  response  is  evident  in  a  comparison  between  Figures  50  and  59  of 
Reference  20. 

The  reasons  for  the  overprediction  of  torsional  response  by  the  calcula¬ 
tion  in  the  present  investigation  have  not  been  definitely  established. 

One  difference  between  the  present  set  of  comparisons  and  the  previous  ones 
is  the  higher  torsional  natural  frequency  to  rotational  frequency  ratio  of 
the  blades  used  in  this  Investigation.  Table  XLVI1X  furnishes  approximate 
sample  values  for  the  nondimensionalized  pitching  velocity  and  accelera¬ 
tion  A  and  B,  for  the  untwisted  dynamically  scaled  blade.  The  values 
given  for  the  variable  inflow  case  shown  exceed  the  unsteady  aerodynamic 
data  tabulated  range  of  A  *  i.Ob  and  B  *  t.01  by  a  small  amount  in  this 
case.  Nevertheless,  a  more  cosiprehensive  determination  of  how  often  the 
unsteady  aerodynamic  data  table  range  is  exceeded  should  be  carried  out. 

The  effect  on  the  blade  torsional  response  calculation  of  these  out-of- 
range  points  should  of  course,  also  be  determined. 

As  pointed  out  earlier,  the  analysis  also  has  a  tendency  to  overpredict 
the  retreating  blade  torsional  response  frequency.  At  the  stall  boundary 
the  experimental  data  indicates  a  moderately  large  blade  torsional  re¬ 
sponse  which  is  much  lover  than  the  torsional  natural  frequency.  This 
was  not  predicted  by  the  calculations  with  unsteady  data  utilized  herein, 
and  is  believed  to  be  caused  by  bending  torsion  coupling.  The  similar 


lower  frequency  response  appearing  in  Figure  30  of  Reference  22  for  steady 
aerodynamic  data  and  zero  moment  coefficient  supports  this  contention. 

At  the  higher  lift  conditions,  the  experimental  data  shows  the  expected 
retreating  blade  response  at  slightly  above  the  torsional  natural  frequency. 
The  calculation  predicts  a  response  in  the  same  region  with  a  considerably 
higher  frequency.  This  tendency  has  also  been  noted  in  the  previous  in¬ 
vestigations  of  Reference  15  and  20.  The  predicted  torsional  frequency  is 
about  35  percent  higher  than  the  experimental  value  for  the  fiber  glass 
spar  dynamically  scaled  blade  and  about  23  percent  higher  for  the  aluminum 
spar  three-times -sc ale  stiffness  blade.  These  frequenc;  ratios  imply  an 
effective  stiffness  82  percent  higher  than  the  elastic  r^iffness  for  the 
fiber  glass  spar  dynamically  scaled  blade  and  51  percent  higher  for  the 
aluminum  spar  three -times -scale  stiffness  blade.  Thus,  the  apparent  added 
stiffness  is  about  twice  as  high  in  absolute  terms  for  the  three-times- 
scale  stiffness  aluminum  spar  blade.  Since  the  apparent  stiffness  appears 
at  relatively  low  torsional  amplitude,  it  is  probably  not  a  result  of 
exceeding  the  unsteady  aerodynamic  table  range.  Further  investigation  is 
required  to  establish  the  source  of  the  apparent  torsional  stiffness. 

The  effect  of  variable  inflow  on  torsional  response  with  the  two  assumed 
vortex  core  radii  is  shown  in  Figure  89.  For  this  unstalled  case,  the 
retreating  blade  torsional  oscillations  occur  at  approximately  the  torsional 
natural  frequency.  The  current  theory  overpredicted  the  magnitude  of  the 
torsional  oscillations,  as  it  did  for  the  stalled  operating  conditions. 

The  low-frequency  torsional  response  of  the  trailing-edge  flapped  blade  was 
underpredicted  with  the  existing  unsteady  aerodynamic  tables,  while  the 
high-frequency  retreating  blade  torsional  response  was  overpredicted. 

This  was  primarily  due  to  the  lack  of  proper  unsteady  aerodynamic  data 
for  the  flapped  airfoil,  as  discussed  earlier. 
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CONCLUSIONS 
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1.  Sample  comparisons  of  model  and  full-scale  test  data  shov  qualitative 
and  quantitative  similarities.  Therefore,  it  is  concluded  that  analy¬ 
tical  methods  for  full-scale  rotor  performance  and  blade  response  can 
be  verified  or  improved  through  correlation  with  model  data. 

2.  The  blade  set  with  the  5-degree  trailing- edge  flap  is  observed  to  have 
substantially  better  performance  than  the  unflapped  blade  sets.  At  an 
advance  ratio  of  u  *  .3,  at  equal  lift  and  propulsive  force,  rotor 
torque  for  the  flapped  blade  set  is  typically  about  80  percent  of  the 
torque  for  the  blade  set  with  -8  degrees  of  built-in  twist. 

3.  Under  stalled  conditions,  the  aluminum  spar  three-times-scale  stiff¬ 
ness  blade  set  has  less  favorable  performance  than  the  other  blade 
sets.  At  an  advance  ratio  of  u  *  .3,  at  equal  lift  and  propulsive 
force,  rotor  torque  for  the  aluminum  spar  three-times-scale  stiffness 
blade  set  is  typically  about  110  percent  of  the  torque  for  the  similar 
fiber  glass  spar  dynamically  scaled  blade  set. 

b.  The  bending  moment  amplitudes  on  the  blade  set  with  the  5-degree  t rail¬ 
ing-edge  flap  are  similar  in  magnitude  to  those  on  the  set  with  -8 
degrees  of  built-in  twist. 

5.  The  torsional  moment  amplitudes  on  the  blade  set  with  the  5-degree 
tralling-edge  flap  are  much  higher  than  those  on  the  unflapped  blades, 
except  at  the  more  highly  loaded  condition*.  At  an  advance  ratio  of 

U  ■  .3,  the  torsional  moments  below  stall  on  the  flapped  blade  set  are 
typically  about  two  to  three  times  the  torsional  moment  on  the  unflapped 
blades.  Blade  s+'JLling,  however,  causes  the  torsional  moments  on  the 
unflapped  blades  to  be  higher  at  the  more  heavily  loaded  conditions. 

6.  The  bending  moment  amplitudes  on  the  aluminum  spar  three-times-scale 
stiffness  blades  are  substantially  higher  than  those  on  the  fiber  glass 
spar  scale  stiffness  blades.  At  an  advance  ratio  of  u  *  .3,  the 
bending  moments  on  these  blades  are  typically  about  two  to  three  times 
those  on  the  fiberglass  spar  scale  stiffness  blades. 

7.  Blade  torsion  moment  amplitudes  on  the  aluminum  spar  three-times-scale 
stiffness  blades  are  similar  to  those  on  the  other  blades  below  stall. 
Above  stall,  at  an  advance  ratio  of  u  ■  .3,  the  aluminum  spar  three- 
times-scale  stiffness  blade  torsional  moments  are  typically  about  25 
percent  to  50  percent  greater  than  those  on  the  fiber  glass  spar  scale 
stiffhess  blades. 

8.  It  appears  that  blade  torsional  vibration  can  be  aggravated  when  wake 
vortex  passage  frequency  relative  to  the  blade  is  close  to  the  torsional 
natural  frequency.  However,  this  effect  is  dependent  on  rotor  blade 
configuration  and  loading  condition. 

9.  In  order  to  obtain  reasonably  good  performance  correlation  between  the 
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model  rotor  test  data  and  the  normal  modes  aeroelastic  analysis,  the 
existing  full-scale  unsteady  aerodynamic  data  tables  must  be  scaled  to 
a  steady-state  maximum  lift  coefficient  which  is  higher  than  the 
measured  static  two-dimensional  model  blade  airfoil  data. 

10.  Correlation  between  calculated  rotor  performance  and  test  data  is  best 
for  the  -8-degree  twist  blade  with  the  variable  inflow  assumption. 

With  this  assumption,  deviations  in  drag  and  torque  are  on  the  order 
of  10  percent  of  typical  operating  values.  The  use  of  the  uniform 
inflow  assumption  results  in  typical  drag  and  torque  deviations  about 
twice  as  large.  The  deviations  in  drag  and  torque  at  a  given  condition 
are  compensating,  but  inconsistent  in  sign  for  different  conditions. 

11.  Correlation  between  calculated  rotor  performance  and  test  data  for  the 
untwisted  blades  shows  compensating  deviations  in  drag  and  torque  that 
are  on  the  order  of  20  to  30  percent  of  typical  operating  values.  The 
use  of  variable  inflow  in  the  calculations  does  not  consistently  improve 
the  accuracy  of  the  drag  and  torque  calculations. 

12.  Because  of  compensating  deviations  in  drag  and  torque,  correlation  of 
rotor  performance  with  test  data  on  a  drag  versus  torque  plot  at  con¬ 
stant  lift  is  good  for  all  three  blade  configurations  without  flaps. 
Deviations  in  torque  are  smaller  than  10  percent  when  correlation  is 
made  on  this  basis. 

13.  Increments  in  rotor  drag  and  torque  due  to  blade  configuration  are 
more  accurately  predicted  than  the  absolute  values.  For  the  unflapped 
blade  configurations,  the  average  increment  deviation  is  less  than  5 
percent  of  typical  operating  values  for  the  twist  variation  and  less 
than  20  percent  for  the  stiffhess  variation. 

lU.  Correlation  between  test  and  calculated  performance  for  the  blade  with 
the  tralling-edge  flap  is  poor,  because  of  overpredicted  rotor  torque 
at  a  given  lift.  This  is  attributed  to  the  use  of  the  existing  symme¬ 
trical  unsteady  airfoil  data.  These  data  do  not  reflect  the  combined 
beneficial  effects  of  unsteady  aerodynamics  and  camber. 

15.  Correlation  between  test  and  calculated  blade  vibratory  bending  moment 
amplitudes  is  reasonably  good  for  the  -8-degree  twist  blade  at  moderate 
loadings  when  the  variable  inflow  calculation  is  used.  Under  these 
conditions,  the  sample  flapwise  moment  amplitude  is  overpredicted  by 

6  percent,  the  chordwise  amplitude  is  overpredicted  by  15  percent,  and 
the  torsional  amplitude  by  17  percent.  Correlation  is  less  favorable 
with  higher  loadings  and  the  uniform  inflow  assumption. 

16.  For  the  untwisted  blade  configurations,  correlation  between  test  and 
calculated  vibratory  bending  moment  amplitudes  is  fair  to  poor  with  the 
variable  Inflow  calculation,  and  poor  with  the  uniform  inflow  assump¬ 
tion. 

With  variable  inflow,  untwisted  blade  calculated  flapwise  moment  ampli¬ 
tudes  are  typically  about  80  percent  of  test  values.  Chordwise  calcu- 
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lated  moment  amplitudes  are  typically  50  percent  of  test  values.  Tor¬ 
sional  moments  are  typically  150  to  200  percent  of  test  values. 

With  uniform  inflow  assumed,  the  calculated  untwisted  blade  flapwise 
moment  amplitude  is  typically  50  percent  of  the  test  value.  The  cal¬ 
culated  chordwise  moment  amplitude  is  typically  about  30  percent  of 
the  test  value,  and  the  calculated  torsional  moment  about  120  percent 
of  the  test  value. 

17.  The  percentage  change  in  flapwise  bending  moment  amplitude  due  to  twist 
is  typically  overpredicted  by  the  calculation.  The  calculation  predicts 
an  increase  in  chordwise  bending  amplitude  due  to  blade  twist,  instead 
of  the  decrease  actually  observed  in  the  test  data.  The  agreement  of 
these  bending  moment  ratios  with  test  is  improved  by  using  variable 
inflow  in  the  calculation.  The  percentage  change  in  torsional  moment 
due  to  twist  is  typically  underpredicted  with  the  uniform  inflow  assump¬ 
tion,  and  overpredicted  with  the  variable  inflow  assumption. 

18.  The  percentage  change  in  bending  moment  amplitude  due  to  the  5-degree 
trailing-edge  flap  is  usually  overpredicted  with  the  uniform  inflow 
assumption  and  underpredicted  with  the  variable  inflow  assumption.  The 
percentage  change  in  torsional  moment  amplitude  is  typically  under¬ 
predicted  with  both  inflow  assumptions. 

19.  The  percentage  change  in  bending  moment  amplitude  due  to  the  blade 
stiffness  variation  was  typically  underpredicted  with  the  uniform 
inflow  assumption  and  overpredicted  with  the  variable  inflow  assump¬ 
tion.  The  percentage  change  in  torsional  moment  due  to  the  blade  stiff¬ 
ness  variation  was  typically  overpredicted  with  both  inflow  assumptions. 

20.  Calculated  moment  amplitude  is  moderately  sensitive  to  assumed  wake 
vortex  core  sizes.  Changing  from  an  assumed  radius  of  .3c  to»6c 
lowered  the  flapwise  moment  amplitude  by  30  percent. 
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RECOMMENDATIONS 


1.  Correlation  of  the  model  data  with  full-scale  data  should  be  continued 
for  an  extended  range  of  operating  conditions,  to  further  evaluate  the 
relationship  of  model  and  full-scale  results.  The  correlation  should 
include  peak-to-peak  and  time  histories  of  blade  equivalent  bending 
moment  or  stress,  as  veil  as  rotor  performance  parameters. 

2.  Correlation  of  the  model  test  data  with  calculated  data  should  be  con¬ 
tinued  for  other  operating  conditions  and  analytical  assumptions,  to 
provide  a  more  comprehensive  evaluation  of  rotor  performance  and  blade 
response  calculations. 

3.  In  general,  rotor  performance  and  blade  response  calculations  should 
include  consideration  of  vake-induced  variable  inflow. 

U.  The  analytical  representation  of  rotor  blade  excitation  by  the  vake 
should  be  improved. 

5.  The  analytical  representation,  data  and  associated  scaling  techniques 
for  considering  unsteady  aerodynamics  in  rotor  performance  and  blade 
response  calculations  should  be  augmented  to  include  characteristics 
of  other  symietrical  and  unsymnetrical  airfoils,  other  Reynolds 
number  variations,  and  any  appreciable  effects  of  yawed  flow. 

6.  The  two-dimensional  model  blade  section  aerodynamic  data  presented 
herein  should  be  compared  with  similar  data  collected  in  a  closed 
throat  wind  tunnel. 

7.  Alternate  methods  for  calculating  blade  bending  moment  should  be 
studied  for  possible  incorporation  in  the  normal  modes  aeroelastic 
analysis . 

8.  The  adequacy  of  the  nodal  representation  of  the  rotor  blade  should  be 
studied,  with  reference  to  the  current  limits  on  numbers  of  modes  and 
radial  blade  aegments. 
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TABLK  I.  H-3I4  MODEL  ROTOR  BLADE  INERTIA  AND  CENTER-OE- 
GRAVTTY  DATA 


(lb/ in.)  (lb-sec2 
xlO-  3) 


Bfl 

(lb-sec2 
xlO"  3) 


Inch  YCG 
(lb-sec2  In. 
xlO-  3) 


.23 
.52 
.28 
.20 
.10 
.17 
.11* 
.11* 
.62 
.17 
.07 
.27 
.12 
.22 
.28 
.51* 
.56 
.65 
•  125 
.1*95 
3.630 
.270 

1*0.  m 


3.00 

3.23 

.11*7 

3.23 

3.75 

.11*7 

3.75 

1*.03 

.1*75 

1*  .03 

l*.23 

.11*9 

1*.23 

1* .  33 

.162 

1* .  33 

1* .  50 

.180 

1*.50 

1*  .61* 

.273 

l».6l* 

1*.78 

.31*3 

l*.78 

5. ho 

.31*3 

5.1*0 

5.57 

.31*3 

5.57 

5.61* 

.295 

5.6U 

5.91 

.112 

5.81 

6.03 

.169 

6.03 

6.25 

.169 

6.25 

6.53 

CO 

O 

6.53 

7.07 

.0822 

7.07 

7.63 

.0588 

7.63 

8.28 

.01*11* 

8.28 

8.1*05 

.0176 

8.1*05 

8.900 

.0081 

8.90 

12.530 

.0091* 

12.53 

12.800 

.0106 

12.80 

52.935 

.0121 

52.935 

53.655 

.0202 

53.655 

53.835 

.0178 

53.835 

5h. 155 

.0276 

51*. 155 

55.000 

.001*1* 

.181 

.0036 

.177 

t  * 

.31 

.10 

.33 

.181 

.0036 

.177 

.30 

.331* 

O 

-Cr 

O 

.251* 

.27 

| 

| 

| 

.1(1 

.331* 

.0801 

.251* 

.11* 

.359 

.0928 

.266 

.ll 

.178 

.0892 

.0892 

0.0 

.178 

.0892 

.0892 

.129 

.061*5 

.061*5 

.0351* 

.0177 

.0177 

.0392 

.0193 

.0199 

.0206 

.0111 

.00955 

.00921* 

.001*2 

.00501* 

.00351 

.00037 

.00311* 

.00293 

.00025 

.00268 

.0087 

.00013 

.0086 

.0090 

.00011* 

.ooe9 

.0116 

.oooie 

.0111. 

.0221* 

.00031. 

.0221 

.0221* 

.00031* 

.0221 

.0221* 

.00031* 

.0221 

.0058 

.00009 

.00571 

NOTES:  1.  Mass  moments  of  Inertia  are  with  respect  to  the  local  chord  line  and  in 
axis  normal  to  it  at  the  feathering  axis. 

2.  Mass  moments  of  inertia  include  parts  at  the  root  that  flap,  lag,  and  pitch 
with  the  blade.  Parts  that  flap  and  lag  but  do  not  pitch  are  excluded. 

3.  The  flapwise  position  of  the  blade  center  of  gravity  is  on  the  chord  line 
for  all  blade  stations. 


TABLE  II.  H-3l»  FIBER  GLASS  MODEL  ROTOR  BLADE  STIFFTIESS 
PROPERTIES 


dr 

(In.) 

rln 

(in.) 

rout 

(in.) 

Iflp 
(in.  ) 

Xchd 

(in.4) 

J 

(in 

.-) 

.23 

3.00 

3.23 

.1*66 

2.68 

2. 

37 

•  52 

3.23 

3.75 

.578 

3.51* 

.28 

3.75 

1* .  03 

1.870 

5.32 

.30 

1*  .03 

1*.  33 

.260 

..*.8 

.1*5 

1* .  33 

1».70 

1.600 

i.iOO 

•7 

.62 

1* .  78 

5.1*0 

1.030 

1.0  j0 

1. 

36 

.63 

5.1*0 

6.03 

.011* 

.811* 

1.36 

.22 

6.03 

6.25 

1.120 

1.120 

1.99 

.20 

6.25 

6.53 

.330 

.359 

.5. 

.51* 

6.53 

7.07 

.33c 

.332 

50 

.56 

7.07 

7.63 

.1907 

.11*97 

36. 

.65 

7.63 

8.28 

.0711 

.0712 

0U0 

.125 

8.28 

8.1*05 

.00311*0 

.02190 

0061*1* 

.1*95 

0.1*05 

8.900 

.00196 

.01730 

00393 

3.630 

8.900 

12.530 

.001236 

.0121*0 

0031*0 

.270 

12.53 

12.800 

g 

Css 

r~S 

O 

O 

.0121*0 

00316 

1*1.035 

12.80 

53.035 

.0010*0 

.01152 

00293 

.320 

53.035 

5**.i55 

.00057 

.00625 

. 

00171* 

.81*5 

51*. 155 

55.000 

.00057 

.00625 

.0017'* 

BOTE:  The  following  value*  of  E  and  G  *hould  be  used  with  I„  ,  I  ...  end  J  respectively 
to  obtain  the  correct  bending  and  twisting  stlffnessesror  tne  fiberglass  model 
blade  configurations . 

E  ■  2.50  x  10*  Ib/in? 

0  ■  1.01  x  10s  lb/in? 

The  values  of  J  for  the  sepsents  between  r^  *8.28  and  r  ■  53.835  have  been 
corrected  to  reflect  the  stiffness  test  resdlts  for  the  f?Mrglass  blades. 


TABLE  III.  H-31*  ALUMINUM  MODEL  ROTOR  BLADE  STIFFNESS  PROPERTIES 


Ar 

(in.) 

rin 

(in.) 

r 

out 

(in.) 

IflE 

(in.*) 

*chd 

(in.1*) 

J 

(in. 

» 

.23 

3.00 

3.23 

.117 

.670 

.600 

.52 

3.23 

3.75 

.1**5 

.885 

.28 

3.75 

1* .  03 

.1*68 

1.330 

.30 

*♦ .  03 

**.33 

.0670 

.0670 

.**5 

*♦.33 

**.78 

.1*00 

.1*00 

.600 

.62 

L.78 

5.1*0 

.258 

.258 

.3U1* 

.63 

5.U0 

6.03 

.201* 

.201* 

.31*1* 

.22 

6.03 

6.25 

.280 

.280 

.503 

.28 

6.25 

6.53 

.0825 

.0898 

.132 

.5** 

6.53 

7.07 

.081*0 

.0830 

.127 

.56 

7.07 

7.63 

.01*77 

.0371* 

.0911 

.65 

7.63 

8.28 

.051*56 

.0571* 

.0367 

.125 

8.28 

8.1*05 

.0021*15 

.0177 

.00583 

.*♦95 

8.1*05 

8.90 

.001501* 

.011*0 

.00356 

3.630 

8.90 

12.53 

.00091*8 

.0100 

.00308 

.270 

12.53 

12.800 

.00091*8 

.0100 

.00286 

*♦1.035 

12.80 

53.835 

.000813 

.00929 

.00266 

.320 

53.835 

5**.  100 

.0001*37 

.00501*1 

.0011*5 

.81*5 

5*1.155 

55.000 

.0001*37 

.00501*1 

.0011*5 

NOTE:  The  following  values  of  E  and  S  should 

be  used  with  I 

I  u...  end  j  respectively 

to  obtain  the  correct  bending  and  twist in*  stiffnesses  for  the  alusinun  aodel 

blade  configuration 

E  ■  10.0  *  106 

lb/in.2 

0  ■  *t.00  *  106 

lb/in ? 

The  values  of  I„  ,  end  J  for  the  segments  between  rin  *  8*28  am 

i 

rout  “ 

53.835  have  been 

corrected  to  reflect  the  stiffness  test  results 

for  the 

altoainuB  blade  configurations. 

TABLE  V.  H-3^  !!ODEL  ROTOR  BLADE  MISCELLANEOUS  DATA 


Parameter  or  Description 

Vai  ue 

Uni  ts 

Radius 

55.0 

in . 

Blade  Airfoil  Chord  outboard  of  8.0“  Radial  Station 

2.69 

in. 

Blade  Shank  Cho.-u  inboard  of  8.9"  ^adial  Station 

Blade  Linear  Twist  from  Center  of  Rotation  to  Tip 

1.19 

*  •*  • 

Untwisted  Blades 

0.0 

Twisted  Blade 

Blade  Shear  Center  Locations 

-f . ) 

t  .g 

Beamwise  with  Respect  to  Feathering  Axis 

0.0 

in. 

Chordwise  with  Respect  to  Feathering  Axis 

Location  of  Blade  Hinges 

Radial  Station  of  Coincident  Flapping  and 

0.0 

in. 

Lead-Lag  Hinge 

3.0 

in . 

Feathering  Bearing  Outboard  End  Radial  Station 

5.60 

in. 

Feathering  Bearing  Inboard  End  Radial  Stations 
(The  feathering,  flapping  and  lag  bearing  *xes 
Intersect  at  the  same  point) 

Location  of  Lead-Lag  Damper 

Rotary  Damper  on  Blade  Lead-Lag  Hinge 

L.6L 

in. 

Pitch  Flap  Coupling  Ratio 

0.0 

Damping  Coefficient  of  Lead-Lag  Damper 

(Expressed  as  a  damping  moment  about  tie  hinge) 

Tabs  or  Flaps 

Flapped  Blade  Set  plain  flap  deflection  ever  the  aft 

17.0 

inr lb-sec 

20%  of  the  airfoil  chord 

Pitch  Control  Geometry 

At  0°  collective  and  cyclic  input: 

Distance  aft  along  flapping  hinge  axis  to  pushrod 

5.0 

deg 

upper  end. 

1.1*3 

ir.. 

Radius  on  swashplate  of  pushrod  lower  end 

Distance  of  plane  of  swashplate  below  plane  of 

1.95 

in. 

flapping  hinges 

Angle  on  swashplate  between  pushrod  lover  end 

6.31 

in. 

and  plane  of  shaft  and  lead-lag  axis 

1*0.5 

deg 

Pitch  Control  Spring  Rate 

Airfoil  Section  -  IACA  0012 

21,000. 

ir..-lb/rad 

TABLE  VI.  MODEL  BLADE  SECTION  TWO-DIMENSIONAL  AERODYNAMIC  DATA 
(HACA  0012  Airfoil,  6-  *  0.0  deg,  c  *  2.69  in., 

M  «  .10,  RN  «  .155  X  iO^) 


(deg) 

1 

(dec< 

4 

(dec) 

l 

• 

0.0 

.00 

0.0 

.0000 

0.0 

0.0 

5.4 

.45 

1.0 

.0084 

4.0 

.003 

7.0 

.59 

2.5 

.0103 

8.0 

.008 

8.0 

.675 

4.0 

.0150 

8.6 

.007 

8.3 

.71 

6.5 

.0295 

9.0 

.004 

9.4 

.74 

8.5 

.054 

11.0 

-.035 

10.0 

.73 

10.0 

.100 

12.0 

-.051 

12.0 

.71 

11.5 

.140 

13.0 

-.063 

14.0 

.68 

12.5 

.170 

14.0 

-.071 

16.0 

.67 

14.0 

.200 

16.0 

-.079 

13.0 

.67 

29.0 

.530 

22.0 

-.097 

20.0 

.69 

45.0 

1.150 

24.0 

-.106 

23.0 

.73 

55.0 

1.420 

26.0 

-.120 

26.0 

.  3i 

63.0 

1.600 

32.0 

-.177 

33.0 

1.00 

70.0 

1.730 

41.0 

-.250 

39.0 

1.07 

75.0 

1.800 

55.0 

-.325 

45.0 

1.06 

80.0 

1.850 

72.0 

-.395 

50.0 

1.03 

85.0 

1.680 

90.0 

-.460 

60.0 

.65 

89.0 

1.890 

107.0 

-.500 

75.0 

.50 

93.0 

1.880 

117.0 

-.513 

92.0 

0.00 

100.0 

1.850 

125.0 

-.510 

120.0 

-.72 

105.0 

1.8o 

133.0 

-.495 

i30.0 

-.69 

112.0 

1.71 

140.0 

-.465 

137.0 

-.94 

120.0 

1.55 

147.0 

-.420 

143.0 

-.94 

130.0 

1.28 

155.0 

-.360 

148.0 

-.90 

145.0 

0.75 

158.0 

-.290 

*56.0 

-.70 

150.0 

.52 

162.0 

-.270 

159.0 

-.67 

155.0 

.40 

165.0 

-.275 

163.0 

-.66 

160.0 

.32 

171.0 

-.310 

167.0 

-.69 

164.3 

.28 

173.0 

-.290 

171.0 

-.75 

170.0 

.14 

178.0 

-.150 

1'5.3 

-.60 

174.0 

.06 

180.0 

0.0 

100.0 

0.00 

180.0 

.02 

TABLE  VII 


MODEL  BLADE  SECTION  TWO-DIMENSIONAL  AERODYNAMIC 
DATA  (IJACA  0012  Airfoil,  6  =  0.0  deg,  c  =  2.C9  in., 

M  =  .18,  Kl  =  .279  X  10°)  h 


t —  ■  ■  - - 

C1 

(dec) 

Cd 

a 

( deg ) 

c 

a 

0.0 

0.0 

0.0 

.0080 

0.0 

0.0 

8.1* 

.72 

0.5 

.008 2 

8.0 

.012 

8.8 

.74 

1.5 

.0088 

9.6 

.012 

9.3 

.76 

2.5 

.0100 

9.2 

.008 

9.8 

.76 

5.0 

.0175 

9.5 

0.0 

n.o 

.72 

7.5 

.036 

9.8 

-.014 

12.9 

.68 

10.2 

.10 

10.4 

-.028 

15.3 

.65 

10.7 

.115 

11.0 

-.040 

17.1 

.64 

11.8 

.145 

12.0 

-.054 

19.2 

.645 

13.0 

.170 

13.7 

-.067 

21.4 

.67 

15.0 

.210 

l6.o 

-.074 

25.0 

.745 

19.0 

.290 

19.6 

-.083 

33.0 

.97 

30.0 

.60 

22.0 

-.092 

37.0 

.99 

41.0 

1.00 

24.0 

-.110 

41.0 

.99 

49.0 

1.25 

27.6 

-.132 

46.5 

.965 

56.0 

1.40 

32.0 

-.176 

52.0 

.92 

61.5 

1.49 

45.0 

-.240 

64.0 

.70 

69.0 

1.57 

58.0 

-.314 

79.5 

.38 

73.0 

1.615 

70.0 

-.365 

115.0 

-.60 

80.0 

1.66 

90.0 

-.433 

122.0 

-.75 

89.0 

1.69 

113.0 

-.49 

128.5 

-.845 

96.0 

1.685 

119.0 

-.485 

134.0 

-.88 

105.0 

1.65 

125.0 

-.482 

140.0 

-.89 

115.0 

1.55 

132.0 

-.470 

144.5 

-.85 

125.0 

1.36 

139.0 

-.450 

148.0 

-.79 

132.0 

1.19 

145.0 

-.410 

153.0 

-.66 

137.5 

1.02 

150.0 

-.365 

155.0 

-.62 

151.5 

.49 

156.0 

-.295 

159.0 

-.60 

155.0 

.40 

161.0 

-.270 

162.5 

-.62 

160.5 

.30 

163.5 

-.266 

157.5 

-.67 

165.5 

.23 

167.0 

-.275 

172.0 

-.65 

171.0 

.11 

170.0 

-.295 

175.0 

-.59 

174.5 

•  05 

175.0 

-.240 

180.0 

0.0 

180.0 

.015 

l80.0 

0.0 

TABLE  VIII.  MODEL  BLADE  SECTION  TWO-DIMENSIONAL  AERODYNAMIC 
DATA  (NACA  0012  Airfoil,  5p  =  0.0  deg, 
c  =  2.69  in.,  M  =  .33,  RN  =  .512  X  1()6) 


d  Cl  01 

(deg)  C1  (deg)  'd  (deg) 


0.0 

0.0 

C.O 

.0081 

0.0 

0.0 

7.9 

.715 

.50 

.0082 

8.1* 

.015 

8.5 

.75 

1.5 

.0088 

8.9 

.Oil* 

8.9 

.76 

2.5 

.010 

9.3 

.010 

9.1* 

.75 

3.5 

.012 

9.7 

-.030 

10.8 

.63 

7.5 

.032 

10.1 

-.01*5 

11.1* 

.60 

8.1* 

.01*1* 

10.7 

-.053 

12.0 

.575 

9.0 

.060 

11.1* 

-.058 

13.1 

.56 

9.6 

.10 

18.0 

-.076 

lb.O 

.56 

9.8 

.12 

22.0 

-.092 

15.5 

.565 

10.5 

.11* 

2*».0 

-.102 

lS.O 

.59 

11.5 

.158 

26.0 

-.115 

21.0 

.635 

11*. 0 

.195 

28.8 

-.138 

2l*.7 

.710 

20.0 

.30 

30.0 

-.11*9 

28.6 

.80 

30.0 

.60 

32.0 

.88 
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TABLE 

IX.  MODEL  BLADE  SECTION 
(NACA  0012  Airfoil, 
M  *  .1*7,  RN  «  .729  > 

TWO-DIMENSIONAL  AERODYNAMIC  DATA 

6p  »  0.0  deg,  c  =  2.69  in., 

106) 

a 

(deg) 

C1 

a 

(deg) 

Cd 

a 

(deg) 

c 

B 

0.0 

0.0 

0.0 

.0080 

O.C 

0.0 

5.7 

.59 

1.0 

.0081 

1.6 

.002 

6.0 

.615 

2.0 

.0087 

3.1* 

.006 

6.5 

.63 

3.0 

.010 

5.2 

.012  ( 

7.0 

.635 

5.5 

.018 

6.1* 

.018 

7.6 

.62 

6.25 

.0235 

6.8 

.018 

9.0 

.57 

6.8 

.030 

7.2 

.016 

11.0 

.535 

7.2 

.01*0 

8.1* 

-.031* 

12.6 

.53 

7.7 

.070 

9.0 

-.01*7 

1U.0 

.53 

8.0 

.080 

9.8 

-.055 

15.6 

.51* 

8.6 

.10 

10.8 

-.058 

17.3 

.56 

9.5 

.125 

13.0 

-.063 

19.0 

.567 

11.0 

.151* 

li*.6 

-.070 

21.0 

.557 

12.5 

.182 

16.3 

-.076 

22.6 

.525 

IV. 5 

.2 18 

18.0 

-.085 

21*.  0 

.1*75 

16.5 

.25 

20.0 

-.092 

20.25 

.30 

22.0 

-.096 

22.5 

.31 

2l*.0 

-.10 

25.0 

.285 

28.0 

-.102 
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TABLE  X.  MODEL  BLADE  SECTION  TWO-DIMENSIONAL  AERODYNAMIC  DATA 
(NACA  0012  Airfoil,  6_  *  0.0  deg,  c  *  2.69  in., 

M  *  .56,  RN  *  .368  X  lO6) 


(deg) 

C1 

(deg) 

Cd 

a 

(deg) 

c 

m 

0.0 

0.0 

0.0 

.0080 

0.0 

0.0 

3.6 

.1*0 

.50 

.0081 

1.1* 

.0010 

U.O 

.1*3 

1.5 

.0087 

2.7 

.001*0 

l*.9 

•  **55 

2.0 

.0091* 

t.5 

.0090 

5-7 

.1*6 

2.75 

.011 

5.2 

.010 

6.2 

.1*55 

3.5 

.015 

5.8 

.010 

7.** 

.1*1* 

5.5 

.01* 

6.5 

.007 

8.5 

.1*35 

7.75 

.10 

7.0 

.002 

10.7 

.1*1*5 

8.25 

.112 

7.1* 

-.001* 

13.1 

.1*5 

9.0 

CO 

C\J 

8.0 

-.020 

H..6 

.1*1*5 

10.6 

.11*8 

8.8 

O 

o 

• 

16.2 

.1*35 

12.0 

.16 

9.5 

o 

1A 

o 

• 

17.6 

.1*2 

13.25 

.175 

10.3 

VO 

O 

1 

19.0 

.395 

15.0 

.187 

11.3 

-.060 

19.6 

co 

17.5 

.20 

13.1 

-.062 

20.0 

.20 

16.0 

$ 

O 

1 

21*.  0 

-.082 

TABLE  XI. 

MODEL 
AC  A 
K  =  . 

BLADE  SECTION 
0012  .drfoil, 
10,  Ki  =  .155  X 

TWO-DIMENSIONAL 
=  5.0  deg,  c 
IQ6) 

AERODYNAMIC 
=  2.6 9  In., 

DATA 

a 

a 

ci 

(deg) 

Cd 

(deg) 

c 

a 

-180.0 

.19 

-180.0 

.025 

-180.0 

.065 

-176.0 

.59 

-176.0 

.050 

-175.0 

.275 

-170.0 

7' 

-168.0 

.21 

-162.5 

.236 

-164.5 

.59 

-158.0 

.34 

-146.5 

.392 

-143.0 

.91 

-127.0 

1.375 

-137.0 

.44 

-136.0 

.92 

-109.5 

1.72 

-125.0 

.47 

-123.0 

.75 

-  93.5 

1.855 

-110.0 

,46^ 

-  62.0 

-.77 

-  84.0 

1.S35 

-  85.0 

.4? 

-  47.0 

-.985 

-  74.0 

1.7  2 

-  60.0 

.31 

-  39.0 

-.995 

-  59.0 

1.43 

-  30.0 

.128 

-  32.0 

-.93 

-  30.0 

.  56 

-  26.0 

.095 

-  24.0 

-.695 

-  13.6 

.17 

-  20.0 

.066 

-  17.5 

-  615 

-  12.0 

.145 

-  12.0 

.040 

-  13. 4 

-.54 

-  11.0 

.11 

-  10.3 

-.040 

-  11.6 

-.585 

-  10.5 

.050 

-  8.4 

-.031 

-  10.3 

-.69 

-  5.3 

.014 

-  5.3 

-.030 

-  8.0 

-.575 

-  2.25 

.0066 

-  1.4 

-.u*.2 

0.0 

.20 

-  2.0 

.0081 

0.0 

-.040 

5.0 

.64 

0.0 

.0090 

3.0 

-.020 

7.8 

.80 

8.25 

.058 

8.0 

-.011 

8.9 

.79 

9.0 

.120 

10.0 

-.081 

11.1 

.63 

10.0 

.153 

18.0 

-.103 

14.2 

.66 

14.0 

.220 

24.0 

-.135 

17.6 

.69 

30.0 

.70 

30.0 

-.185 

32.0 

1.03 

52.0 

1.40 

60.0 

-.360 

36.5 

1.08 

68.0 

1.715 

85.0 

-.450 

45.0 

1.07 

82.5 

1.85 

110.0 

-.495 

61.0 

.81 

90.0 

1.87 

120.0 

-.499 

116.0 

-.67 

104.0 

1.77 

135.0 

-.468 

132.0 

-.94 

120.0 

1.49 

146.5 

-.40 

142.0 

-.94 

157.5 

.34 

154.0 

-.30 

158.5 

-.64 

163.5 

.245 

158.5 

-.285 

168.5 

-.78 

175.5 

.035 

171.5 

-.315 

180.0 

.19 

180.0 

.025 

180.0 

.065 

TABLE  XII. 

MODEL  BLADE  SECTION 
OIACA  0012  Airfoil, 
M  =  .18,  XI  =  .<r79 

T’.VO-D  I  MENS  I ONAL  AERODYU  AMI  C 
=  5.0  deg,  c  *  2.69  in., 

(  10tJi 

DATA 

_ 

Cl 

a 

a 

(deg) 

ci 

(deg) 

Cd 

( deg) 

C 

m 

-180.0 

.25 

-180.0 

.025 

-180.0 

.063 

-177.0 

.58 

-175.0 

.070 

-175.0 

.275 

-172.5 

.63 

-169.5 

.175 

-165.0 

.21* 

-16I*.0 

•  57 

-158.0 

.35 

-159.5 

.255 

-11*1.0 

.85 

-138.5 

1.00 

-11*5.0 

.375 

-128.0 

.71* 

-111.5 

1.1*7 

-130.0 

.1*2 

-  69.0 

-.51 

-  97.0 

1.62 

-110.0 

.1*25 

-  1*9.5 

-.83 

-  87.5 

1.65 

-  90.0 

.38 

-  1*0.0 

-.88 

-  82.0 

1.635 

-  60.0 

.29 

-  32.0 

-.81 

-  70.0 

1.535 

-  30.0 

.118 

-  23.0 

-.605 

-  55-0 

1.30 

-  22.0 

.061* 

-  16.9 

-.51* 

-  30.0 

.52 

O 

CD 

r-4 

I 

.01*8 

-  12.9 

-.505 

-  13.3 

.150 

-  11*.  5 

.01*7 

-  11.3 

-.675 

-  11.8 

.113 

-  13.0 

.01* 

-  9.3 

-.59 

-  11.3 

.01.7 

-  12.2 

.027 

-  7.6 

-.1*8 

-  1**5 

.0106 

-  11.6 

-.036 

0.0 

.23 

-  2.8 

.0085 

-  10.8 

-.01*0 

5.1* 

.68 

-  1.5 

.0082 

-  I*.0 

-.01*5 

7.9 

.92 

0.0 

.0090 

0.0 

-.038 

8.6 

.79 

1.5 

.0121 

7.0 

-.016 

10.8 

.62 

8.3 

.050 

8.0 

-.016 

13.6 

.61* 

9.1* 

.130 

9.6 

-.073 

17.0 

.61* 

10.5 

.16 

11.5 

-.007 

22.9 

.735 

30.0 

.65 

17.0 

-.093 

27.3 

.89 

1*9.0 

1.15 

23.0 

-.12 

33.5 

.96 

66.0 

1.1*5 

30.0 

-.168 

38.5 

.965 

81*.  0 

1.595 

60.0 

-.313 

1*1*.  0 

.9*. 

100.0 

1.635 

90.0 

-.1*0 

69.0 

.61 

110.5 

1.56 

110.0 

-.1*23 

120.0 

-.71*5 

137.0 

1.00 

130.0 

-.1*1 

11*1.0 

-.85 

155.5 

.35 

11*5.0 

-.36 

161.5 

-.615 

165.5 

.20 

159.0 

-.21*5 

173.0 

-.71*5 

173.0 

.065 

171.0 

-.315 

180.0 

.25 

180.0 

.025 

180.0 

.063 

TABLE  XIII.  MODEL  BLADE  SECTION 
(NACA  0012  Airfoil, 
M  =  .33,  RN  =  .512 

Ta'O-DI MENS  I OL’AL  AER  DYNAMIC  DATA 

*  5.0  deg,  c  =  2.6 9  in., 

X  106) 

— 

C1 

a 

(deg) 

Cd 

a 

(deg) 

c 

m 

-30.0 

-.725 

-30.0 

.50 

-30.0 

.106 

-25.0 

-.605 

-17.5 

.22 

-22.0 

.066 

-22.0 

-.55 

-13.1* 

.16 

-18.0 

.05 

-l8.6 

-.53 

-12.3 

.11* 

-16.0 

.01*5  | 

-15.0 

-.1*85 

-11.6 

.125 

-11*. 6 

.01*3 

-13.3 

-.U9 

-11.3 

.10 

-12.1* 

.035  j 

-11.7 

-.51*5 

-11.0 

.06 

-11.8 

.022 

-10.  J* 

-.65 

-10.8 

.05 

-11.1 

-.022 

-  9.2 

-.59 

-  9.5 

.035 

-10.7 

-.030 

-  8.2 

-.53 

-  1*.5 

.010 

-  9.6 

-.035 

0.0 

.22 

-  3.5 

.0088 

-  8.0 

-.037 

5.2 

.65 

-  2.0 

.0080 

-  6.0 

-.036 

6.6 

.71*5 

-  1.0 

.0082 

0.0 

-.028 

8.1 

.81 

0.0 

.0088 

2.6 

-.025 

8.9 

.70 

6.3 

.035 

5.0 

CO 

H 

O 

1 

10.0 

.62 

7.5 

.01*5 

7.0 

-.007 

11.6 

.60 

8.3 

.055 

8.0 

-.001. 

13.1 

.62 

0.5 

.069 

8.6 

1 

0 
0 
— J 

15.3 

.605 

8.8 

.110 

9.0 

-.026 

17.1 

.625 

9.3 

.135 

9.2 

-.056 

20.8 

.68 

10.0 

.156 

10.0 

-.071* 

28.0 

.815 

12.0 

.19 

11.0 

-.080 

30.0 

.85 

16.8 

.27 

12.0 

-.082 

30.0 

.61* 

15.0 

-.087 

19.0 

-.103 

27.0 

-.11*2 

30.0 

-.15 

TABLE 

XIV.  MODEL 
(NACA 
M  =  . 

BLADE  SECTION 
0012  Airfoil, 
1*6,  RN  *  .713 

TWO-DIMENSIONAL  AERODYNAMIC  DATA 
*  5.0  deg,  c  =  2.69  in., 
x  io6) 

a 

(deg) 

C1 

(deg) 

Cd 

a 

(de<) 

C» 

-30.0 

-.71 

-30.0 

.60 

-30.0 

.126  ) 

-21.3 

-53 

-23.5 

.36 

-22.0 

.071* 

-17.9 

-.1*9 

-18.5 

.23 

-18.0 

.053 

-1U.8 

-.1*8 

-15.0 

.175 

-15.0 

.01*5 

-11.6 

-.51 

-12.0 

.11*2 

-13.0 

.01*3 

-10.1* 

-.51* 

-11.3 

.130 

-11.6 

.035 

-  9.U 

-.615 

-10.5 

.110 

-10.6 

.016 

-  7.9 

-.50 

-10.3 

.100 

-  9.9 

-.020 

0.0 

.25 

-  9.0 

.050 

-  9.3 

-.033 

3.1 

.51*5 

-  7.9 

.030 

-  8.6 

-.036 

3.5 

.55 

-  5.0 

.0115 

-  5.0 

-.03; 

l*.7 

.51 

-  i*.o 

.0091* 

0.0 

-.027 

6.1 

.515 

-  2.8 

.0082 

5.0 

-.02* 

7.0 

.55 

-  1.5 

.0080 

6.5 

-.020 

8.1 

.58 

-  .60 

.0084 

7.0 

-.050 

9.7 

.59 

0.0 

.0090 

7.8 

-.066 

13.7 

.565 

5.0 

.060 

9.6 

-.077 

18.0 

.55 

6.0 

.092 

15.0 

-.089 

23.0 

.1*8 

6.8 

.110 

21*. 0 

-.121 

30.0 

.38 

8.0 

.135 

30.0 

-.11*3 

10.0 

.165 

13.5 

.210 

17.5 

.26 

30.0 

.1*5 

I  J 


\ 
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TABLE  XVI.  CONDITION  DESCRIPTIONS 


Condition 

Numbers 

1,  1>*,  35 

2,  3,  T,  15,  19-21 

U-6.  16-18,  36,  37,  50-53 

8-13,  22-31,  38-1*5 

32-3U ,  U6-l*9 

5**-6l 

62-67 

68-71,  76-90 
72-75 
91,  93 
92 


Description 

Parameter  Variation  Base  Point 
Collective  Variation 
Flapping  Variation 

Shaft  Angle  and  Collective  Variation 
Blade  and  Wake  Intersection 
Blade  Resonance* 

High  u 
Blade  Stall 

Blade  Stall  with  Flapping  Variation 
Hover,  Zero  Lift,  B  R  «  700  FT  /SEC 
Hover,  Zero  Lift,  flgR  ■  lOO  FT /SEC 


•Run  for  3  *  scale  stiffness  alimlnum  blades  only 


1  ABLE  XVII.  BLADE  CONFIGURATIONS,  TRANSDUCERS ,  AND  CONDITIONS  WITH  VALID  BLADE  RESPONSE 
DATA  (Blade  Configuration:  Q  =  0  deg,  =  0  deg,  S„  =  1) 
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TABLE  XXII.  H-3U  MODEL  BLADE  PERFORMANCE  DATA  (Configuration:  0  =  0  deg 

5„  =  0  deg,  S  *  l) 
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TABLE  XXIII.  H-3U  MODEL  BLADE  PERFORMANCE  DATA  (Configuration.  0  =  -8  deg 

6p  »  0  deg,  Sp  *  1) 
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TABLE  XXIV.  H-31*  MODEL  BLADE  PERFORMANCE  DATA  (Configuration:  0.  =  0  deg 
6y  *  5  deg,  SP  *  1) 
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TABLE  XXV.  H-3l*  MODEL  BLADE  PERFORMANCE  DATA  (Configuration:  0=0  deg 
6p  *  0  deg,  Sp  =  3) 
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TABLE  XXXIV.  TEST  AND  CALCULATED  ROTOR  PERFORMANCE  (Condition  1, 
01  *  0  deg,  6p  =  0  deg,  Sg  =  l) 
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TABLE  XL.  SUMMATION  OF  MODAL  CONTRIBUTIONS  TO  FLAPWISE  MOMENT 
RESPONSE  AT  .20R  (Condition  7t  ip  =30  deg.  Uniform 
Inflow,  G1  =  0  deg,  6F  =  0  deg,  S£  =  l) 


Bending 

Modal 

Modal 

Curvature 

Summations 

Mode 

Amplitude 

Curvature 

Contribution 

Curvature  Moment 

Number 

(in. -lb) 

i 

Y  " 

.'wi 
( .20R) 

Wv"i 

.009171 

10.63 

1  -.0011*59 

-51**82 

• 

1  .0001520 

91.02 

. 

-.0003857 

-111*.  3 

• 

.097U8 


•  xi  f 

.1913 

.2351* 
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TABLE  XLI . 


SUMMATION  OF  MODAL  CONTRIBUTIONS  TO  CHORDWISE  MOMENT 
RESPONSE  AT  .20R  (Condition  7,  p  =  50  deg.  Uniform 
Inflow,  0^  =  0  deg,  6^  =  0  deg,  Sj.  =  1 ) 


Bending 

Mode 

Number 

Modal 

Amplitude 

Modal 

Curvature 

Curvature 

Contribution 

Summations 
Curvature  Moment 

(in. -lb) 

J 

*vj 

*vJ 

(.20R) 

qvJYvJ 

SlTd  y" 

Tyvjvj 

I 


-.0011*12 

’  .0003062 

-1* 

-.01111* 

"  • 

-.011*05 

”  a 

TABLE 

XLII .  SUMMATION  OF  MODAL  CONTRIBUTIONS 

RESPONSE  AT  .65R  (Condition  7»  ^ 

Uniform  Inflow,  0=0  deg,  <$  = 

X  F 

TO  FLAP WISE 
=  120  deg, 

0  deg,  S£  = 

MOMENT 

1) 

Bending 

Modal 

Modal 

Curvature 

Summations 

Mode 

Humber 

Amplitude 

Curvature 
( .65R) 

Contribution 

Curvature 

Moment 
( in. -lb) 

i 

qvi 

»» 

Yvi 

a  .y" 

-vi  vi 

Z'lwi^vi 

Elya  / . 

VI 

.UO1970 

11. 2h 

.02223 

.02223 

1.07 

-.001073 

30.51 

-.03271* 

-.0105 

-.503 

.000191*1 

10.81 

.002101* 

-.0081* 

-.1*03 

mm 

-.00001*31*8 

-110.81 

.001*818 

-.00358 

-.172 

TABLE 

XLIII.  SUMMATION  OF  MODAL  CONTRIBUTIONS 

RESPONSE  AT  .65R  (Condition  7,  * 

Uniform  Inflow,  0*0  deg,  6  = 

X  r 

TO  FLAPWISE  MOMENT 
=  280  deg, 

0  deg,  Sg  =  l) 

Bending 

Modal 

Modal 

Curvature 

Summations 

Mode 

Amplitude 

Curvature 

Contribution 

Curvature 

Moment 

Number 

(.65R) 

(in. -lb) 

1 

“vi 

If 

Yvi 

Wvi 

SVCt 

El  Vn  Y" 
R^vivi 

1 

.009089 

11. 2l* 

.1022 

.1022 

U.90 

2 

-.0003008 

30.51 

-.009177 

.0930 

U.U6 

3 

-.00008789 

10. 8U 

-.0009527 

.0920 

k.kl 

1* 

-.0001682 

-110.81 

+.0186U 

.1106 

5.31 
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TABLE 

XLIV.  SUMMATION  OF  MODAL 

CONTRIBUTIONS 

TO  FLAPWISE 

MOMENT 

RESPONSE  AT  .65R  (Condition  7»  <P 

=  120  deg. 

Uniform 

Inflow,  0^ 

=  -8  deg,  6p  = 

0  deg,  3^ 

-  1) 

Bending 

Modal 

Modal 

Curvature 

Summations 

Mode 

Amplitude 

Curvature 

Contribution 

Curvature 

Moment 

Number 

(.65R) 

(in. -lb) 

i 

Sri 

Y". 

wi 

Srt^i 

tJVU 

1 

-.01059 

11.21* 

-.1190 

-.1190 

-5.71 

2 

-.001702 

30.51 

-.0519 

-.1709 

-8.20 

3 

.OOOU925 

10.81* 

.0053 

-.1656 

-7. 91* 

U 

-.00011*37 

-110.81 

.oi;9 

-.11*97 

-7.18 

TABLE  XLV.  SUMMATION  OF  MODAL  CONTRIBUTIONS  TO  FLAPWISE  MOMENT 
RESPONSE  AT  .65R  (Condition  7,  <(,  *  280  deg. 

Uniform  Inflow,  *  -8  deg,  6p  ■  0  deg,  Sg  *  l) 


Bending 

Mode 

Number 

Modal 

Amplitude 

Modal 

Curvature 

(.65R) 

Curvature 

Contribution 

Summation 

Curvature  Moment 

( in. -lb) 

i 

qwi 

Y  " 

'  vi 

VYvi 

ZqwiYwi 

1 

.009286 

11. 2l* 

.101*1* 

.101*1* 

5.011 

2 

-.000321 

30.51 

-.0098 

.091*5 

1* .  51* 

3 

.0000001*1*21 

10.81* 

.0000 

.091*5 

1* .  51* 

1* 

-.0001095 

-110.81 

.0121 

.1067 

5.12 
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Figure  1.  Helicopter  Model  With  H- Model  Blades. 
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Figure  5.  Blade  Deflection  Test  and  Calibration  Stand 
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Figure  6.  H-31*  Model  Rotor  Blade  Flapvise  and  Chordvise  Deflections,  Blade  Clasped 

at  6.25-Inch  Radius,  Three-Times-Scale  Stiffness  Aluminum  Spar. 
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Figure  8.  H-3^  Model  Rotor  Blade  Flapvise  and  Chordvise  Deflections,  Blade  Clamped 

at  6.25-Inch  Radius,  Scale  Stiffhess  Fiber  Glass  Spar. 
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Figure  9.  H-31*  Model  Rotor  Blade  Torsional  Deflections,  Blade  Clanped  at  6.25-Inch 

Radius,  Scale  Stiffness  Fiber  Glass  Spar. 


r’igure  10.  Blade  Natural  Frequency  Test  Setup. 
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Figure  11.  Natural  Frequencies  Versus  Rotational  Tip  Speed  for 
Three -Tines -Scale  Stiffness  Aluminum  Spar  H-3k  Model 
Blades . 


Figure  12.  Natural  Frequencies  Versus  Rotational  Tip  Speed  for 
Dynamically  Scaled  Fiber  Glass  H-31*  Model  Blades. 
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Figure  13.  Two-Dimensional  blade  Section  V.'ind  Tunnel  Test  Setu 


Airfoil  Section. 
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Figure  lU.  Concluded. 


COEFFICIENT  OF  LIFT 


-30  -20  -10  0  10  20  30 

ANGLE  OF  ATTACK, o  ,  DEG 


(a)  Coefficient  of  Lift  -30°  <  a  <  30° 


Figure  15.  Two-Dimensional  Aerodynamic  Coefficient*  for  the  2.69-Inch 
Chord  Flapped  Airfoil  Section. 
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Figure  15.  Continued. 


COEFFICIENT  OF  DRAG. 


ANGLE  OF  ATTACK,  a  ,  OEG 


(c)  Coefficient  of  Dreg  -30°  <  a  <  30° 


Figure  15 .  Continued. 
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Figure  15.  Continued 


COEFFICIENT  OF  PITCHING  MOMENT 


Figure  16.  Comparison  of  Lov  Reynolds  Number  Two-Dimensional  Airfoil  Test  Data. 


5  x  10" 


(b)  Lateral  Response 


14x10**- 


Q I  I  I  I  I  I  1  1  1  1  1  1  ■  A .  1  J 

24  28  32  36  40  44  48  52 

FORCING  FREQUENCY  ,  uf  ,  CPS 


(d)  Torsional  Response 


Figure  17.  Concluded. 
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PHASE  ANGLE  .  *UA  .  DEG 


LATERAL  INERTIA  ,  **tfry  .  LB  SEC*/ IN 


(a)  Effective  In-Plane  Inertia 

Figure  18.  Calculated  Effective  Ebmanlc  Propertle*  at  tbe  Hot  or  Hub. 
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FORE-AFT  INERTIA  .  MCfrx  .  LB-SEC*/IN 


(b)  Effective  In-Plane  Spring  Rate 

Figure  18.  Continued. 
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FORE-AFT  SPRING  CONSTANT  ,  KEFFX  ,  LB/ IN. 


LATERAL  DAMPING  COEFFICIENT  ,  CEFFY  .  LB -SEC/ IN. 


i - 1 - 1 - 1 - 1 - 1  | 

25  30  35  40  45  50 

FORCING  FREQUENCY  ,  wF  ,  CPS 


(c)  Effective  In-Plane  Damping 


Figure  18.  Continued. 
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FORE-AFT  DAMPING  COEFFICIENT  ,  CEFFX  ,  LB-SEC /IN. 


FORCING  FREQUENCY  ,  wF  ,  CPS 

(d)  Effective  Vertical  and  Torsional  Inertia 

Figure  18.  Continued. 
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TORSIONAL  INERTIA  ,  MEFFT  ,  IN. -LB-SEC2 


FORCING  FREQUENCY  ,  uf  ,  CPS 


(e)  Effective  Vertical  and  Torsional  Spring  Rate 


Figure  18.  Continued. 


VERTICAL  DAMPING  COEFFICIENT  ,  CEFFV  ,  LB  SEC /IN. 


(f)  Effective  Vertical  and  Torsional  Damping 
Figure  18.  Concluded. 
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TORSIONAL  DAMPING  COEFFICIENT  ,  C£FFT  ,  IN.-LB-SEC 
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TYPICAL  CR 


Figure  ?0.  H-3^  Model  Blade  Instrumentation. 
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Figure  21.  Concluded 


ROTOR  DRAG  COEFFICIENT  -  SOLIDITY  RATIO 
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■004 
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-.008 
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— -  -  -[■  — — — - 1 - J - j - 1 — \ 

— o-  8,  *0* .  8f  »0*  ,  FIBER  GLASS  SPAR  MODEL  N 
- O  8,  *0*  .  8f«  5*  ,  FIBER  GLASS  SPAR  MODEL 

— A-  a,  «-8#,8f*0#,  fiber  glass  spar  model 
— 0-  a,  *0*  ,  8f  * 0*  ,  aluminum  spar  , 

3X  SCALE  STIFFNESS  MODEL 

-  -O"  a,  «0*  .  FULL  SCALE  H-34  REF  ( 3  )  - 

-  -Ar-  a,  «-8#  ,  FULL  SCALE  H-34  REF(  3  ) 


ROTOR  LIFT  COEFFICIENT  -  SOLIDITY  RATIO,  CL/<r 


Figure  22.  Drag  vs.  Lift  Performance  for  Model  Rotor  Configurations; 
Vg  *  kn,  u  8  0.3,  a^_  8  b  =  0  deg,  ag  8  0  deg. 
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ROTOR  DRAG  COEFFICIENT  -  SOLIDITY  RATIO 


ROTOR  DRAG  COEFFICIENT  -  SOLIDITY  RATIO  , 


ROTOR  LIFT  COEFFICIENT  •  SOLIDITY  RATIO,  C  L/tr 


(b)  u  ■  0.1»,  V,  »  166  kn 


Figure  2U.  Continued. 
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ROTOR  DRAG  COEFFICIENT  - 


Cl/* 


(c)  w  •  0.5.  V§  ■  208  kn 


Figure  2*».  Concluded. 
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ROTOR  LIFT  COEFFICIENT  -  SOLIDITY  RATIO,  C^/r 


(a)  v  •  0.3,  Vg  »  12U  kn 


Figure  2 5.  Torque  vs.  Lift  Performance  of  Model  Rotor  at  Vari 
Advance  Ratios;  6=0  deg,  5_  *  0  deg,  S  *  1,  a. 


ROTOR  TOROUE  COEFFICIENT  - 
SOLIDITY  RATIO  ,  C Q/<r 


ROTOR  TORQUE  COEFFICIENT 
SOLIDITY  RATIO  ,  Cq/o- 


.016 


.012 


0  .02  .04  .06  .06 

ROTOR  LIFT  COEFFICIENT  -  SOLIDITY  RATIO  , 

C\_/v 


(c)  u  ■  0.5,  V,  ■  208  kc 


Figure  25.  Concluded. 
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Figure  27.  Cyclic  Pitch  to  Trim  Flapping  vs.  Collective  Pitch  at  Various  Advance  Ratios; 
a.  ■  0  deg,  5_  *  0  deg,  S_  *  1,  a.  *  h.  *  0  deg. 


930*  S,V  '39NIH  IV  ONIddVIi 
TVd31Vl  Midi  01  INSNOdWOO  HDild  911019 


J3 


w 

H 

■ 


176 


Figure  27.  Continued. 
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Figure  28.  Continued. 


COLLECTIVE  PITCH  AT  .75R  ,  Bc  ,  DEG 


29*  Pi  Mentation  and  Correlation  of  Model  Rotor  Performance  vs.  Collective  Pitch; 
Vg  -  121*  kn,  9  -  -8  deg,  6  «  0  deg,  S  «  1,  a  =  0  deg,  V  =  0.3. 


(b)  Rotor  Drag  Coefficient-Solidity  Ratio 


(c)  Rotor  Torque  Coefficient  Solidity  Ratio 
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Figure  30.  Continued 


i  aoioa 


COLLECTIVE  PITCH  AT  .75 R 


(d)  Cyclic  Pitch  to  Trxai  Flapping 


Figure  31.  Presentation  and  Correlation  of  Model  Rotor  Performance  vs.  Collective  Pitch* 
V8  «  12h  kn,  flx  -  0  deg,  5p  -  5  deg,  S  -  1,  o  «  0  deg,  v  »  0.3. 
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0  NORMAL  MOOES  AE  ROE  LAST  1C  ANALYSIS 
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(c)  Rotor  Torque  Coefficient-Solidity  Ratio 


Figure  31.  Continued. 
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(u)  Cyclic  Pitch  to  Trim  Flapping 


COLLECTIVE  PITCH  AT 


COLLECTIVE  PITCH  AT  75R, 


.020 
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(d)  Cyclic  Pitch  to  Trim  Flapping 


FLAPW1SE  BENDING  MOMENT  AT  65  RADIUS  ,MBF65R  ,  IN.- LB 


(a)  Flapwlse  Banding  Moment 


Figure  33.  Range  of  Blade  Response  vs.  Lift  Coefficient-Solidity  Ratio 
for  Model  Blade  Configurations;  Vg  ■  12b  kn,  u  ■  0.3,  ag  * 

0  deg,  alg  -  b  lg  «  0  deg. 
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ROTOR  LIFT  COEFFICIENT  -  SOLIDITY  RATIO  ,  CL/<r 


(b)  Chordwlse  Bending  Maeent 


Figure  33.  Continued. 
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ROTOR  LIFT  COEFFICIENT  -  SOLIDITY  RATIO  ,  CL /<r 


(c)  Torsional  Moment 


Figure  33.  Concluded. 
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EQUIVALENT  FLAPWISE  STRESS 


EQUIVALENT  FLAPWISE  STRESS  ,  LB/ IN. 
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(a)  CONTINUED. 


Figure  3b.  Continued. 
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EQUIVALENT  TORSIONAL  STRESS  ,  LB /IN. 


0  40  10  120  ISO  ZOO  Z40  ZIO  320  3S0  400 

HZIMUTM.  PSI .  OEG 


(a)  Condition  31*,  V  «  0.352,  CL/o  -  0.0790,  C D/o  ■  0.00222 


Figure  35.  Average  Revolution  Time  History  Data;  .65R  Flapvise  Bending, 
91  *  0  deg,  ■  0  deg,  S£  ■  1. 
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(b)  Condition  25,  V  ■  O.VO,  C^/a  -  0.07®9,  Cp/o  »  -0.00113 


Figure  35>  Continued. 


AZIMUTH.  MSI .  OCG 


25,  U  -  0.399,  CL/o  -  0.079>»,  Cp/ 


Figure  36.  Continued. 


AZIMUTH 


(e)  Condition  1*9,  w  *  0.502, 


Figure  36. 
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AZIMUTH.  PCI.  OEG 


(b)  Condition  25 »  w  ■  0.399,  CL/a  -  0.082U,  CD/a  «  -0.0011*2 


Figure  37.  Continued. 
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(c)  Condition  49,  »  -  0.502,  CL/o  -  O.O826,  CD/o  -  0.00545 


AZIMUTH.  PS  I.  DEC 


(a)  Condition  3*»»  y  *  0.332,  C^/o  ■  0.0775*  C^/o  ■ 


Pigure  38.  Average  Revolution  Tine  History  Data;  .65R 
Bending,  0,  •  0  deg,  «_  ■  0  deg,  Sw  -  3. 
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(b)  Condition  25,  v-  -  0.399,  Cj/o  ■  0.0752,  Cp/o  ■ 


0.00153 


Figure  38.  Continued. 
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Figure  38.  Concluded. 


(a)  Condition  3kt  u  -  0.352,  CL/o  -  0.0790,  Cp/o  ■  0.00222 


Figure  39.  Average  Revolution  Time  History  Data;  .65R  Chordwise 
Bending,  01  *  0  deg,  6p  »  0  deg,  Sg  *  1. 
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(b)  Condition  25,  w  ■  O.kO,  CL/  o  ■  0.0789,  Cp/o  ■  -0.00113 
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Flfure  39.  Continued. 
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fUIMUTH.  PSI.  DEG 


(c)  Condition  1»9.  v  -  0.50,  CJ  a  -  0.0835,  Cp/o  »  0.00l»l»$ 


Figure  39.  Concluded. 


219 


.(SR  CMOROWISE  BENDING.  M0C.(5R  .  IN.-LB 


(a)  Condition  3U,  y  ■  0.352,  Cg/o  ■  0.0796,  Cp/o  »  0.00328 


Figure  UO.  Average  Revolution  7t»e  Hi*tory  Data;  .65R  Chordvise 
Bending,  0^  ■  -8  deg,  6^  ■  0  deg,  Sg  »  1. 
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AZIMUTH.  PCI.  DCS 


(c)  Condition  1»9*  u  -  0.502,  C L/o  -  0.0810,  Cp/o  ■  0.00t33 


400 


Figure  to.  Concluded. 


(a*  Condition  3fc,  u  -  0.352,  C Jo  -  0.0795,  Cp/o  -  0.0W31 


Figure  1*1.  Average  Revolution  Tiae  History  Data;  .65R  Chordvlse 
Bending,  0^  ■  0  deg,  6p  ■  5  deg,  SE  ■  1. 
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40  0C  IfO  1(0  COO  C40  200  ICO  3(0 

AZIMUTH.  PSI.  KB 


(b)  Condition  25,  U  ■  0.399,  CL/c  -  0.082U,  C p/o  ■  -0.0011*2 


Figure  1*1.  Continued. 


(c)  Condition  1*9,  u  B  0.502,  CL/0  «  0.0826,  Cp/o  ■  0.0051*5 


Figure  1*1.  Concluded. 
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IN-LB 


(a)  Condition  3*».  w  ■  0.352,  C Jo  ■  0.0775,  C D/o  ■  0.00157 


Figure  U2.  Average  Revolution  Tiae  History  Data;  .65R  Chordvise 
Bending,  ■  0  deg,  6f  ■  0  deg,  Sg  ■  3. 
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. S5R  CHOROWI 


(c)  Condition  1*9»  W  ■  0.502,  C^/o  ■  0.0822,  C^/o  ■  O.OOVfl* 


Figure  1*2.  Concluded. 
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0  40  10  ICO  ICO  200  240  200  320  ICO  4C 

AZIMUTH.  PSI.  DEG 


(a)  Condition  3*»,  u  -  0.352,  CL/o  ■  0.0790,  C^/a  «  0.00222 


Figure  1»3.  Average  Revolution  Tine  History  Data;  .20R  Torsion, 


8,  a  0  deg,  *6j;  »  0  deg,  Sp  *  1. 


.20R  TORSION.  HT.20R  .  IN-LB 
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1(1 _ l _ I _ _ 1 - 1 - 1 - 1 - 1 - 1 - 1 

0  40  10  120  ISO  200  240  210  320  3S0  40) 

rtZIMUTH.  PS I,  DEG 


(b)  Condition  25,  U  «  0.1)0,  CL/o  «  0.0789,  CD/a  »  -0.00113 


Figure  1)3.  Continued. 
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AZIMUTH.  PSf.  0E6 


(c)  Condition  1*9.  v  •  0.50,  CL/o  -  0.0835,  Cp/o  »  0. 001*1*5 


Figure  1*3.  Concluded. 


AZIMUTH,  PSI.  DEC 


(a)  Condition  3**.  u  -  0.352,  CL/o  ■  0.0796,  Cp/o  -  0.00328 


Figure  UU.  Average  Revolution  Tine  History  Data;  .20R  Torsion 
•,  ■  -8  deg.  d_  *  0  deg,  S_,  ■  1. 


(c)  Condition  1*9,  v  •  0.502,  CL/o  ■  0.0810,  Cp/e  ■  0.001*33 


Figure  1*1*.  Concluded. 
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(a)  Condition  35,  W  ■  0.352,  CL/o  ■  0.0795,  C D/o  *  0.00131 


figure  U5.  Average  Revolution  Tioe  History  Data;  .20R  Torsion, 
6X  ■  0  deg,  fip  ■  5  deg,  S£  «  1. 
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-Jfl - 1  I _ l  l  1  i 

0  40  00  ttl  1(0  too  (40 

AZIMUTH.  PCI.  KB 


NO  120  ICO 


(b)  Condition  25,  u  ■  0.399,  CL/o  ■  O.O82I4,  Cp/o  -0.001U2 


Figure  U5 .  Continued. 
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(c)  Condition  1*9,  u  -  0.502,  CL/o  -  0.0026,  C^/a  ■  0.005l»5 


Figure  1*5 .  Concluded. 
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40  10  120  ICO  too  240  210  320  SCO 


AZIMUTH.  PS  I.  KG 


(a)  Condition  3U,  u  »  0.352,  CL/o  ■  0.0775,  CD/o  -  0.00157 


Figure  U6.  Average  Revolution  Tins  History  Data;  .20R  Torsion, 
6j^  ■  0  deg,  4p  ■  0  deg,  Sg  ■  3. 
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0  40  M  ItO  160  tOO  140  tM  J?0  160  40l 

AZIMUTH.  PSI.  DEG 


(c)  Condition  1*9,  i*  *  0.502,  C^/o  ■  0.0822,  C^/o  «  O.OOl*7l* 


Figure  1*6.  Concluded. 
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Coefficient 


SECTION  MOMENT  COEFFICIENT, 


(b)  Moment  Coefficient 


Figure  U7.  Concluded. 


2 1*2 
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and  Model  Data. 


ANGLE  OF  ATTACK  ,  a  ,  DEG 


2h6 


Figure  U8.  Continued 


2k9 


Figure  1*9.  Continued. 


Figure  1»9.  Continued. 


Figure  U9 .  Concluded . 


•*  '£W 


Figure  50.  S*nple  of  9te»dy  end  Unsteady  Aerodynemic  Regions  of 
Application  for  the  Flapped  Blade. 


253 


Figure  >2.  Cc-rtinued . 


- 1 


Figure  52.  Concluded. 


LATERAL  CYCLIC  TORQUE  COEFFICIENT  -  L  IFT  COEFFICIENT  - 

PITCH,  A,]  ,  OEG  SOLIDITY  RATIO  ,  CQ/<r  SOLIDITY  RATIO  ,  C, /<r 


Figure  53.  Correlation  of  Test  and  Calculated  Data  for  Condition  7; 

Vg  *  1?U  kn,  u  =  0.30,  ag  ■  0  det ,  a^g  =  b  =  0  deg. 
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CHORDWISE  BENDING  “LAPWISE  BENDING 

MOMENT  ,Mbc,  IN.- L3  MOMENT ,  M#F  ,  IN. -LB 


a  a 


DRAG  COEFFICIENT-  LIFT  COEFFICIENT - 

K.IOITV  RATIO.  C_/,  SOLIDITY  RATIO.  C./v 


0*0  0 


LATERAL  CYCLIC 

PITCH  ,  A„  , OEG  SOLIOITY  RATIO.CpM/* 


C HOMOWISE  BENOING  FLATWISE  BENOING 

MOMENT,  M.e,  IN  *  LB  MOMENT ,  M-.  ,  IN- LB 


O  tot 

□  UNIFORM  INFLOW 

a  VARIABLE  NTLOW,  UNSTEADY  ACNO  SCALING  CMC  9 
A  VARIABLE  INFLOW,  UNSTCAOT  AERO  SCALING  CASE  I 


Figure  51*.  Concluded. 
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LATERAL  CYCLIC  TOROUE  COEFFICIENT  -  LIFT  COEFFICIENT* 

PITCN  .  A„  ,  OEG  SOLIOITY  RATIO  .  CQ/r  SOLIOITY  RATIO  .  CL/r 


LATERAL  CYCLIC  TORQUE  COEFFICIENT  -  LIFT  COEFFICIENT  - 

PITCH  .  A,.  .  DEG  SOLIDITY  RATIO  ,  C Q/<r  SOLIDITY  RATIO  ,  C,  /<r 


Figure  56.  Correlation  of  Test  and  Calculated  Data  for  Condition  80; 

V  »  1PU  kn,  u  *  0.30,  a  ■  -5  def,  a,  »  b.  •  0  deg. 

8  S  18  18 
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Figure  56.  Concluded. 
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LATERAL  CYCLIC  TOROUE  COEFFICIENT  -  LIFT  COEFFICIENT  - 

PITCH  .  A,«  .  OEG  SOLIOITY  RATIO.  C0 SOLIDITY  RATIO  ,  Ct  fv 
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Figure  57.  Correlation  of  Test  and  Calculated 

V  *  12U  kn,  u  ■  0.30,  a  ■  -5  deg, 
s  s 
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J5R  FlHPHISE  BENOING. 

MBF.15R  .  IK-LB  ?0R  FlARWIBC  KNDIN6.  M0f  !0«  .  iNrLi 


(a)  .20R  Flapwise  Bending,  Condition  7 


(b)  .35R  Flapwise  Bending,  Condition  7 


Figure  59.  Correlation  of  Test  and  Calculated  Flapwise  Bending  With 
Uniform  Inflow;  0^^  ■  0  deg,  6p  ■  0  deg,  S£  «  1. 
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(d)  .0OR  Flapvise  Bending,  Condition  7 


Figure  59.  Continued. 


FLflPMtSC  BENOING  20R  FLRPMISE  BENDING 

«BF  4JR  .  IN-LB  HBF  20R  .  IN.-LB 
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(e)  .20R  Flapvise  Bending,  Condition  36 


(f)  .*oR  Flapvise  Bending,  Condition  36 


Figure  59.  Continued. 
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(g)  .65H  Flapwise  Bending,  Condition  36 


(h)  .35R  Flapwise  Bending,  Condition  71 


Figure  59.  Continued. 
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IVNI  '  MS!  j m,  'ONION 98  XIHriWU  lit  *  M*  J0M  'ONION 38 


U 


(i)  .65R  Flapvise  Bending,  Condition  71 


(j)  .80R  Flapvise  Bending,  Condition  71 


Figure  59.  Continued. 
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(b)  .65R  Flapvise  Bending*  Condition  71 


Figure  60.  Correlation  of  Test  and  Calculated  Flapvise  Bending  With 
Variable  Inflow;  8.  ■  0  deg,  f_  ■  0  deg,  Sr  «  1. 


*5R  FLRPNISE  BENDING.  HBF  EJR  .  IN.-LB  *5R  FLRPHISE  BENOING.  MBF.S5R  .  INrLB 


(c)  .65R  Flapwise  Bending,  Condition  71 


(d)  .65R  Flapwise  Bending,  Condition  80 


Figure  6l. 


Continued. 
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(e)  .65R  Flapvise  Bending,  Condition  82 


Figure  6l.  Concluded. 
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S5R  FLRPW1SE  BENOING.  MBF.J5R  .  IN.-LB  .S5R  FLRPMISE  BENDING.  MBF.S5R 


CALCULATED 


(a)  .65R  Flapvise  Banding,  Condition  36,  Unsteady  Aerodynamic 

Scaling  Case  3 


RZINUTH.  PS  I.  DEG 


(b)  .65R  Flapvise  Bending,  Condition  36,  Unsteady  Aerodynamic 

Scaling  Case  1 


Figure  62.  Correlation  of  Test  and  Calculated  Flapvise  Bending  With 
Variable  Inflow;  6.  -  -8  deg,  5„  «  0  deg,  S  ■  1. 


S5R  FLflPHISE  BENOING.  MBF  C5R  ,  IN-LB 


(c)  .65R  Flapvise  Bending,  Condition  71,  Unsteady  Aerodynamic 

Scaling  Case  1 


Figure  62,  Concluded. 
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.(SR  FLflPW I SE  BEJOING.  M8F.S5R  .  IN.-LB  *S»  ElRRWIU  KN0IM6.  Oft  .  tNrLl 


(a^  .65R  Flapwise  Bending,  Condition  7 


XX I _ I _ I _ 1 -  - 1 - 1 - 1 - 1 - 1 - 1 - 1 

0  40  10  120  ISO  200  240  210  320  3(0  400 


RZINUTM,  PS  I,  OEG 


(b)  .65R  Flapwise  Bending,  Condition  36 


Figure  63.  Correlation  of  Test  and  Calculated  Flapwise  Bending  With 
Uniform  Inflow;  0=0  deg,  6„  ■  5  deg.  S_  ■  1. 
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AZIMUTH.  PSl.  OCG 


(b)  .65R  Flapvise  Bending,  Condition  71 


Figure  6i».  Correlation  of  Test  and  Calculated  Flapvise  Bending  With 
Variable  Inflow;  0=0  deg,  6„,  ■  5  deg,  S„  ■  1. 
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Figure  65 


40  00  120  1(0  200  240  2(0 

AZIMUTH,  RSI.  DEG 


(b)  .65R  Flapvise  Bending,  Condition  36 


Correlation  of  Tent  and  Calculated  Flapvise  Bending  with 
Uniform  Inflow;  6,  ■  0  deg,  6p»  0  deg,  S£  «  3. 
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SSR  FLAPMISE  BEN01N6.  MBF.C5R  .  tN.-LB  65R  FLRPHISE  BENDING. 
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(a)  Uniform  Inflow 


PZIHUTH.  PSI.  0E6 


(b)  .3  Chord  Vortex  Core  Size 


Figure  67.  Effect  of  Variable  Inflow  Vortex  Core  Size  Assumption  on 
Correlation  of  Test  and  Calculated  .65R  Flapvise  Bending 
0^  ■  0  deg,  fip  ■  0  deg,  S£  ■  1,  Condition  1. 
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35R  CHOROUISE  BENOING.  NBC  J5R  .  I Nr LB  20R  CHOROMISE  BENDING.  MBC.20R  .  I Nr LB 


(a)  .20R  Chordvise  Bending,  Condition  7 


(b)  .35R  Chordvise  Bending,  Condition  7 


Figure  68.  Correlation  of  Test  and  Calculated  Chordvise  Bending  With 
Uniform  Inflow;  0^  «  C  deg,  *  0  deg,  S£  ■  1 
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CALCULATED 


(e)  .20R  Chorawise  Bending,  Condition  36 


— 


psi.  deg 


(f)  ,1»5R  Chordvise  Bending,  Condition  36 


Figure  68.  Continued. 


ISC  BENOING.  MBC.tOft  .  IN.-LB  .C5R  CHOROMISE  BENOING.  HBC.85R  .  INrLB 


<0 


(I)  .65R  Chordwise  Bending,  Condition  71 


0  40  00  120  1(0  200  240  200  120  3(0  400 

BZINUTH.  PSI.  KG 


(j)  ,80R  Chordwise  Bending,  Condition  71 


Figure  68.  Continued. 
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AZIMUTH.  PSI.  DEG 


(l)  .65R  Chordvise  Bending,  Condition  82 


Figure  68.  Concluded 


65R  Chordvise  Bending,  Condition  36 


nzihuth.  ps I,  DEB 


(b)  .65R  Chordvise  Bending,  Condition  71 


Figure  69.  Correlation  of  Test  and  Calculated  Chordvise  Bending  With 
Variable  Inflov;  0  »  0  deg,  6„  *  0  deg,  SF  »  1. 


Figure  70.  Concluded. 


G5R  CHOROHISE  BENDING.  MBC.S5R  .  IN.-LB  S5R  CHOROHISE  BENDING.  HOC  t5R 


(a)  .65R  Chordvise  Bending,  Condition  36,  Unsteady  Aerodynamic 

Scaling  Case  3 


AZIMUTH,  RSI.  DEG 


(b)  .o5R  Chordvise  Bending,  Condition  36,  Unsteady  Aerodynamic 
Scaling  Case 


Figure  71.  Correlation  of  Test  and  Calculated  Chordvise  Bending  With 
Variable  Inflow;  fl,  ■  -8  deg,  6„  «  0  deg,  S  «  1. 
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(c)  .65R  Chordvise  Bending,  Condition  71»  Unsteady  Aerodynamic 

Scaling  Case  1 


.(SR  CHOROMISC  BENDING.  H9C.(5R  .  IK  .9  .(SR  CHOROWISE  BENDING.  NBC  (SR  .  INrLB 


(a)  .65R  Chordvlae  Bending,  Condition  7 


0  40  to 


110  1(0  tOO  t40  710  170  1(0  400 

RZINUTM.  PSI,  OEG 


(b)  .65R  Chordwise  Bending,  Condition  36 


Figure  72.  Correlation  of  Test  and  Calculated  Chordvise  Bending  With 
Uniform  Inflow;  6^0  deg,  6^  *  5  deg,  Sg  «  1. 
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(a)  .65R  Chordwise  Bending,  Condition  36 


(c)  .65R  Chordwise  Bending,  Condition  71 


(d)  .65R  Chordwise  Bending,  Condition  8 2 


Figure  7I*.  Concluded. 
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AZIMUTH,  PSl.  DEG 


(b)  .3  Chord  Vortex  Core  Size 


Figure  76.  Effect  of  Variable  Inflow  Vortex  Core  Size  Assumption  on 
Correlation/Test  and  Calculated  .65R  Chordvise  Bending; 

0  ■  0  deg,  6-  ■  0  deg,  S„  «  1,  Conditionl. 


TEST 

CALCULATED 


LAGGING  AT  HINGE.  ZETA  .  DEG  8LA0E  FLAPPING  AT  HINGE.  BETA 


(a)  Blade  Flapping  at  Hinge,  Condition  1 


AZIMUTH.  PS I.  DEG 

(b)  Blade  Lagging  at  Hinge,  Condition  1 


Figure  77.  Correlation  of  Test  and  Calculated  Blade  Hinge  Motions  With 
Uniform  Inflow;  ei  =  0  deg,  5p  ■  0  deg,  S£  a  1. 
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Figure  77.  Continued. 


(h)  Blade  Lagging  at  Hinge ,  Condition  71 


Figure  77.  Continued. 
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LAGGING  AT  HINGE.  ZETA  .  DEG  BLAOE  FLAPPING  AT  HINGE.  BETA  ,  OEG 


AZIMUTH,  psi.  deg 
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(e)  Blade  Flapping  at  Hinge,  Condition  71 
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(f)  Blade  Lagging  at  Hinge,  Condition  71 


Figure  78.  Concluded. 
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(a)  Blade  Flapping  at  Hinge,  Condition  82 


(b)  Blade  Lagging  at  Hinge,  Condition  82 


Figure  79.  Correlation  of  Test  and  Calculated  Blade  Hinge  Motions 

With  Uniform  Inflow;  6  «  -8  deg,  6_*  0  deg,  S  =  1. 

J.  r  E 


318 


)  40  10  It!  ISO  too  140  210  120  1(0  400 

AZIMUTH.  PS  I,  DEG 

(b)  .20R  Torsion,  Condition  36 


Figure  8l.  Correlation  of  Test  and  Calculated  Torsion  With  Uniform 
Inflow;  8.  «  0  deg,  6_.  ■  0  deg,  Sr  «  1. 


HOZ  1U  'NOlSdOl  80Z 


TOWS  I  ON*  NT.tM  .  IN-lB  .20R  TORSION.  MT.20R 


(b)  .20R  Torsion,  Condition  71 


Figure  8?.  Correlation  of  Test  and  Calculated  Torsion  With  Variable 
Inflow;  ei  *  0  deg,  6p  ■  0  deg,  S£  =  1. 
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Figure  83.  Continued. 


20*  TORSION.  HT.20R  .  INrLS  ,t»R  TORSION.  NT. (OR 


(a)  .20R  Torsion,  Condition  36,  Unsteady  Aerodynamic  Scaling  Case  3 


(b)  .20R  Torsion,  Condition  36,  Unsteady  Aerodynsmic  Scaling  Case  1 


Figure  8U.  Correlation  of  Test  and  Calculated  Torsion  With  Variable 
Inflow;  6X  *  -8  deg,  6p  *  0  deg,  »  1. 
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20ft  TORSION.  NT. 20ft  .  IN-LB 


(c) 


•20R  Torsion, 


Condition  71 


t  Unsteady 


Figure  8U.  Concluded 
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•  20R  TORSION,  NT.20R  .  INrLB  .20R  TORSION,  MT.20R  .  IN  LB 


(c)  .20R  Torsion,  Condition  ?1 


(d)  .20R  Torsion,  Condition  82 


Figure  85.  Concluded. 
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AZIMUTH,  psi.  ocg 


(b)  .20R  Torsion,  Condition  71 


Figure  86.  Correlation  of  Test  and  Calculated  Torsion  With  Variable 
Inflow;  6*0  deg,  *  5  deg,  S_  *  1. 
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(b)  .20R  Torsion,  Condition  36 


Figure  87.  Correlation  of  Test  and  Calculated  Torsion  With  Uniform 
Inflow;  0  ■  0  deg.  6-  ■  0  deg,  S„  «  3. 
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(ft)  .20  Torsion,  Condition  36 
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•  «0  10  lit  1(0  200  240  200  120  ISO  400 


RZINUTH.  RSI.  0E6 


(b)  .20  Torsion,  Condition  71 


Figure  88.  Correlation  of  Test  ftnd  Calculated  Torsion  With  Variable 
Inflow;  6:  =  0  deg,  *  0  deg,  S£  =  3. 
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AZIMUTH.  PS  I.  DEG 


Figure  89 


(b)  .3  Chord  Vortex  Core  Size 


.  Effect  of  Variable  Inflow  Vortex  Core  Size  Assumption 
Correlation  of  Test  and  Calculated  .  20R  Torsion;  0  ■ 

6„  =  0  deg,  Sr  »  1,  Condition  1. 


APPENDIX  I 


FULL-SCALE  H-3^  BLADE  DATA 


TABLE  XLIX.  H-3U  MAIN  ROTOR  BLADE 

MISCELLANEOUS  DATA 

PARAMETER 

VALUE 

UNITS 

Airfoil  section  NACA  0012 

Airfoil  chord 

16. it 

in. 

Radial  location  of  airfoil  Inboard  end 

53.7 

in. 

Elastic  axis  location  aft  of  leading  edge 
(for  radial  stations  outboard  of  5**  in.) 

3.9 

in. 

Rotor  radius 

336.0 

in. 

Flapping  and  lead-lag  hinge  radial  location 

12.0 

in. 

Young's  modulus  (E ) 

10xl06 

lb/in? 

Shear  Modulus  (C) 

ltxlO6 

lb/in? 

Feathering  axis  location  aft  of  leading  edge 

k.l 

in. 
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TABLE  L.  H-31*  MAIN  ROTOR  BLADE  DISTRIBUTED  WEIGHT,  CHORDWISE  CENTER - 
OF-GRAVITY  POSITION,  AND  TORSIONAL  INERTIA  (STANDARD  BLADE) 


*  Distance  of  center  of  gravity  forward  of  feathering  axis. 


TABLE  LI.  H-3U  MAIN  BDTOR  BLADE  DISTRIBUTED  WEIGHT,  CHOHDWISE 
CENTER-OF-GRA. VTTY  POSITION,  AND  TORSIONAL  INERTIA 
(AS  INSTRUMENTED  FOR  PRESSURE  MEASUREMENTS) 

Ar 

(in) 

rtn 

(in) 

rout 

(in) 

V 

(lb/in) 

Ycg* 

(in) 

ie 

(lb-sec2) 

1.00 

12.00 

13.00 

9.00 

0.00 

0.270 

11.00 

13.00 

24.00 

5.50 

0.00 

0.129 

4.00 

24.00 

28.00 

4.75 

0.00 

0.068 

4.00 

28.00 

32.00 

1.64 

0.00 

0.043 

2.00 

32.00 

34.00 

1.57 

0.00 

0.034 

2.00 

34.00 

36.00 

1.29 

0.00 

0.024 

2.00 

36.00 

38.00 

1.07 

0.00 

0.015 

3.00 

38.00 

41.00 

0.45 

0.00 

C.008 

13.00 

41.00 

54.00 

0.30 

0.00 

0.005 

23.00 

54.00 

77.00 

0.35 

-0.55 

0.012 

1.00 

77.00 

78.00 

0.40 

0.20 

0.014 

2.00 

78.00 

80.00 

0.50 

0.20 

0.018 

3.00 

80.00 

83.00 

0.45 

0.20 

0.016 

7-00 

83.00 

90.00 

0.55 

0.20 

0.018 

36.00 

90.00 

126.00 

0.45 

0.20 

0.016 

13.00 

126.00 

139-00 

0.55 

0.20 

0.018 

35.00 

139.00 

174.00 

0.45 

0.20 

0.016 

12.00 

174.00 

186.00 

0.55 

0.20 

0.019 

60.00 

186.00 

246.00 

0.45 

0.2C 

0.016 

11.00 

246.00 

257.00 

0.55 

0.20 

0.019 

24.00 

257-00 

& 31.00 

0.45 

0.20 

0.016 

24.00 

281.00 

305.00 

0.55 

0.20 

0.019 

12.00 

305.00 

317.00 

0.45 

0.20 

0.016 

6.00 

317.00 

323.00 

0.55 

-0.27 

0.018 

4.00 

323.00 

327.00 

0.80 

-0.01 

0.028 

2.00 

327.00 

329.00 

0.95 

-0.01 

0.034 

2.00 

329.00 

331.00 

0.45 

-0.97 

0.016 

2.00 

331.00 

333.00 

0.90 

0.0 

0.031 

2.00 

333.00 

335.00 

0.15 

-2.3 

0.0050 

1.00 

335.00 

336.00 

0.25 

-2.3 

0.0080 

*  Distance  of  center  of  gravity  forward  of  feathering  axis. 


TABLE  LI I .  H-3I4  MAIN  ROTOR  BLADE  STRUCTURAL  AREA  AND  MOMENTS  OF 

INERTIA 

Ar 
( io) 

,rln. 

( in. ) 

rout 
(  In. ) 

A® 

(in?) 

('% 

!flp 
(  in!*) 

1.50 

12.0C 

13.50 

78.33 

81.76 

7.80 

2.50 

13.50 

16.00 

55.00 

50.00 

7.80 

2.00 

16.00 

18.00 

55.00 

17.80 

50.00 

8.00 

18.00 

26.00 

53.33 

U6.20 

50.00 

1.50 

26.00 

27-50 

It?.  50 

29.20 

29.20 

1.50 

27.50 

29.00 

26.77 

9.20 

11.  lit 

3.50 

29.00 

32.50 

16.30 

3;  .00 

23.10 

2.50 

32.50 

35.00 

lit.  70 

30.35 

15.83 

1.00 

35.00 

36.00 

12.90 

26.33 

10.87 

It. 00 

36.00 

Uo.oo 

5.59 

23.66 

7.07 

It. 00 

U0.00 

It  It. 00 

3.59 

l8.3»t 

1.90 

9.00 

It  It.  00 

53.00 

3.k5 

17-91 

1.85 

2  It.  00 

53.00 

77.00 

3.21 

17.18 

1.71 

259.00 

77.00 

336.00 

2.92 

15.91* 

l.it? 
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TABLE  LIII.  H-3l*  MAIN  ROTOR  BLADE  STRUCTURAL  AREA 
TORSIONAL  RADIUS  OF  GYRATION 


At 

(ia) 

finx 

(ia) 

rout 

(ia) 

res 
( ia) 

l.GC 

12  00 

13.00 

3.41 

n.oo 

13.00 

24.00 

3.02 

4.00 

24.00 

28.00 

2.35 

4.00 

28.00 

32.00 

3.20 

2.00 

32.00 

34.00 

2.88 

2.00 

34.00 

36.00 

2.71 

2.00 

36.00 

38.00 

2.35 

3.00 

38.00 

41.00 

2.50 

13.00 

41.00 

54.00 

2.41 

23.00 

54.00 

77.00 

2.43 

1.00 

77.00 

78.00 

2.46 

2.00 

78.00 

80.00 

2.46 

3.00 

80.00 

83.00 

2.45 

7.00 

83.00 

90.00 

2.37 

36.00 

90.00 

126.00 

2.45 

13.00 

126.00 

139.00 

2.37 

178.00 

139.00 

317.00 

2.45 

6.00 

317.00 

323.00 

2.37 

13.00 

323.00 

336.00 

2.45 
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APPENDIX  II 


COMPUTER  PROGRAM  DOCUMENTATION 


USER’S  INSTRUCTIONS  FOR  DATA  RETRIEVAL  AID  DISPLAY  PROGRAM 
(Refer  to  the  report  section  "TEST  DATA  OBTAINED" ,  p.30. 
Input 

A.  Tape  for  one  of  the  blade  configurations 

B.  Input  Cards 

(l)  First  Card 


Field 

Columns 

Input  Item: 

1 

1-1* 

Operating  condition  number 

2 

5-8 

Physical  measurement  (channel)  number 

3 

9-12 

First  revolution  number  in  a  time  history  saaple 

k 

13-16 

Last  revolution  number  in  a  time  history  sample 

5 

17-20 

First  revolution  in  a  time  history  sample  for 

21-21* 

frequency  analysis 

6 

Last  revolution  in  a  time  history  sample  for 

25-31* 

frequency  analysis 

7 

Lowest  frequency  component  in  analysis 

8 

35-1*1* 

Frequency  increment  in  analysis 

9 

1*5— 1*8 

Humber  of  frequencies  in  analysis 

10 

1*9-52 

Control  vord  for  summary  option  (mean,  rms,  max,  min) 

11 

53-56 

Control  vord  for  azimuth-averaged  revolution  time 

57-60 

history 

12 

Number  of  harmonics  in  harmonic  analysis 

13 

61-61* 

Harmonic  table  control  number 

Fields  7  and  8  are  real  positive  decimal  numbers .  All  other  fields  are 
right-adjusted  integers.  Fields  3  and  5  oust  be  nunbers  from  1  to  1*0  to 
exercise  the  corresponding  options.  Fields  1*  and  6  must  be  equal  to  or 
greater  than  3  and  3  respectively,  and  no  greater  than  1*0  to  exercise 
the  corresponding  options.  Fields  10,  11  and  13  oust  be  1,  2  or  3  to 
exercise  those  options.  If,  on  the  first  card,  any  of  the  fields  3,  5,  10, 
11,  12,  13  contain  999  or  are  blank,  the  corresponding  option  vill  not  be 
exercised  for  that  condition-channel  combination.  The  options  may  be 
turned  on  for  succeeding  condition-channel  combinations  by  inputs  in  those 
fields. 


(2)  Succeeding  Cards 

The  foraat  is  the  same  as  the  first  card  format.  If  any  fields  are  0  or 
blank,  the  inputs  vill  remain  as  they  were  for  the  preceding  card.  If  any 
of  fields  3,  5,  10,  11,  12,  and  13  contain  999,  that  option  vill  not  be 
exercised. 
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If  condition  number  9999  is  input,  the  body  of  data  to  be  retrieved  is 
considered  to  be  complete. 

Input  for  a  new  body  of  data  can  now  be  entered,  otherwise  the  computer 
run  will  be  terminated. 

A  maximum  of  100  request  cards  can  be  input  for  each  body  of  data.  If 
more  are  desired  in  a  computer  run,  9999  can  be  entered  a*i  a  condition 
number,  and  a  new  body  of  data  started. 

Output 

Time  history  stapling  will  be  carried  out  in  response  to  legal  entries  in 
fields  3  and  U.  The  indicated  revolutions  will  be  tabulated  in  vertical 
columns.  The  table  heading  will  contain  the  blade  configuration  descrip* 
tion,  condition  number,  and  rotor  operating  parameters.  A  one-page 
printer  plot  of  the  time  history  data  will  appear  for  each  revolution. 

The  condition  number  and  rotor  operating  parameters  will  be  printed  be¬ 
side  the  plot. 

Frequency  analysis  will  be  carried  out  in  response  to  proper  entries  in 
fields  5,  6,  7,  8,  and  9>  The  amplitude  consents  A.,  B.,  and  R.  will 
be  tabulated  and  plotted  for  each  frequency.  The  nonharmonic  analysis 
frequencies  should  range  between  .15  and  15.  cycles  per  revolution  if 
all  40  revolutions  of  data  are  used.  Ah  input  of  1,  2,  or  3  in  field 
10  will  print  out  average,  root-mean-square  about  the  average,  maximum 
value,  and  minima  value  in  the  data  for  that  condition  and  channel. 

The  tabulation  ia  continued  by  specifying  additional  conditions  and/or 
channels  on  succeeding  cards.  If  1  is  input  in  field  10,  the  tabulation 
will  proceed  by  channels  for  a  condition.  The  tables  will  be  output  in 
ascending  condition  number.  The  channels  will  appear  in  the  same  order 
as  they  appeared  in  the  input  deck  for  the  specified  body  of  data.  If  3 
is  input  in  field  10,  the  tabulation  will  proceed  by  conditions  for  a 
channel.  The  tables  viU  be  output  in  ascending  channel  number.  The 
conditions  will  appear  in  the  same  order  as  they  were  entered  in  the 
input  deck.  If  2  is  input  in  field  10,  both  table  types  will  be  output. 

An  input  in  field  11  will  print  out  average  revolution  time  histories. 

The  tabulation  of  similar  data  for  succeeding  channels  and  conditions 
will  proceed  by  channels  for  a  condition  if  1  is  placed  in  field  11. 

Hie  tabulation  will  proceed  by  conditions  for  a  channel  if  3  is  placed 
in  field  11.  A  2  in  field  11  will  cause  the  printing  of  both  table  types. 
Time  history  printer  plots  of  each  average  revolution  will  also  be  printed. 

Similarly,  inputs  in  fields  12  and  13  provide  a  harmonic  analysis  which 
is  tabulated  by  channels  for  a  condition  if  a  1  is  input  in  field  13,  by 
conditions  for  a  channel  if  3  is  input,  and  both  tables  if  2  is  input. 


LIST  OF  SYMBOLS 


AS 

Al. 

A1.R 

*i 

a(-j) 

*1. 


B 

bj 


B 


IsR 


b 


bU'j) 


nondimens ional  pitch  rate 
flexibility  matrix 

indicated  angle  of  attack  for  test  point  (i),  deg 

even  (cosine)  part  of  frequency  component  from  numerical  fre¬ 
quency  analysis,  engineering  units 

cosine  part  of  nonharmonic  amplitude  component,  engineering 
units 

2 

blade  spar  structural  crossection  area,  in. 

2 

tunnel  test  section  crossectional  area,  ft 

lateral  cyclic  pitch,  blade  pitch  dovn  over  tail,  deg 

indicated  value  of  A^s,  deg 

cosine  part  of  ith  modal  response  to  unit  cosine  force,  in. 
analytical  counterpart  of  Aj 

longitudinal  first-harmonic  flapping  component,  blade  dovn 
over  tail,  deg 

r.ondimens  ional  pitch  acceleration 

odd  (sine)  part  of  frequency  component  from  numerical  fre¬ 
quency  analysis,  engineering  units 

sine  part  of  nonharmonic  amplitude  component,  engineering  units 

longitudinal  cyclic  pitch,  blade  pitch  dovn  over  right 
(advancing)  side,  deg 

indicated  value  of  ,  deg 

number  of  rotor  blades 

sine  part  of  1th  modal  response  to  unit  cosine  force,  in. 

lateral  first-harmonic  flapping  component,  blade  dovn  over 
right  side,  deg 

analytical  counterpart  of 

section  coefficient  of  drag 


MS 


'd(i) 


'dR(i) 


V° 


'DP 


'EFF 


'EFFT 


'EFFV 


'EFFX 


'EFFY 


'(i) 


'l(i) 


CL/o 


■  (  i ) 

CmR(i) 

Cmn 

CMREG1 

CMREG2 


'ot 


corrected  coefficient  of  drag  at  the  ith  two-dimensional  test 
point 

uncorrected  coefficient  of  drag  at  the  ith  two-dimensional 
test  point 

rotor  drag  coefficient  -  solidity  ratio,  DR/pwQ2R4o 
coefficient  of  pth  viscous  damper 
effective  rotor  head  damping 
torsional  damping  coefficient,  in. -lb-sec 
vertical  damping  coefficient,  lb-sec/in. 
fore-aft  damping  coefficient,  lb-sec/in. 
lateral  damping  coefficient,  lb-sec/in. 
indicated  chordwiae  force  at  test  point  (i),  lb 
section  coefficient  of  lift 

coefficient  of  lift  at  the  ith  two-dimensional  teat  point 
rotor  coefficient  of  lift 

rotor  lift  coefficient-solidity  ratio,  i^/pwfl2R4o 
section  coefficient  of  pitching  moment 

corrected  coefficient  of  pitching  moment  at  the  ith  two-dimen¬ 
sional  test  point 

uncorrected  coefficient  of  pitching  moment  at  the  ith  two-dimen¬ 
sional  test  point 

balance  conversion  matrix  elements 

upper  limit  of  local  angle  of  attack  at  which  unsteady  aerody¬ 
namic  pitching  moment  data  is  used  in  blade  response  calcula¬ 
tion,  deg 

lower  limit  of  local  angle  of  attack  at  which  steady  aerodynamic 
pitching  moment  data  is  used  in  blade  response  calculation,  deg 

value  of  Cm  in  unsteady  aerodynamic  data  table 

section  normal  force  coefficient 


546 


CNREG1 


CNREG2 


C  (M) 
nss 

C 

nsst 

Cnt 

CPM/0 

Cq/° 

CRM/o 

CT 

Cy/O 

c 

M 

M 

db 

DR(i) 


(lWj 


BCi 


BFi 


upper  limit  of  local  angle  of  attack  at  which  unsteady  aerody¬ 
namic  pitching  moment  data  is  used  In  blade  response  calcula¬ 
tion,  deg 

lower  limit  of  local  angle  of  attack  at  which  steady  aerody¬ 
namic  pitching  moment  data  ia  used  in  blade  response  calcula¬ 
tion  deg, 

value  of  steady  state  C  at  Mach  number  M  and  a  (M) 

J  n  sn 

value  of  steady  state  C  at  a  .in  unsteady  aerodynamic  data 
table  n  8nt 

value  of  in  unsteady  aerodynamic  data  table 
rotor  pitching  moment  coefficient-solidity  ratio,  Mg/pnn2R5o 
rotor  torque  coefficient-solidity  ratio,  -Ng/pwfl2R5o 
rotor  rolling  moment  coefficient  -  solidity  ratio  Lg/pnfi2R5o 
rotor  thrust  coefficient,  -Zs/pxft2R4 

rotor  side-force  coefficient  -  solidity  ratio,  Yg/p*n2R4o 

airfoil  chor*'. ,  in.  or  ft 

physical  damping  matrix 

generalized  damping  matrix 

rotor  drag  force,  lb 

indicated  drag  force  at  test  point  (i),  lb 

2 

Young's  modulus,  lb/in. 

arbitrary  offset  constant,  physical  units 

rotor  force  factor,  lb 

Jth  analysis  frequency,  cycles/rev 

initial  analysis  frequency,  cycles/rev 

experimental  nonrotating  natural  frequency  for  the  ith  chord- 
wise  mode,  cps 

experimental  nonrotating  natural  frequency  for  the  ith  flap- 
wise  mode,  cps 


fTi  experimental  nonrotating  natural  frequency  for  the  ith  tor¬ 

sional  mode,  cps 
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GA(  i ) 


AE 


AO 


AS 


DE 


JD(i) 


DO 


DS 


LE 


JL(i) 


LO 


LS 


ME 


JM(i) 


°MS 

H 

»L 

«N 


H 


x 


2 

shear  modulus,  lb/ln. 

angle  of  attack  transducer  reading  at  ith  two-dimensional  test 
point,  gage  units 

angle  of  attack  transducer  end  zero,  gage  units 

angle  of  attack  transducer  average  zero,  gage  units 

angle  of  attack  transducer  start  zero,  gage  units 

balance  drag  component  end  zero,  gage  units 

balance  drag  component  reading  at  ith  two-dimensional  test 
point,  gage  units 

balance  drag  component  average  zero,  gage  units 
balance  drag  component  start  zero,  gage  units 
balance  lift  component  end  zero,  gage  units 

balance  lift  component  at  ith  two-dimensional  test  point,  gage 
units 

balance  lift  component  average  zero,  gage  units 

balance  lift  component  start  zero,  gage  units 

balance  p?tch  moment  component  end  zero,  gage  units 

balance  pitch  moment  component  at  ith  two-dimensional  test 
point,  gage  units 

balance  pitch  moment  component  average  zero,  gage  units 

balance  pitch  moment  component  start  zero,  gage  units 

2 

total  pressure,  lb/ft 

aerodynamic  plus  gravity  tare  moment  in  L-..  direction,  gage 
units 

aerodynamic  plus  gravity  tare  moment  in  direction,  gage 
units 

aerodynamic  plus  gravity  tare  moment  in  N_.c  direction,  gage 
units 

aerodynamic  plus  gravity  tare  force  in  force  direction, 
gage  units 


34S 


^ 


H 


Y 


{>} 

[*] 

:chd 

Jflp 

^mch 

Imfl 

xe 

j 

K 


"r 

Vi) 


aerodynamic  plus  gravity  tare  force  Y^^,  force  direction, 
gage  units 

aerodynamic  plus  gravity  tare  force  in  Z force  direction, 
gage  units 

unit  column  matrix 
unit  diagonal  matrix 

U 

chordvise  blade  bending  stiffness  area  moment  of  inertia,  in. 

U 

flapwise  blade  bending  stiffness  area  moment  of  inertia,  in. 

chordvise  mass  moment  of  inertia  of  blade  per  unit  span,  lb-sec 

2 

flapwise  mass  moment  of  inertia  of  blade  per  unit  span,  lb-sec 
torsional  mass  moment  of  inertia  of  blade  per  unit  span,  lb-sec 
blade  torsional  stiffness  constant,  in. 

number  of  digital  data  values  in  a  time  history  segment  for 
frequency  analysis 

stiffness  matrix 

angle  of  attack  transducer  conversion  constant,  deg'1 

effective  stiffness,  lb/in. 

torsional  spring  constant,  in. -lb 

vertical  spring  constant,  lb/in. 

fore-aft  spring  constant,  lb/in. 

lateral  spring  constant,  lb/in. 

p 

pressure  transducer  conversion  constant,  lb/ft 

indicated  right  side  down  rolling  moment  in  rotor  shaft  axis 
system,  gage  units 

static  offset  correction  for  L^^,,  gage  units 
rotor  lift  force,  lb 

indicated  lift  force  at  ith  two-dimensional  test  point,  ib 

rotor  right  side  down  rolling  moment  in  shaft  axis  system, 
ft -lb 


M 


Mach  number 


[«] 

M»(i) 

V, 

mbc.xr 

mbp.xr 

Mc 

MCB( i ) 

mcfs 

mcrm 

v? 

“efft 

^FFV 

“efpy 

«F 

Mns 

mfrm 

M(i) 

MGi 

Vi) 

Vi) 

M 

o 

V) 


mass  matrix 

aerodynamic  pitching  moment  for  the  ith  test  point,  in. -lb  or 
ft-lb 

indicated  nose-up  pitching  moment  in  rotor  shaft  axis  system, 
gage  units 

chordvise  bending  at  .X  radius,  leading  edge  in  tension,  in. -lb 
upward  flapvise  bending  at  .X  radius,  in. -lb 
fiber  glass  model  chordvise  bending  moment,  in. -Id 

pitch  moment  correction  at  ith  two-dimensional  test  point  for 
normal  force  and  chordvise  deflection,  in. -lb 

full-scale  chordvise  bending  moment,  in. -lb 

replica  model  chordvise  bending  moment,  in. -lb 

2 

effective  mass,  lb-sec  /in. 

2 

torsional  inertia,  in. -lb-sec 
vertical  inertia,  lb-sec  /in. 
fore-aft  inertia,  lb-sec  /in. 
lateral  inertia,  lb-sec '/in. 

fiber  glass  model  flapvise  bending  moment,  in. -lb 
full-scale  flapvise  bending  moment,  in. -lb 
replica  model  flapvise  bending  moment,  in. -lb 
Mach  number  at  ith  two-dimensional  test  point 
generalized  mass  of  the  ith  natural  mode 

pitch  moment  correction  at  ith  tvo-dimensional  test  point  for 
normal  deflection,  in. -lb 

pitch  moment  correction  at  ith  tvo-dimensional  test  point  for 
airfoil  weight  and  normal  deflection,  in. -lb 

static  offset  correction  for  gage  units 

indicated  pitch  moment  at  the  ith  tvo-dimensional  test  point, 
in . -lb 
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*r 

mtfs 

^‘tRM 

^t.xb 

^(i) 


M 


1,90 


BAC 


FD 


FE 


FS 


"(i) 

*PRV 

S 

SFL 

SFM 

SFN 

*SFX 

*SFY 

*SFZ 


N, 


N 


N 


P 

<1 


w 


'vj 


fiber  glass  model  torsion  moment,  in. -Id 

full-scale  torsion  moment,  in. -lb 

replica  model  torsion  moment,  in. -lb 

torsion  moment  at  .X  radius  nose  up,  in. -lb 

pitch  moment  correction  at  ith  two-dimensional  test  point  for 
airfoil  static  unbalance,  in. -lb 

blade  tip  Mach  number  at  90°  azimuth 

static  offset  correction  for  gage  units 

indicated  nose  right  yawing  moment  in  rotor  shaft  axis  system, 
gage  units 

arbitrary  constant  for  discontinuity  detection 

ending  revolution  for  time  history  segment  for  frequency 
analysis 

starting  revolution  for  time  history  segment  for  frequency 
analysis 

indicated  normal  force  at  ith  two-dimensional  test  point,  lb 

digital  time  history  data  values  per  revolution 

rotor  nose  right  yaving  moment  about  shaft,  ft-lb 

scale  factor  for  Lg^,  gage  units/lb 

scale  factor  for  Mg^,  gage  unlts/lb 

scale  factor  for  Ng^,  gage  units/lb 

scale  factor  for  gage  units/lb 

scale  factor  for  Yg^,,  gage  units/lb 

scale  factor  for  Zg^,  gage  units/lb 

2 

static  pressure,  lb/ft 

2 

dynamic  pressure,  lb/ft 
column  matrix  of  modal  amplitudes 
amplitude  of  Jth  chordwiae  mode 

SSI 


1 


Sri 

R 


amplitude  of  ith  flapvise  mode 
rotor  radius,  ft 

physical  equivalent  of  electrical  calibration  signal,  engineer¬ 
ing  units 


G 

J 


RN 

R 


NH 


GS 


in 


out 


r 

vc 

r(wj) 

S 


t 

U 

V 


v*u) 

VD(i> 

VL(i) 

Vi) 


parameter  for  compressible  adiabatic  flow  equations.  Op-  sec^/ft^ 

resultant  frequency  ccnponent  from  numerical  frequency  analy¬ 
sis,  engineering  units 

Reynolds  number 

resultant  of  nonharmonic  amplitude  component,  engineering  units 

radius  of  gyration  of  blade  spar  structural  crossection,  in. 

radial  distance  of  inboard  end  of  blade  segment  from  center  of 
rotation,  in. 

radial  distance  of  outboard  end  of  blade  segment  from  center  of 
rotation,  in. 

vortex  core  size,  ft 

analytical  counterpart  of  Rj 

airfoil  span,  in. 

ratio  of  blade  elastic  stiffness  to  dynamically  scaled  value 
geometric  scale  factor 
kinetic  energy,  in. -lb 

o_ 

wind  tunnel  stagnation  temperature,  F 
time,  sec 

resultant  velocity  with  respect  to  blade,  ft /sec 
flow  velocity,  ft/sec  or  kn 

angle  of  attack  transducer  net  output  at  test  point  (i) 

balance  drag  component  net  output  at  test  point  (i),  gage  units 

balance  lift  component  net  output  at  test  point  (l),  gage  units 

balance  moment  component  net  output  at  test  point  (i),  gage 
units 
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simulated  full-scale  equivalent  speed,  ft/sec  or  kn 

momentum  theory  induced  inflow  velocity,  ft /sec 

digital  data  value,  computer  units 

blade  response  data  value  in  engineering  units 

weight  of  airfoil  for  two-dimensional  test,  lb 

average  data  value  on  tape  in  zero  reference  condition,  com¬ 
puter  un^ts 

average  data  value  on  tape  in  "RCAL"  records,  computer  units 

kth  digital  value  in  selected  time  history  segment  for  fre¬ 
quency  analysis 

average  data  value  on  tape  in  "ZCAL"  records,  computer  units 
blade  weight  per  unit  span,  lb/ in. 

indicated  longitudinal  force  normal  to  rotor  shaft  on  model 
internal  balance,  gage  units 

column  of  modal  deflections  at  head  in  direction  of  force 

Jth  data  value  in  a  time  history  record 

interpolated  quantity  for  Jth  data  value 

static  offset  correction  for  Xg^,,  gage  units 

rotor  longitudinal  force  forward  normal  to  shaft,  lb 

resultant  of  head  response  to  unit  cosine  force  or  moment,  in. 
or  rad 

sine  part  of  head  response  to  unit  cosine  force  or  moment,  in. 
or  rad 

cosine  part  of  head  response  to  unit  cosine  force  or  moment, 
in.  or  rad 

indicated  side  force  to  right  and  normal  to  rotor  shaft  on  model 
internal  balance,  gage  units 

distance  of  mass  center  of  gravity  of  blade  section  aft  of  the 
feathering  axis,  in. 

static  offset  correction  for  gage  units 


JSS 


Ye  rotor  side  force  to  right  normal  to  shaft,  lb 

O 

Zn._  indicated  vertical  force  downward  along  rotor  shaft,  gage  units 

i>AU 

Zq  static  offset  correction  for  Zg^,  gage  units 

Zg  rotor  thrust  force  down  along  shaft,  lb 

a  corrected  angle  of  attack,  deg 

ae  rotor  control  angle  of  attack,  <*8-B2s 

(a^,)  indicated  model  fuselage  nose  up  angle  of  attack,  deg 


°0(i) 


“»(M> 

asmt 


a 

a 


sn 

snt 


(M) 


°T 

0 

Y 


vi 

Vj 


airfoil  midspan  geometric  angle  of  attack  at  the  ith  test 
point,  deg 

start  and  end  angle  of  attack,  deg 

corrected  rotor  shaft  nose  up  angle  of  attack,  deg 

angle  of  attack  for  pitching  moment  stall  at  Mach  number  M,  deg 

angle  of  attack  for  pitching  moment  stall  in  unsteady  aerody¬ 

namic  data  table,  deg 

angle  of  attack  for  normal  force  stall  at  Mach  number  M,  deg 

angle  of  attack  for  normal  force  stall  in  unsteady  aerodynamic 

data  table,  deg 

angle  of  attack  for  unsteady  aerodynamic  data  table  look-up, 

deg 

blade  upward  flapping  at  hinge,  deg 

specific  heat  ratio  parameter  for  compressible  adiabatic  flow 
equations 

deflection  of  the  pth  viscous  damper  in  the  1th  mode 

mode  shape  column  matrix 

1th  flspwise  mode  shape 

slope  of  ith  flapwlse  mode  shape 

curvature  of  1th  flapwlse  mode  shape 

curvature  of  Jth  chordvise  mode  shape 


3S4 


Ar 

Aa 

6F 

r, 

CBCi 

CBFi 

^Ti 

n 

®c 

®e 

6(i) 


X 

X 

\ 

u 

0 

P 

< 

a 

o. 


CFS 


frequency  increment,  cycles/rev 

pressure  transducer  output  for  total  minus  static  pressure, 
gage  units 

radial  length  of  segment  along  blade,  in. 

vail  effect  correction  to  fuselage  angle  of  attack,  deg 

plain  flap  downward  deflection,  deg 

blade  lag  in  drag  direction  at  hinge,  deg 

experimental  st-uctural  damping  for  the  nonrotating  ith  chord- 
vise  natural  mode,  percent  of  critical 

experimental  structural  damping  for  the  nonrotating  ith  flap- 
wise  natural  mode,  percent  of  critical 

experimental  structural  damping  for  the  nonrotating  ith  tor¬ 
sional  natural  mode,  percent  of  critical 

structural  damping  coefficient 

collective  pitch  at  75%  radius 

elastic  blade  torsional  deflection,  rad  or  deg 

airfoil  midspan  torsional  elastic  deflection,  deg 

built-in  blade  linear  twist,  from  axis  of  rotation  to  blade 
tip,  deg 

ratio  of  inflow  velocity  to  ftR 

ratio  of  inflow  velocity  to  flR  from  prescribed  wake  calculation 

rotor  advance  ratio,  V/flH 

air  density  in  flow,  slugs/ft^ 

air  density  at  zero  velocity,  slugs/ft^ 

rotor  solidity,  be /*R 

2 

full-scale  chordvise  bending  stress,  lb/in. 

fiber  glass  mo^el  chordvise  bending  full-scale  equivalent 
stress,  Ib/in.*" 

O 

full-scale  flapvise  bending  stress,  lb/in. 


sss 


°FFS 

°T 

°T.15R 

°TFS 


°TFS.15R 

♦UA 

* 

^SP 

Q 

QS 

“F 

“j 

“n 

“NH 

“ni 


fiber  glass  model  flapvise  bending  full-scale  equivalent 
stress,  lb/in. ^ 

2 

full-scale  torsion  stress,  lb/in. 

full-scale  torsion  stress  at  15?  rotor  radius,  lb/in. c 

fiber  glass  model  torsional  full-scale  equivalent  stress, 
lb/in.2 

fiber  glass  model  torsional  full-scale  equivalent  stress  at 
15?  .otor  radius,  lb/4.r.2 

rotor  head  phase  response  to  unit  vibratory  force,  deg 
azimuth  angle,  deg 
svashplate  phase  angle*  deg 
rotor  rotational  speed,  rad/sec 

simulated  full-scale  rotor  rotational  speed,  rad/sec 

forcing  frequency,  cycles  or  rad/sec 

ratio  of  analysis  frequency  to  rotational  frequency 

natural  frequency,  cycles  or  rad/sec 

nonharmonic  frequency,  cycles/rev 

natural  frequency  of  ith  natural  mode,  cycles  or  rad/sec 
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